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This missile system with nuclea 
capability is on the way to becoming 
a standard Army weapon. To be 
produced by Sperry-Rand, it i: 
scheduled ultimately to replace the 
CORPORAL artillery missile system. 
The Jet Propulsion Laboratory. 
developed SERGEANT is invulner- 
able to any known enemy electronic 
counter-measures. And, thanks to a 
simplified Thiokol solid propellant 
rocket motor, the missile is extremely 
mobile — can be transported on mili- 
tary vehicles. 

When SERGEANT becomes oper- 
ational, the Army will possess one 
of the most maneuverable and accu- 
rate artillery weapons in its history 
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new projects present challenging 
problems and a chance for 
greater responsibility. 


Thiokol 


CHEMICAL CORPORATION 


TRENTON, N. J. * ELKTON, MD. 
HUNTSVILLE, ALA. * MARSHALL, TEXAS 
MOSS POINT, MISS. * BRIGHAM CITY, UTAH 

DENVILLE, N. J. * BRISTOL, PA. 


® Registered trademark of the 
Thiokol Chemical Corporation for its liquid 
polymers, rocket propellants, plasticizers 
and other chemical products. 


a 
q =f & By: 
In 
j Jt 
ag 


High-volume production techniques... perfected by Bulova to maintain their 
leading competitive position in consumer markets... assure repeatability and 
reliability of industrial and defense products. 


To date, Bulova has manufactured more than 4,000,000 arming devices. Although 
reliability was crucial, only approved sampling techniques were required for 
inspection. 


Experience in precision design and precision manufacture is the Bulova tradition, 
the Bulova capability. It has been for over 80 years. For more information write — 


Industrial & Defense Sales, Bulova, 62-10 Woodside Ave., Woodside 77, N.Y. 
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Scope of ARS JOURNAL 


This Journal is devoted to the advance- 
ment of astronautics through the dissemina- 


tion of original papers disclosing new scientific 


knowledge and basic applications of such 
knowledge. The sciences of astronautics are 
understood here to embrace selected aspects 
of jet and rocket propulsion, space flight 
mechanics, high-speed aerodynamics, flight 
guidance, space communications, atmospheric 
and outer space physics, materials and struc- 
tures, human engineering, overall system 
analysis, and possibly certain other scientific 
areas. The selection of papers to be printed 
will be governed by the pertinence of the topic 
to the field of astronautics, by the current or 
probable future significance of the research, 
and by the importance of distributing the in- 
formation to the members of the Society and 


to the profession at large. 
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presentation of the ideas will allow. x- 
clusion of dispensable material and concise- 
ness of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of references, 
1 page of abstract and 12 illustrations. 
Fewer illustrations permit more text, and vice 
versa. Greater length will be acceptable 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes or Technical Comments. 
They may be devoted to new developments 
requiring prompt disclosure or to comments 
on previously published papers. Such manu- 
scripts are usually published within two 
months of the date of receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication, according to the above-stated 
requirements as to subject, scope or length, 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal 
or as a special supplement. if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as a 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be double spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ 
positions and affiliations in a footnote on the 
first page. Equations and symbols may be 
handwritten or typewritten; clarity for the 
printer is essential. Greek letters and unusual 
symbols should be identified in the margin. If 
handwritten, distinguish between capital and 
lower case letters, and indicate subscripts and 
superscripts. References are to be grouped at 
the end of the manuscript and are to be given 
as follows: For journal articles: authors first, 
then title, journal, volume, year, page numbers 
for books: authors first, then title, publisher, 
city, edition and page or chapter numbers. 
Line drawings must be clear and sharp to make 
clear engravings. Use black ink on white 
paper or tracing cloth. Lettering should be 
large enough to be legible after reduction. 
Photographs should be glossy prints, not 
matte or semi-matte. Each illustration must 
have a legend; legends should be listed in 
order on a separate sheet. 

Manuscripts must be accompanied by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests 
with the author. 

Submit manuscripts in duplicate (original 
plus first carbon, with two sets of illustra- 
tions) to the Editor, Martin Summerfield, 
Professor of Aeronautical Engineering, Prince- 
ton University, Princeton, Preprints 
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General I Chemical « offers 
Helpful Technical Data 
on ROCKET FUEL 

OXIDIZERS 


uming 


Liquid Fluorine 


Chiorine Trifluoride and 
other Halogen Fluorides 


To assist specialists working in the field of 
rocket and missile propulsion, General Chemical now offers 
a valuable set of product information bulletins packed with 
physical and chemical data, tables and graphs. This infor- 
mation is based on General Chemical’s extensive experience 
with these products as America’s leading producer of 
fluorine and fluorine-based chemicals, and one of the 
nation’s primary producers of nitric and other mineral 
acids. The questions our customers ask most often about 
these products are answered here. 


Here are the specific bulletins offered: 


“Nitric Acid, Fuming’’—29 pages—Including detailed prod- 
uct description of both red and white fuming nitric acid, 
chemical and physical properties, corrosion data, directions 
for storage and handling. 


‘*Fluorine’’—21 pages—Including extensive data on the chemi- 
cal and physical properties of this most powerful oxidizing 
agent, a description of materials and equipment for han- 
dling gaseous and liquid fluorine, safety precautions. 


*‘Chlorine Trifluoride’’-—35 pages—Including the chemical 
and physical properties of chlorine trifluoride and other 
halogen fluorides, recommended materials for use in halo- 
gen fluoride systems, directions for safe handling. 


Also available are bulletins on: 


Handling Elemental Fluorine Gas In The Laboratory 
Liquid Fluorine Unloading Procedure 
Chlorine Trifluoride Vapor Pressures 


Mail coupon now for your free copy of any or all 6 of these 
vafuable technical bulletins. 
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Bill Knight is a rocket chemist. His field of interest is 
research on advanced solid and liquid propellants to meet 
the stringent requirements of today’s military programs 
(and tomorrow’s). 

The continued success of the American missile program 
depends on men like Bill Knight. 

But creation takes talent. The kind of talent that top- 
notch scientists and engineers can use to the utmost on 


Laight 


Aerojet’s rapidly expanding programs in chemistry and 
other fields. 

There are many projects at Aerojet that intrigue the 
scientific and engineering investigator. Your resume is 
invited. The address: AEROJET-GENERAL CorRP., P.O. Box 
296D, Azusa, California or P.O. Box 1947D, Sacramento, 
California. Attention: Director of Scientific and Engineer- 
ing Personnel. 
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HE name “magnetohydrodynamics” has been given to 

the branch of fluid mechanics that deals with motion of 

electrically conducting fluids in the presence of magnetic or 
electrical fields. 

The name is in some ways a poor one, for one is usually not 
concerned with water but with gases—typically, with gases 
moving at supersonic or hypersonic speeds. Thus, my 
colleague Professor Resler and I have suggested the names 
“magneto-aerodynamics” (1)? and ‘‘magneto-gasdynamics” 
(3). But the name magnetohydrodynamics, sometimes 
abbreviated as MHD, seems to be more popular, and perhaps 
should be accepted, since “hydrodynamics” is, after all, a 
synonym for ‘fluid dynamics” in a broad sense. 

This field is not the same as particle physics, although the 
two have much in common. It is distinguished by the fact 
that the medium in motion is relatively dense, so that the 
differential equations of fluid mechanics describe the flow, 
with new terms being added to them to account for electro- 
magnetic effects. The Maxwell equations of electrodynam- 
ics also enter, but they must be solved together with the 
fluid mechanical equations. From the microscopic point of 
view, one may say that the word magnetohydrodynamics 
applies to cases in which there are many and frequent colli- 
sions between the constituent particles. The equations that 
are solved, like the usual fluid mechanical equations, are state- 
ments about properly averaged quantities. 

Let us now review, from the macroscopic point of view, the 
new phenomena that appear in fluid motions if the fluid is a 
conductor and there are electric and/or magnetic fields pres- 
ent. 

The first and perhaps the most interesting of these phe- 
nomena is a body force equal to j X H per unit volume, where 
j denotes the current density vector and H the magnetic field 
vector. We shall consider briefly the origin of this force later, 
but now let us simply recognize this as a familiar effect, 
namely the force on any conductor carrying a current in a 
magnetic field—the force that makes an electric motor oper- 
ate. Engineers working in fluid mechanics are particularly 
interested in this body force, for, except for gravity, they do 


Presented at the ARS 13th Annual Meeting, New York, 
Nov. 17-21, 1958 

1 This research was supported jointly by the U. S. Air Force 
through the Office of Scientific Research of ARDC and the Me- 
chanics Branch of the Office of Naval Research. Reproduction 
in whole or in part is permitted for any purpose of the United 
States Government. 

2 Numbers in parentheses indicate References at end of paper. 
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not usually have such forces to work with, and even more 
rarely do they have a body force that may be controlled. 

The second electrical phenomenon is the induced electro- 
motive force q < H, where q is the fluid velocity vector. We 
shall postpone detailed discussion of this effect; again it is 
familiar in electrical work, the emf induced by moving any 
conductor across a magnetic field—the phenomenon that 
makes a generator function. In our work it means that 
electrical potentials will be set up, in general, when conducting 
fluids flow in the presence of magnetic fields; thus currents 
will flow, and the body force will be present. 

The third effect to be considered is also familiar in elec- 
tricity, namely Joule heating. It is given by j?/o in energy 
units per unit time per unit volume if o is the conductivity 
of the fluid. o is a material property, the reciprocal of the 
resistance of a unit cube of the fluid medium; thus its units 
are (for example) mhos per em. We shall discuss this fluid 
property later and try to give some idea of its magnitude in 
aeronautical situations. At this point let us simply remark 
that it is the increase of o from extremely small values for 
common gases at room temperature and pressure to appre- 
ciable values under flight conditions that makes a science of 
magneto-aerodynamics possible. Everyone knows that at- 
mospheric air is a very good insulator under ordinary condi- 
tions. On the other hand, an electric arc is a pretty good 
conductor, because the air in it is hot enough to be ionized. 
We are told that the temperatures reac hed in some problems 
of hypersonic aerodynamic heating approach those of electric 
ares; therefore we recognize the possibility that regions of 
appreciable electrical conductivity may exist in the vicinity 
of flying objects. 


Equations of Motion 


The discussion above permits us to write down the equa- 
tions of magnetohydrodynamics; they are: 
Continuity 
dp/dt + div (pq) = 0 = 2am [1] 
Momentum 
Dq 


p= 


in Ae x H) + (viscous terms) [2] 


= 
=, 
a a4) dk J ‘ 
Dt 
907 
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Here D/Dt is the familiar convective derivative, so that 
Equation [3] is a statement of the heat energy added to an 
element of the fluid as it flows. @Q denotes all forms of heat 
added to the element, such as by viscous dissipation, heat 
transfer, combustion or radiation, exclusive of Joule heating. 
Thus a simplifying assumption has been made in this equa- 
tion, namely that there is no other transfer of energy from 
electromagnetic to thermodynamic forms. Although this is 
not correct for fluids whose magnetic permeability is affected 
by temperature or field strength, for example, it is a good ap- 
proximation for our purposes. The viscous terms in Equa- 


tions [2 and 3] have not been written out because they are ~ 


somewhat lengthy and are not actually used in this paper. 

In the absence of the magnetohydrodynamic terms, this 
set of equations is rendered determinate by inclusion of an 
equation of state. Here we employ the equation appropriate 
for perfect gases 


p = ReT [4] 


We shall assume that ® is constant, i.e., that the number of 
particles per unit mass does not vary. There are flight con- 
ditions (where dissociation and ionization of air occur) where, 
to be more precise, ® must be considered a variable. 

We also have a generalized Ohm’s Law, which describes 
macroscopically the flux of electrical charges under the in- 
fluence of an electric field Z and the induced emf mentioned 
above 


j=o0(E+qx [5] 


This formula will be discussed later but also requires some 
explanation at this point. It is a correct approximation only 
for conductors that are electrically neutral; it is a very good 
approximation for ionized gases (as well as for metallic con- 
ductors), because the electrons and positive ions are never 
separated by more than an extremely small distance. The 
vector E includes the contribution due to this separation of 
charges, and this contribution is not necessarily small even 
though the separation distance is minute. One may interpret 
Equation [5] as follows: At any point in the fluid there is 
applied to a small element an emf equal to g X Hi plus that 
part of E that is due to effects exterior to that element, i.e., 
voltage applied by neighboring elements and external elec- 
trodes. This total, local emf is opposed by: (a) the local 
“TR drop” j/o minus (b) the part of E that is due to charge 
separation. The fact that we admit charge separation while 
we describe the fluid as everywhere neutral in charge means 
that we neglect the charge separation distance and consider 
the charges as forming doublets. The fluid may be said to be 
“polarized.” 

It might also be remarked that £ is not invariant to transla- 
tion of coordinates, and must be measured in the same frame 
of reference as the velocity q. 

At this point our set of equations is still indeterminate, but 
we have not used the so-called Maxwell equations. In the 
forms suitable for our problems these are (in electromagnetic 
units and for fluids whose permeability is equal to that of 
empty space) : 


Ampére’s Law = curl H [6] 


Faraday’s Law curl E = —dH/at pee, [7] 


“Displacement currents” have been neglected in Equation 
[6]. Two effects are usually included under this name (see 
(2)), and both of these are negligible in most magnetohydro- 
dynamic problems, as we have pointed out in (3). 


Dimensionless Parameters 


Although the set of equations [1 to 7] is now determinate, 
it is indeed formidable and usually requires some simplifica- 
tion in order to provide solutions, except in cases of particu- 


larly simple geometry. To make such simplifications, one 
resorts to the familiar technique of identifying the significant 
dimensionless parameters and considering the effects of their 
magnitudes. For example, in (1) we have shown that the 
typical ratio of the electrical body force to the dynamic o: 
inertial force is given by the parameter 


N =08 .*L/(1/2)poV [8 


Moreover, we evaluated this number for typical high-spee« 
aeronautical cases and concluded that the ratio is larg: 
enough to produce important effects.* 

The typical ratio of body force to viscous forces is th« 
product of this number N and the Reynolds number, whicl 
is twice the square of the Hartmann number (5). Thus w 
concluded that at high Reynolds numbers it is reasonable t 
neglect viscosity, except in boundary layers and wakes 
as in ordinary fluid mechanics. We shall therefore omit the 
viscous terms in Equations [2 and 3] for the purposes of this 
discussion. There are, in contrast, some important investi- 
gations of boundary layer flow (6 to 8) and channel flow (9 t 
11) with magnetohydrodynamic effects. 

Perhaps the most important dimensionless parameter is the 
so-called “magnetic Reynolds number” 


R, =oVL [9] 


As we have pointed out (1), this can be interpreted as the ratio 
of the motion-induced magnetic field strength to the applied 
magnetic field strength. When R,, is large, the applied field 
may often be neglected, and one is led to the concept of the 
ideal conductor or fluid of infinite conductivity, as we shall 
see. On the other hand, for small R,, one may consider 
the magnetic field to consist of the applied field only slightly 
perturbed. This case will also be treated later. 

Still another dimensionless parameter is found useful, 
namely 


m= (V/H «) 
This is the ratio of the flow speed V tothespeed => 
a=H./V 4p [11] 


which is the speed of propagation of certain magnetohydro- 
dynamic waves—the Alfvén waves. Thus m is analogous to 
a Mach number, i.e., it is the ratio of a characteristic flow 
speed V to a propagation speed of wavelike disturbances. It 
might well be called the “Alfvén number” or the “magnetic 
Mach number.” 

As usual, additional dimensionless parameters can be 
formed by combining those already mentioned. Moreover, 
it is interesting to note that the number N defined in Equa- 


tion [8] is a combination of R,, and m 

N = 12 


so that the Hartmann number is proportional to the square 
root of ReR,,m-?, ete. 

Another interesting combination is the ratio of the mag- 
netic Reynolds number to the Reynolds number, which is, 
of course, equal to the product ov, where v is the kinematic 
viscosity. As we have shown in (15), 47R,, plays the same 
part in defining certain “magnetohydrodynamic boundary 
layers” as does Re in defining the viscous boundary layer. 
Thus the dimensionless number 47rov is an indication of the 
relative thicknesses of these two kinds of layers; actually it is 
the square of the ratio of viscous to magnetohydrodynamic 
boundary layer thicknesses. This is particularly interesting 
because 4rcv is a material property; for example, it has a 
value of about 10~7 for mercury at room temperature. This 


3 Unfortunately we erred in plotting Fig. 2 of (1) where the 
values given are 100 times too small. hus our conclusion re- 
garding the possibility of appreciable body forces might have 
been considerably more emphatic. See (4). 
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parameter might be called the ‘magnetohydrodynamic 


Prandtl number.”’ 


Magneto-Gasdynamic Channel Flow 


In (1 and 3) we considered a problem of particular interest 
in aeronautics: Quasi-one-dimensional flow of a gas in a 
channel in the presence of magnetic and electrical fields. 
\Ve neglected viscosity and assumed all the pertinent param- 
eters, such as cross-section area and applied field strengths, 
to vary only slowly with the distance along the channel z. 
Che equations for this type of flow are (where primes denote 
lifferentiation with respect to x) 


puA = constant [13] 
puu' +p’ —uH)H [14] 


pucpT’ + wu’ = o(E — uH)E +Q 
ogether with Equation [4]. 

These equations cannot be integrated in general, but only 
or special geometries A(x). We can calculate the acceler- 
ition u’, however, in terms of the conditions at any given 
station 2; the result is 


[15] 


(17) 


The first and third right-hand terms in Equation [16] are 
familiar in gasdynamics. The first term describes the effects 
of convergence or divergence of the channel, which are op- 
posite for subsonic and supersonic flows. The third term 
likewise describes the opposite effects of heat addition at | 
subsonic and supersonic flows, namely that heat addition i 
to accelerate subsonic flow but decelerates supersonic flow. 

The remaining right-hand term in Equation [16] is the 
magnetohydrodynamic effect, a combination of the effects 
of body force and Joule heating. For supersonic streams the 
net magnetohydrodynamic effect is acceleration only if u lies 
between the two significant speeds wu; and u; defined in Equa 
tion [17]. This is understandable, for when u > uz the body 
force has changed sign and tends to decelerate the flow, as 
does the Joule heating. « is the speed at which the acceler- 
ation due to body force is canceled by the deceleration due to 
Joule heating, as can easily be confirmed by means of Equa- 
tions [13 to 15]. The situation at subsonic speeds is some- 
what anomalous, for the interval u, < u < uz is the range of 
deceleration, although both the body force and the Joule 
heating would seem to be accelerating the flow. In contrast, 
there is acceleration for u > us, although here the body force 
has changed sign. 

The importance of the Mach number, as well as the fact of 
the singularity at M = 1, leads us to calculate also the deriv- 
ative M’; it may be important, in a magnetohydrodynamic 
channel designed to accelerate (or decelerate) a fluid flow, to 
avoid the sonic regime, because the singularity there clearly 
suggests the phenomenon of “choking.” The result is 


—] 
1+ M? te — Us) (u — ue) — 
2 ap 


Y 2 ap 


JUNE 1959 


CODE: 
dufax > 0 
aM/ax > O 
or A AsympToTE As 


“TUNNEL 


y j CHOKING 


° 


Fig. 1 Magneto-gasdynamic acceleration of one-dimensional 
flow in a constant-area channel. Regimes of positive acceleration 
and increasing Mach number are designated by cross-hatching 
as shown. Arrows denote typical behaviors of flows in the var- 
ious regimes; arrowheads show direction of increasing x 
(downstream ) 


— 


where a denotes 4/yRT and 


= 


2+(y —- 
Since 0 < uz < wu, at subsonic speeds, and w < ue < u3 at super- 
sonic speeds, it is clear that there are regimes where u and M 
have opposite tendencies as the result of magnetohydrody- 
namic effects. 

These tendencies can conveniently be summarized in a 
u, M diagram as shown in Fig. 1. In order to isolate mag- 
netohydrodynamic effects, this diagram is drawn for a con- 
stant-area channel without heat addition. Such a diagram 
contains much useful information, such as the various asymp- 
totic behaviors indicated. The asymptotic approach to the 
point (uw, 0) from subsonic speeds means that the temperature 
increases without limit, as a result of Joule heating. The 
asymptotic approach to us; from supersonic speeds, however. 
is due to the gradual disappearance of the electric current as 
uz is approached. The value uz can be made as large as de- 
sired, 

Of particular interest in Fig. 1 are the two “tunnels” 
through the choking barrier at 1 = 1, one for increasing 
speeds at (wm, 1), the other for decreasing speeds at (u3, 1). 
These are cases for which the numerators as well as the de- 
nominators in Equations [16 and 18] (for A’ = 0 = Q) van- 
ish. Thus these are constant-area magnetohydrodynamic 
analogs of the Laval nozzle and its reverse (diffuser), respec- 
tively. Their stability has not been investigated. 

Equations [13 to 15] with constant o, FE and H have been 
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integrated for constant-area and also for a constant-tempera- 
ture duct by Rosa at Cornell University (13), and for a con- 
stant-density duct by Resler and Sears (3). The constant- 
density case is of course only one of many incompressible- 
flow solutions of Equations [13 and 14], but it is the particu- 
lar one that has gasdynamic significance, since it satisfies 
Equation [15] as well. In (1) we also used the results of 
Equations [16 and 18] to design a constant-area duct for 
maximum acceleration of a supersonic stream, again assum- 
ing constant ¢; this is one along which E/H is adjusted to 
maintain the equality 


(We neglected to point out in (1) that this value produces 
maximum deceleration of the fluid at subsonic speeds.) 


Magnetohydrodynamic Containment 


It has already been remarked that air is a poor conductor at 
the temperature of common structural materials, and only 
becomes a good conductor if considerably higher temperatures 
can be attained. This will pose serious difficulties in the de- 
sign of magnetohydrodynamic channels. There is a possi- 
bility, however, of using a magnetohydrodynamic phenome- 
non to “contain” the fluid, and thereby to keep it away from 
duct walls. This is the famous “pinch effect,” which occurs 
in electric arcs and there confines the ionized gas at high tem- 
perature and pressure. The phenomenon is sketched in Fig. 
2. It isa straightforward result of the body force 7 X H, as 
indicated, and has been analyzed in various places, e.g., (14). 

The simplified model analyzed is one in which the con- 
ducting gas is entirely confined within a cylinder of radius 7, 
i.e., the conductivity for r > r; is zero. Assuming steady 
axisvmmetric flow without viscosity, Equations [2 and 6] 
state, forr <r; 


dr 


= — rH) — [22] 


where j, is the current-density component in the axial direc- 
tion and H¢ the tangential magnetic field component. Ohm’s 
Law (Eq. [5]) tells us that 7, is equal to of, and is there- 
fore a constant proportional to the voltage difference between 
the ends of the cylinder. It is easily verified that the solution 


jz = constant = jo ts 
P = Po — Tr*jo [24] 
He = 2mrjo [25] 
where p» denotes the pressure at the axis of the cylinder. 
For r > 17; 
p = constant = po — mr17j0o? [26] 
He = [27] 


To be “confined” by the magnetic field, then, the pressure 
mr,7jo? must be of the same order of magnitude as the center- 
line pressure pp. This result establishes the magnitude of 
the axial current component for given py and required 1. 


Effects of Varying Magnetic Reynolds Number © 


It has already been mentioned that the parameter R, is 


one of the most important of the dimensionless magnitudes __ 
that determine the character of magnetohydrodynamic flows. 
To illustrate some of its effects, we shall consider in this sec- 


tion a particularly simple class of flows and see how their 
character changes with changing values of R,,. 

Before doing so, it will be useful to simplify our equations, 
Equations {1 to 7], by eliminating both jand E. This is done 
by means of Equations [6 and 7]. Neglecting viscous terms, 
Equation [2] becomes 


bq! - WewH 2 


and Equation [3] becomes 


1 1 2 

Ce all = = — ( curl H) [29] 
Dt Di po \4r 


By taking the curl of both sides of Equation [5], Ohm’s Law, 
and using Equations [6 and 7], one has (noting that div H=0) 


= [—(0H/at) + H-Vq —9:VH — HA divg) [30] 


Thus, our unknown functions are now q, p, p, T and H, and 
the set of equations to be solved consists of Equations [1, 4, 
28, 29 and 30]. 


Small-Perturbation Plane Incompressible Flow: 
Crossed Fields 


Consider now a class of flow problems that constitutes a 


certain extension of the familiar small-perturbation flows of 
classical aerodynamics; namely, suppose that a uniform stream 
is perturbed by a thin cylinder in the presence of a magnetic 
field which, in the undisturbed flow, is also uniform but is 


directed normal to the fluid stream. Thus 

q=(V + uw’, v’, 0) [31] 
and 

H = (hz, Hw + hy, 0) [32] 


where wv’ and v’ are small compared to V, and h, and h, are 
small compared to H.. We assume, and can confirm, that 
h:/H ~ and h,/H ~ are at most of the same order of magni- 
tude as u’/V and v’/V. We consider, in particular, steady 
incompressible flow. If we introduce dimensionless vari- 
ables, where velocity components are divided by V, the pres- 
sure by pV?, magnetic field components by H ., and coordi- 
nate distances by a characteristic length L, we have upon 


linearization 
og -2 oH 
Vi) 


and 


where m and R,, are the quantities defined in Equations [10 
and 9]. 
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These are put into an interesting form by operating with 
the curl operator. Let 2 denote curl g, and let € denote curl 
H (or 41j). Both of these vectors are perpendicular to the 
ry plane. Let the corresponding dimensionless quantities 
Q and = (or 42j) be defined as (L/V)Q and (L/H .)é, re- 
spectively. The results are 


ind 

or oy 4rR, Ox 


We have treated these equations under various assumptions 

egarding the magnitude of #,,._ First consider the case R,, = 
», ie., the case of an ideal conductor. According to Equa- 
tions [37 and 38], both the vorticity Q and the current density 
£/4 satisfy the two-dimensional wave equation. Thus, if 
the flow is undisturbed at y = + © these quantities are con- 
stant along lines having the inclination dy/dzx = m~!. The 
analogy between this geometry and that of familiar (Ackeret) 
small-perturbation two-dimensional supersonic flow is strik- 
ing. On the other hand, the velocity components themselves 
do not satisfy the wave equation, but rather are made up of 
rotational parts distributed in this way and irrotational parts 
distributed according to Laplace’s equation. All of these 
components, as well as the corresponding magnetic field com- 
ponents, are finally determined by the boundary conditions 
and continuity conditions at the solid evlinder. 

The wavelike disturbances running outward and down- 
stream from the wall are recognized as Alfvén’s electromag- 
netic waves (12). It is seen that the disturbance pattern 
consists of irrotational components together with these waves, 
which propagate outward with the speed mentioned in 
Equation [11] and are carried downstream. 

One of the cases treated in detail in (15) is that of an infinite, 
sinusoidal, wavy wall constructed of an insulating substance. 
The calculated pressures at the wall are presented in Fig. 3 in 
the form of a “vector diagram.’ This diagram gives the 
amplitude of the pressure fluctuation and its phase relation- 
ship to the wall geometry for various values of m. The par- 
ticular curve being discussed here is the one marked R,, = o. 

For large but finite values of R,,, the solution just discussed 
can be corrected by the same technique that is used for sound 
Waves propagating in a medium of small viscosity. We 
proceed by assuming that the solutions of Equations [37 and 
38] are now damped sinusoids; for example 


= constant X [39] 


It is then found without difficulty that » =-O(R,.~'), where 
in this case R,, is based on the wave length 27/A. Neglect- 
ing higher order terms, one finds 


= (1/4)m(1 + m?)AR,,~ [40] 


Thus, one effect of introducing small electrical resistance 
in the fluid is to cause attenuation of the Alfvén waves that 
propagate outward from the wall. This is worked out in 
(15), and in (16) McCune has calculated the corresponding 
effect on the wall surface pressure. His result for the sinus- 
oidal wall is shown in Fig. 3 for Rn = 5. In (15 and 16) 
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Fig. 3 Vector diagram showing relative magnitude and phase of 

pressure on sinusoidal wall for various f,, and m. Crossed 

fields. Results for infinite /7,, in comparison with McCune’s ap- 

proximate results for large and small F,,, respectively. fF, is 
based on wall wave length 


these methods were also applied to small-perturbation flows 
around slender cylinders (airfoils). 

Let us now consider, by way of contrast, approximate 
solutions appropriate for small values of R,,. In this case 
we express all the perturbation quantities in power series in 


gs i.e. 
= fo + Rafi + +... 


where f denotes any one of the perturbation quantities and 
terms carrying the subscript 0 are the values for conventional 
flow without magnetohydrodynamic effects. 

Upon introducing such power series expressions into Equa- 


[34], one quickly finds that —, = 4af, or, in dimensional 
form, to the first order in RP, 
(A o/VL)Rmuo = CH [42] 


Equation [33] or Equation [2] is then used to construct fur- 
ther details of the first-order flow field. Thus the first-order 
result corresponds to neglecting the perturbation of the mag- 
netic field in calculating the current density. This is clearly 
consistent with the interpretation given of R,, earlier in this 
paper. 

This theory has also been worked out by McCune, but is 
unpublished as yet.*_ He finds that the successive terms can 
be caleulated by a process of iteration, and moreover that 
the differential equations to be solved at each step of the pro- 
cedure are all linear. His results for the pressure on an in- 
finite sinusoidal wall, using both the linear and quadratic 
terms in R,, are again plotted in Fig. 3, for the particular 
vases R,, = 3. 

It may be of interest to point out that Alfvén waves are 
completely absent from flows with small Rn. 

Finally, for the particular case of the infinite wavy wall, 
Resler and McCune have succeeded, in an investigation still 
unpublished,* in solving the equations for arbitrary Rp. 
Thus, for this case we have results free from any assumptions 
about the size of R,,, with which the approximate results for 
small R,, and large R,, can be compared. For example, in 
Fig. 4 we have plotted the drag per wave length against R,,. 
For the classical case R,, = 0, of course, the pressure pertur- 
bation is in phase with the wall height, and there is no drag. 
The terms proportional to R,, contribute to the drag, how- 


4 These results will be included in ‘The Proceedings of the 


Third Lockheed Symposium on Magnetohydrodynamics,”’ D. 
Bershader, Ed., Stanford University Press, 1959. 
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flows is seen in Equation [30], namely that the expression 


Ry vay within brackets in Equation [30] must vanish. Hayes (17) 
shows that this leads to the conclusion 
| that the number of lines passing any area being 
convened with the fluid is constant. (This is, of course, 
consistent with the Law of Induction, for any change of the 
number of lines intercepted would induce an emf in a circuit 
D > ee the area, and this in turn would lead to infinite 
04 + -_ eurrent flow in an ideal conductor.) We have seen in the 
examples above cases where Equation [49] is satisfied but the 
L field lines are not convected with the fluid. 
In view of Equation [48], one may take the curl of both 
" sides of Equation [45] and conclude immediately that, pro- 
ot vided m + 1 


Fig.4 Coefficient of drag per wave length of sinusoidal wall with 
crossed fields. Resler-McCune theory (solid curve) compared 
with approximate results for small and large R», respectively. 
The drag coefficient is defined as (drag per wave length) + (27? X 
pV*eL~') where the wall shape is given by « cos (27x/L) 


ever, while those proportional to R,.? do not. It is of par- 
ticular interest to note that, at least in this problem, the 
asymptotic value for R,, = © is nearly reached at R,, = 5. 
Tentatively we shall take this to be an indication of what is 
meant by “large” R,,. If this indication is correct, a flow of 
mercury, for example, at large Reynolds number could be a 
case of large magnetic Reynolds number as well, in spite of 
the very small value of the “magnetic Prandtl number” 
4rov mentioned earlier in this paper. i 

i 98 
Small-Perturbation Plane Incompressible Flow: 
Parallel Fields 


In (15) we have applied analogous methods to still another 
class of two-dimensional small-perturbation flows, namely 
those in which the undisturbed magnetic field is uniform and 
parallel to the stream. In this case 


The equations analogous to Equations [33 and 34] are now, 
for incompressible flow 


_ 


The equation for x-components in Equation [45] leads im- 
mediately to the result, familiar in conventional incom- 


pressible ‘small-perturbation cases 
p = constant —pVu’ a ‘ed [47] 


Again let us begin with a discussion of the case R, = ©. 
Equation [46] then leads to 


q-H [48] 


which means that the velocity and magnetic perturbations are 
proportional vectors, and the streamlines and magnetic field 
lines are everywhere parallel to one another. This is an ex- 
ample of a result often stated, that magnetic lines of force are 
convected with the fluid in flows of infinite conductivity. 
Actually this result is not always correct; the general situ- 
ation for ¢ = o for both compressible and incompressible 
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Thus, for infinite R,, the flow is irrotational and current- 
free. The streamline pattern is unaffected by magnetohydro- 
dynamic effects. But the boundary conditions on the mag- 
netic field are not satisfied, in general, by the resulting con- 
figuration without introduction of current sheets at solid 
boundaries. These are sheets of infinite current density, 
and of course they contribute to the force on the boundary. 
Thus, although the streamline pattern and the pressure in the 
flow are not affected by magnetohydrodynamics, the force on 
a solid body is affected. 

In this case the introduction of small resistance (large but 
finite R,,) leads not to the wave equation and wavelike solu- 
tions but to an equation of boundary layer character. This 
equation is Equation [46] (or more generally Eq. [30]), where 
it is seen that the derivative 0?H/05* appears only with the 
small factor R,,~!. It can be proved that there is therefore 
a layer of thickness O(R,,~’/?) in which the current density is 
large and outside of which it becomes negligible. It is clear 
that this magnetic boundary layer replaces the current sheet 
of the infinite-R,, solution. The situation is analogous to 
boundary layer theory for fluids of small viscosity. The 
ideal fluid approximation can be thought of as requiring vor- 
tex sheets along all solid boundaries in order to satisfy the 
boundary conditions of viscous flow. The boundary layer 
then constitutes the next approximation, in which the 
vortex sheet is replaced by a large but finite vorticity, of 
thickness O(Re-"'”*). 

In general our small-perturbation equations are not ade- 
quate for this boundary layer type of calculation, because the 
boundary conditions on H often require large perturbations. 
Equation [30] must then be used to provide the boundary 
layer equation, and it is, of course, nonlinear. The infinite, 
insulating, sinusoidal wall again provides a_ particularly 
simple example, in which the linearized Equation [46] is 
sufficient; this case was worked out in (15). 

For small R,,, on the other hand, the same type of power- 
series solution (Eq. [41]) can be used as for the case of 
crossed fields described above. Once again the appropriate 
form of Ohm’s Law, namely Equation [46], leads immediately 
to the first-order result £; = —47io, or 

J = or 
to first order in R,,. McCune has studied the infinite wavy 
wall under this approximation and finds that positive drag 
appears. Again the first-power terms contribute to drag 
whereas the quadratic terms do not. His results are plotted 
in Fig. 5 in comparison with the value (zero) obtained for 
both large and infinite R,, in (15). Also plotted here are re- 
sults obtained by Resler and McCune for arbitrary R,,. 

It appears that the theory for “large”’ R,, in this case is not 
a good approximation for values of R,, less than several thou- 
sand. The contrast with the results (Fig. 4) for crossed fields 
is striking. We propose, as a conjecture, that the difference 
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is related to the absence of the Alfvén mechanism in the pres- 
ent problem, i.e., that flows producing Alfvénlike waves 
will generally exhibit greater magnetohydrodynamic effects 
than flows, like the present one, in which the boundary layer 


character prevails and magnetohydrodynamic effects are Lap ALL R,, 
carried into the flow by diffusion. Cc —--— SMALL R,, APPROX. 
D 
Small-Perturbation Compressible Plane Flow: i 
Parallel Fields 08 


The analysis of (15) has been extended in still another direc- 
tion by Professor Resler, namely, to plane, steady compres- 
-ible flows of an ideal conductor with parallel undisturbed 
‘elocity and magnetic fields. The extension is relatively 
traightforward and can be sketched here; it is as yet un- 
sublished. 

The equation of momentum, Equation [28], is unchanged, 
recomes Equation [45] upon linearization (where now m? 
lenotes .V?/H and leads once again to a result 
inalogous to Equation [47] 


i) 10 10 10 100 1000 
Rm 


Fig.5 Coefficient of drag per wave length of sinusoidal wall with 
parallel fields. Resler-McCune theory (solid curve) compared 
with approximate results for small 2,,. (Large-F,, approximation 
gives zero drag.) The drag coefficient is defined as in Fig. 4. 
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Pp = Po — PoVu'’ [53] 


3y taking the curl of Equation [45] one again finds a linear 
elation between 2 and & but cannot conclude that either = 
‘anishes 


Q=mE or (54) 


Similarly, Equation [30] still reads, for R, = © 


curl (q X H) = 0 [55] 


which requires g X H to be constant for plane flow, and in 
particular, to be equal to zero in this case. The linearized 


a 8 
form of Equation [55] is L HyPERBOLI 
Vhy = [56] 
Since there is no Joule heating in the ideal conductor, the cf 
flow is isentropic. The linearized form of the z-component 3 oo 
of Equation [30] is v 1 £ 
oh: ou’ Y 
£ 
or, with the linearized equation of continuity ka ~ Sunt . TH 


J a: 
¥ Fig. 6 Diagram showing elliptic and hyperbolic regions for 
Resler’s equation of small-perturbation compressible flow with 
parallel fields 


= 


t- 


In Fig. 6 are shown the regions of an m, M. diagram in 
which Resler’s equation is elliptic and hyperbolic, respec- 
tively. The physical meanings of “elliptic” and “hyper- 
bolic” are related to the existence of finite propagation 
speeds. Thus Equation [60] can be used to determine the 
[58] disturbance-propagation speed; it is V/K where 1 — K? is the | 
coefficient of the first term in Equation [60]; i.e., it is 


by virtue of Equation [53]. After integration of Equation 
[57 |, the following interesting result appears 


Vh, = (1 — ou’ 


The results given in Equations [56 and 58] constitute the 


generalization of Equation [48] for the case of compressible law? + a? — (a%a~?/V?)]'” [61] 


flow. 

Resler has proceeded to introduce the stream function as 
is sometimes done for small perturbation flow (see (18), for 
example) 


oy 
oy 


and has worked out the differential equation for y. He finds 


(1 M )(m 1) [60] 

Mo*+m?—1 dy? 
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where a is the Alfvén speed defined in Equation [11]. This 
is in agreement with an analysis of the speeds of propagation 
of plane magnetohydrodynamie waves. Thus K rather than 
M « is the ratio of the stream speed to the true sound speed, 
given by Equation [61]. When Equation [60] is hyperbolic 
(K > 1), the wave angle of disturbances is sin~! (1/K). 
It is interesting to note that there is a hyperbolic region at 
subsonic flow speeds (see Fig. 6). In this region the waves 
(characteristics) involved are forward-facing; i.e., they slant 
upstream from the point of disturbance. This follows from 
the fact that the wave angle increases with flow speed V. 
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In the elliptic regions, Equation [60] is transformed into 
Laplace’s equation by an affine (‘‘Prandtl-Glauert”’) trans- 
formation. Thus the rotational compressible magnetohydro- 
dynamic flows are related to certain irrotational incom- 
pressible flows. 

Solutions of Equation [60] for hyperbolic regions can be 
written immediately 


=fic —- VK?-1y) + 92+ VK?—1y) [62] 


where f and g are arbitrary functions. This form can easily 
be applied to determine the flow pattern and pressure dis- 
tribution in two-dimensional flows. For a flat plate at angle 
of attack e, for example, at both the upper and lower surfaces 


v' = —Ve [63] 
and 
— 1 lower 


Thus, according to Equation [53], the pressure difference be- 
tween lower and upper surfaces is 
Ap = 2p2V?[~WK? — 1/(Mq? — Ile [65] 


Just as in the incompressible case treated in (15), this solu- 
tion, together with Equations [56 and 58], involves a surface 
current at the plate. Its value J, is given by —1/4z times 


(66) 


An upward force equal to H ..J, is exerted on the plate by vir- 
tue of this current. Thus the total upward load per unit 
area of the plate is 


2p K? — 1)(1 — [67] 


In the incompressible limit, K? — 1 becomes 1, and the result 
given in Equation [67] is reduced to that obtained in (15). 


Some of the examples treated in the preceding sections sug- 
gest aeronautical applications of magnetohydrodynamics. 
There appear to be possibilities of using the phenomena to 
affect the lift and drag of flying bodies, to inhibit the transfer 
of heat from hot gas streams to adjacent bodies, and to pro- 
duce or augment thrust. It is also clear that the magneto- 
hydrodynamic accelerator, with electric potential reversed, 
becomes a generator of electric power; thus there appears 
to be a possibility of using these effects to extract electric 
power directly from a flowing, conducting fluid. 

Although it is still too early to predict which, if any, of these 
possibilities will be found after engineering research and de- 
velopment to be practical and economical, some of the engi- 
neering problems that will require solution to make magneto- 
aerodynamics a reality are easily recognized. The most 
obvious are those related to the great weights of electrical 
power components. It will be necessary to provide both 
electric and magnetic fields at much less cost in weight than 
in conventional devices. 

It appears that the aerodynamicist will also be faced, for 
some time at least, with problems of producing and main- 
taining electrical conductivity in gases. In (1) we presented 
some typical data on the conductivities of pure air and of air 
“seeded” with a trace of an easily ionized substance, namely 
potassium. Although those data were based on the best 
available information at the time, recent shock tube meas- 
urements by Lin (19) show considerably greater values of 
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Fig. 7 Electrical conductivity of air measured in shock tube at 
initial pressure 1 mm Hg (ref. 19), and calculated conductivity of 
air seeded with 0.1 per cent potassium 


o. A new graph based on Lin’s measurements in air and a 
corresponding calculation for seeded air is shown in Fig. 7. 
As pointed out in (1), these values of o do not make air a 
very good conductor of electricity, in comparison with metals, 
for example. 

There are, moreover, some phenomena related to the con- 

Bene of electricity in gases which cause serious loss of con- 
ductivity at low densities. Fundamentally, the conductivity 

_ arises from the forces that act on charged particles in an elec- 
tric field. Such forces act in opposite directions on electrons 
and ions; the electrons, being more mobile, account for most 
of the conductivity. In a dense gas they experience collisions 
with atoms and ions, but an electron drift is set up which 
constitutes a current. 

In empty space, however, a charged particle under the in- 
fluence of simultaneous electric and magnetic fields does not 
travel in the direction of the electric field but tends to circle 
about the magnetic lines of force (14), since it experiences a 
force proportional to the vector product of its velocity and 
the local magnetic vector. 

The electrons and ions of an ionized gas would perform this 
same circling motion if it were not for their collisions with 
other particles. It is, in fact, the anisotropic nature of this 
motion, that, through collisions, gives rise to the body force 
j X H introduced in Equation [2] and mentioned frequently 
in this paper. Now, in a rarefied gas the circling motion be- 
comes more important because of the greater time between 
collisions, and the drift of electrons in the direction of the elec- 
tric field is thereby inhibited. If the electrons can flow in the 
direction perpendicular to both fields, they do so, and they 
produce the Hall current. The theory of this situation is 
worked out in (14) and elsewhere. The results show how the 
electrical conductivity is affected by density and magnetic 
field strength; for the strong fields involved in magnetohydro- 
dynamics, the conductivity would effectively disappear at 
altitudes of about 200,000 ft. 
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Fortunately, the conductivity is partially restored if no 
Hall current can flow; instead, a Hall potential is set up, and 
the electrons again drift in the direction of the original field. 
Thus, there are possibilities of maintaining conductivity at 
high altitudes, but their successful realization must depend 
on details of the geometry. The resulting boundary value 
problems would seem to be very complicated in practical 
cases. Thus the kind of magnetohydrodynamic calculations 
described in the present paper may really be greatly over- 
simplified in comparison with the practical problems that 
will occur in connection with the design of vehicles for flight 
at very high altitudes. 

In conclusion, it can be reported that the field of magneto- 
ydrodynamics in the United States is attracting the atten- 
tion of a growing number of engineers who specialize in fluid 
nechanics. Both theoretical and experimental studies are 
veing undertaken, although there is little indication of how 
1 when this intriguing new branch of fluid mechanics will 
‘ontribute to practical engineering. Successful applications 
may have to wait for new developments in other fields, such 
as lightweight fission or fusion power sources. It seems more 
profitable to continue with basic studies of the field, so that a 
store of knowledge and experience will be built up for future 
use. The history of science seems to show that new phe- 
nomena ultimately find practical applications, and this should 
ipply to the phenomena of conductivity and electrical body 


force in high-speed aerodynamics. 


Nomenclature 


a = V/ yRT 
A = cross-section area of duct 
Cy = specific heat, constant pressure 


Cy = specific heat, constant volume 


E = electric field vector (magnitude E) 

f = perturbation quantity (Eq. |41]) 

fig = arbitrary functions (Eq. [62]) 

hr, hy = perturbation components of magnetic field vector 

H = magnetic field vector (magnitude H) = =< 

He = tangential component of H et 

J = current-density vector 

Sy = components of j 

Jo = see Equation [23] 

SF = surface current (current per unit length) 

K = ratio of V to speed of propagation of i a, 
Equation [61] 

L = typical length 

m = V/a 

M = 

N = dimensionless quantity defined in Equation [8] 

p = pressure 

q = fluid velocity vector (magnitude q) 

Q = heat added in addition to Joule heating 

= radial coordinate (cylindrical ) 

r = value of r at outside edge of conducting cylindrical 
stream 

R = gas constant 

Re = Reynolds number VL/y 

Rm = magnetic Reynolds number VLoe 

t = time 

u,v,w = components of g 

v0" = perturbation velocity components 

U1, Us, Uy = significant speeds defined in Equations [17 and 19] 

V = typical speed, e.g., stream speed 

= Cartesian coordinates 

a = speed of Alfvén waves (Iq. [11]) Are 

Y = Cp/Co 

€ = angle of attack se 

d = 2m X (wave length of sinusoidal wall) ~! 

“ = damping coefficient of waves i. 

= kinematic viscosity 

= curl H (magnitude ¢) 

p = density of fluid 

o = electrical conductivity of fluid 

y = stream function 


vorticity (magnitude 


Subscripts 
co = undisturbed, free-stream, or typical value ; 
0, 1, 2,... = successive approximations as defined in Equation 


[41], except as otherwise noted above 
Bar indicates dimensionless quantities, and V is the operator in- 
volving derivatives with respect to %, 9 and 2. 
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Rocket Thrust Transients 


Chance Vought Aircraft, Inc. 


Dallas, Texas 


Rocket thrust termination transients are of importance in the range control of ballistic missiles 
_ and in the separation of boosters or different missile stages. Items of interest include the nature of 
the thrust decay, the impulse increment of the thrust termination transient and the time constant 
of the transient. Equations are developed for these items for idealized cases with reasonable sim- 
plifying assumptions, and the results plotted in generalized form. This analysis is considered to be 
a good approximation for liquid propellant rockets with instantaneous propellant shutoff and 
“‘sliverless”” burnout of solid propellant rockets. 


Assumptions Gas Flow Through the Nozzle ag a 
N A ROCKET combustion chamber, the rate of burned 
gas generation is equal to the sum of the gas flow rate out The gas flow rate through a choked nozzle during a thrust 
of the combustion chamber through the nozzle plus the rate termination transient may be expressed by the following 
at which gas accumulates in the combustion chamber equation 
my = tie + (d/dt) (1) k 
The following assumptions are made: gRT..\k+1 : 

1 The transient is initiated by instantaneously i 
The mass of gas in the combustion chamber at any instant 

2 The gas properties are uniform throughout the com- ad ¥ 
bustion chamber at every instant of time. 

3 The gas flow through the nozzle may be considered as gRT [3] 
a quasi-steady flow process. 

4 The nozzle is choked throughout the transient. This is Differentiating Equation [3] with respect to time and sub- 
true for expansion into a vacuum, and is true for the most stituting the result and Equation [2] into Equation [1] for 
significant part of the transient for finite values of ambient the case when mm, = 
pressure. 

5 The gas processes in the combustion chamber and nozzle _L* kh 1\GFDAE-1 
are isentropic during the transient. ) 

6 The perfect gas law applies. g RT, » dp. 
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Thrust decay transient for Case 1 


Fig. 1 


pressure during the transient. Integrating Equation [4] ” 
between the limits 4; and ¢, and p,, and 


Thrust Decay During Thrust Termination 
Case Fixed Nozzle Area Ratio and Zero Ambient 
Pressure 


This case is the more practical of those considered because 
most rocket nozzles have fixed geometry, and the ambient 


pressure is very small for the altitudes at which thrust termi- _ 
nation is initiated for range control and missile stage sepa- 


The rocket thrust equation is (1)? 

F = pAiCr (6] 
It may be shown that Cy is constant during the transient 
under the assumed conditions. 

The ratio of thrust at any time (6é) after initiation of 
thrust termination to the initial thrust may be determined by 
combining Equations [5, 6 and A-1] for the time (N6t) set 
equal to (t — ¢,) 


F | —2k/(k—1) 
—) =(14+———N 7 
( 2 ) 


This equation is plotted in Fig. 1 for k = 1.20 and k = 1.25 


ration. 


Case 2 Variable Nozzle Area Ratio and Finite 
Ambient Pressure 

This is an idealized case in which the nozzle area ratio is 
assumed to be continuously variable during the thrust tran- 
sient so as always to provide complete expansion to ambient 


pressure. The nozzle throat area remains constant. The 
thrust equation for this case is (2) 
(k+1)/(k-1) (k—-1)/k 
F = 2k? ( 1 
k-—1/\k+1 Pe 
[8] 


2 Numbers in parentheses indicate References at end of paper. 
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_ Fig. 3 Thrust decay transient for Case 2 and k = 


1.25 


The ratio of thrust at any time (N6t) after initiation of 
thrust termination to the initial thrust may be determined 
in the same manner as for Equation [7] 


=a 2 y) x 
1 


(& — + — 1)/2)}2N? 


[9] 


For zero ambient pressure, Equation [9] reduces to Equation 
[7]. Equation [9] is plotted in Figs. 2 and 3 for k = 1.20 and 
k = 1.25, respectively. 
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Impulse Increment During Thrust Termination 


Case l Fixed Nozzle Area Ratio and Zero Ambient Pressure 


one The impulse during the thrust termination transient is 
Al = Fat [10] 


Substituting Equations [4 and 6] into Equation [10] and 
integrating 


Pe L*A 2 k + 
VkgRT.. \k +1 2 


A reference impulse increment for this case may be obtained 
from Equations [6, A-2 and A-3]. Combining the result with 
Equation [11] 


process during the transient relative to 
that during steady state operation with the same quantity of 
gas as discussed in Appendix 1. 


Case 2. Variable Nozzle Area Ratio and Finite 
Ambient Pressure 


The impulse during the thrust termination transient in this 
case is the following when written in terms of a and z from 
Equations [4 and 8] 


Al 
V/ 2ax z? 
dz 


If the exponent (2k — 1)/(k — 1) in the above equation is 
an integer, the integral can be integrated by means of Equa- 
tion [224] of (3). This exponent is an integer for both k = 
1.20 and k = 1.25, which are reasonable values for rocket 
exhaust gases 


(Al) ¢=1.20 = 6 ( 


Pa L*A t 


yX 
V/ 1.20 gR =) V2a 


8 8 
1 
* 


231a3 
(ADg=125 = 10 V2 (2a — x 
V1.25gRT.) N2a 


= 


9a 21a? 105a3 


2 
{14b] 


The reference impulse increment for this case is obtained 
in the same manner as for Equation [12]. Then equations 
analagous to Equation [12] may be obtained 


(3) = §(2a i) 
81 
[i+ 
lla’ 99a? 


Ol /k=1.25 :105a* 31504 


[15b] 


+ [15a] 


* 1155a‘ * 3465a° 


For complete expansion into a vacuum, Equations [15a and 
15b] reduce to Equation [12] for k = 1.20 and k = 1.25, 
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Fig. 4 Thrust termination impulse increment for Case 2 


respectively. Equations [15a and 15b] are plotted in Fig. 
4 for various values of initial pressure ratio. 


Time Constant of Thrust Termination Transient 


The following approximation for the time constant of the 
thrust termination transient may be made based on the 
suggestion in Appendix 2 and Equation [B-2] 


At = Al/F; [16] 


The points for At indicated on Figs. 1, 2 and 3 show that 
At is a good approximation to the time constant for the thrust 
termination transient. An exponential curve corresponding 
to Equation [B-1] is also plotted for reference, indicating that 
the thrust decay curves for the higher initial pressure ratios 
do not differ appreciably from an exponential function. 


_Pa k = 1 Appendix 1 


Reference Time and Impulse Increments 


The selected reference time may be defined as the average 
“dwell time” of the gases in the combustion chamber during 
steady-state rocket operation as defined by Equation [1] 
when (d/dt) (m,) = 0. The following equ ition may be 


written for this ‘‘dwell time” 


Substituting Equations [2 and 3] modified for steady-state 
operation into Equation [A-1 ] 


L* + wat 
VkgkT. N 


By definition, the mass of the total gas flow through the 
nozzle during the time 6¢ is equal to the mass of the gas con- 
tained within the combustion chamber at any instant during 
steady-state conditions. This is also equal to the mass of the 
gas contained in the combustion chamber at the beginning of 
the thrust transient. 

The reference impulse increment is defined as 


= [A-3 ] 
where F is the steady-state thrust which is equal to the thrust 
at the beginning of the transient. re ee 


In a system which has an exponential output described by 
thefunction 


in response to a negative step input, on Shei eee’ constant is de- 


fined as that time at which the output has decreased to the 
fraction (1/e) or 0.368 of the initial value and occurs when 
t = 1/c (4). 

The area under the above exponential curve may be ob- 
tained by integrating Equation [B-1]. This is equal to the 


Appendix 2 Response to a Step Input 
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area of a rectangle with the height equal to the initial function 
value (6) and the width equal to the response time (1/c). 
This suggests the following means of defining the time con- 
stant fer this application 


Nomenclature 
A = area 
= (3) 


a 

b = constant 

(, = thrust coefficient 
c = constant 

e = base of natural logarithms 

= thrust 

¢ = standard acceleration of gravity 
4/7 = thrust termination impulse increment 
6’ = reference impulse increment 

k = ratio of specific heats 

L* = eombustienc hamber characteristic length 
nm = mass 


indicated. 


NALYSES of chemical rocket systems have shown that 

performance limitations, aside from material temper- 
ature limitations, result from molecular characteristics, such 
as bond energies, molecular weight and molecular dissociation 
energy. 
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Effect of Dissociation on the 
Performance of Working Fluids for 


Recent developments indicate a growing interest in nuclear propulsion. 
enthalpy relations are derived for an idealized diatomic gas to investigate the effect of molecular 
weight, dissociation energy and pressure on the performance of diatomic working fluids. It is 
shown that maximum performance for a given molecular weight should be obtainable from mole- 
cules with bond energies in the range 1.5 to 2.5 ev. 
The results of the study for the diatomic system also indicate desirable ranges of 
- molecular weight and aon energy for more complex molecular systems. 


m = mass rate with time 

N = number of reference time increments 
p = pressure 

R = gas constant 

= time 
At = time constant a 
= reference time increment 

= 

@ = function 


Subscripts 

a = ambient 

b = burned gas 
c = combustion chamber 
t = nozzle throat 


= initial conditions 
lal conditio 
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In this paper, general 


The influence of pressure on dissociation is also 


oo 


With increasing interest in nuclear propulsion, it is desir- 
able to obtain corresponding limits in performance for stor- 
able fluids which could be employed for nuclear engines. 
It would be expected that the same factors limiting perform- 
ance in chemical systems would also limit performance in 
nuclear systems, namely: Bond energy, molecular weight 
and molecular dissociation energy. 

In a nuclear propulsion system, there is essentially a con- 
stant-temperature, infinite energy source. Consequently, 
working fluids should be selected which will absorb the maxi- 
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mum energy. Such systems could include complex poly- 
atomic molecular fluids with large numbers of degrees of 
freedom which are not excited at low temperatures, and sys- 
tems with low values of dissociation energies and bond 
energies. At the temperatures considered for practical 
heat transfer processes, the effects of ionization may be 
neglected since relatively little contribution to enthalpy can 
be expected. 

A second requirement is that the excitation process be re- 
versible in the sense that stable, excited or dissociated (free 
radical) states exist in thermodynamic equilibrium at the heat 
transfer temperature in order that this energy can be returned 
to the system by recombination during the expansion process. 

Complex molecules with large numbers of degrees of free- 
dom would tend to have large molecular weight values and, 
hence, small enthalpy per unit mass values. The most 
promising storable fluids for nuclear applications would, 
therefore, appear to be those with low molecular weights, low 
bond energies and low dissociation energies. 

In this paper, a general analysis is undertaken of the rela- 
tionship between dissociation energy, molecular weight and 
total enthalpy. Enthalpy values are determined for a typical 
homonuclear diatomic fluid based upon the characteristics 
possessed by several common diatomic liquids with low vapor- 
ization temperatures. The enthalpy values are computed 
in the molecular weight range 2 to 80 and the dissociation 
energy range | to dev. 

The results of this study are generalized for more compli- 
cated molecular systems. Detailed calculations for poly- 
atomic systems will be presented in a subsequent paper. 


Assumptions and Approximations 


In the enthalpy calculations, a temperature of 6200 R was 
assumed as a characteristic maximum operating temperature 
for heat transfer applications. C nes ulations were carried 


out for the molecular process 


in the molecular weight range 2 to 80 and the dissociation 
energy range 1 to Sev. 

The equilibrium constant and enthalpy of the individual 
constituents are dependent upon the molecular and atomic 
parameters for the ground and excited states. At a temper- 
ature of 6200 R, only the ground state of the atomic species 
will contribute significantly to the enthalpy. A general 
atomic state 7P was assumed in Table 1. 

The enthalpy contributions of individual atomic and 
molecular components are independent of the statistical 
weight. Statistical weight factors will only effect the equilib- 
rium constant, and even this effect is not significant for large 
values of the equilibrium constant (K > 1). 

At a temperature of 6200 R, the summation over energy 


o,(1 — e~™) 


with a statistical weight factor of 1. Calculations indicated 
that a good approximation to the molecular enthalpy con- 
tribution would be given by 


with 
7 
o, = 26 X 10-4 


un = 0.24 


These values were based upon averages for the molecula- 
parameters of typical diatomic molecules. 

The molecular parameters will significantly affect the cal- 
culation of the degree of dissociation only when the equilib- 
rium constant is of the order of 1. In this paper, primary 
interest is centered on systems which are almost completel) 
dissociated. Consequently, the molecular mole fraction wil! 
be approximately zero, and the mole cular contribution t 
total enthalpy will be small. a 7 


Enthalpy Relations 


For thermodynamic equilibrium, the enthalpy per mole is 
given by (1,2) 


H (x, 


1 
‘a 


RT M 


1 df As *) 
dT 


where X; is the mole fraction and 
afi 
aT 


is the dimensionless enthalpy contribution of the species 7. 
fi is the partition function. Using the molecular and atomic 
approximations, the enthalpy per unit mass becomes 


RT? Da, 3 
— = 3.40 X4. + 2) 
[ x. + + 3.40 Xa, {2} 


The equilibrium constant may be written in the form 


Ks, = 2.11 X 10* po/pM (3) 


Mole fractions can be calculated from the relations 


X4 


1 
(3/4) + (1/4) + 16/Ka 


1— Xa (5) 


Ionization contributions were neglected. Although the 


ionization process will contribute to the enthalpy in the form 


states 


The molecular state was approximated by the state 'Z 


may be approximated by 


Table 1 Atomic statistical weight factors 


Atom State Statistical weight 
Cl 2p 6 
H 28 2 
F 2p 6 
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the mole fraction, X4+, will be extremely small. For ex- 
ample, using an ionization potential of 4 ev, the equilibrium 
constant becomes 

p/poK = 1.913 X 10-5 


Typical ionization values are shown in Table 2. 
Typical dissociation energy values are shown in Table 3. 


Results 


As shown in Fig. 1, the degree of dissociation for the di- 
atomic system is strongly dependent upon the dissociation 
energy of the molecule as well as the density or pressure. 


* Numbers in parentheses indicate References at end of paper. 
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of 10 atm, peak enthalpy values occur at dissociation energy _ 
Table 2 Atomic ionization potentials values between | and 2.5 ev. 
ed ' At pressures between 10 and 100 atm, diatomic molecules 
leila ential potential sun with a molecular weight of 6 and a dissociation energy of 
; —— - . 2 ev approach the enthalpy of hydrogen molecules with 
Cs a weight of 2 and a dissociation energy of 4.5 ev. 
a 43 7 The calculations for a molecular weight of 2 are in reason- 
we) 11.264 ably good agreement with the actual enthalpy values for 
in) 13.595 hydrogen at high pressures. At lower pressure values when 
a, er a 13.614 al the equilibrium constant and fraction of dissociation are more 
N 14.540 dependent upon molecular parameters, poorer agreement is 
obtained as a result of the use of average molecular parameters 
typical of heavier molecules. 
Fig. 3 indicates the effect of different heat transfer temper- 
atures on enthalpy for a gas of molecular weight 10 and pres- 
Table 3 Molecular dissociation energies (D)) aN sures of 1, 10 and 100 atm. 
7 Molecule Dev) Mol. wt. = The square root of enthalpy values shown in F igs. 2 and 3 
; Lis 03 furnish a direct indication of the relative merits of various 
IC] ail, 9-452 diatomic gases for propulsion applications. 
ss BH té‘iéi kD The specific effect of pressure on the degree of dissociation 
Ch  —_—so2.«. AT enters only through the mole fraction and molecular weight 
LH  — 2.8 terms. In enthalpy calculations, the enthalpy contribution 
CIF 2. 616 of the individual species is independent of pressure. When 
Bz ya 3.6 no dissociation occurs, therefore, the enthalpy is pressure 
3.75 independent. Specific pressure dependence of the enthalpy 
. aoe _ £.476 and specific impulse can be demonstrated for small values of 


At high densities or pressures (for density ratios greater than 
1, P/Po > 10) practically no dissociation will occur for dis- a. . TEMP. =3444.4°K 
sociation energies greater than 4 ev unless the molecular weight 01 


approaches 100. ———— P/P,= 10 
As shown in Equation [2], the enthalpy per unit mas Wie P/P.2100 
varies directly with dissociation energy. The variation of 
the square root of enthalpy with dissociation energy for pres- - 220 fs 
sures of 0.1, 10 and 100 atm is shown in Fig. 2. The square : 
root values plotted permit a direct comparison of relative 
specific impulse values. 200K 
Peak enthalpy values, at pressures of {9 atm, occur at _ 
about 3 ev for any given molecular weight. Above pressures 
— 
MOL.WT. =10 = 
- ° 
TEMP. = 3444.4°K 
TEMP. = 3100°K —— 
TEMP. = 2800°K ---- 
08 z H5(P/P, = 100) 
WwW WwW 
oO . 
S 
oO O6F 
° 
= = sot 
4 
< 
= 
re) 
< 
2 3 4 5 2 3 4 5 
DISSOCIATION ENERGY (ev.) DISSOCIATION ENERGY (ev.) 
Fig. 1 Effect of dissociation energy and density on the degree of Fig. 2 Variation of enthalpy with molecular weight and dissocia- 
dissociation for diatomic molecules. A:72A tion energy for diatomic molecules. A,=2A 
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The dissociation contribution to enthalpy can be much 
larger than the molecular contribution. It can be expected 
that as dissociation occurs, the specific impulse will demon- 
strate a strong dependence on pressure of the form 


in regions where dissociation occurs. In regions of no dissoci- 
ation or almost complete dissociation, enthalpy and impulse 
are pressure independent (Fig. 4). As shown in Fig. 5, at 
high pressures specific impulse values between 800 and 1000 
sec should be achievable for diatomic systems having molecu- 
lar weights between 6 and 10 and dissociation energies of 
above 2 ev. 

The relative effects of dissociation energy, pressure and heat 
transfer temperature on peak specific impulse values are 
shown in Fig. 6 for a molecular weight of 10. 

The results of these analyses can be generalized for more 
complic ated molecular systems which possess bond energies 
in the range 1.5 to 2.5 ev and whose diatomic constituents 
also possess dissociation energies in this range. For these 
values, the molecular system upon breakdown will be com- 
posed essentially of free radicals which are stable at the high 
temperatures. The degree of free radical formation will de- 
termine to a large extent the advantages or disadvantages 
of storable propellants for high pressure applications. De- 
tailed calculations for complex molecular systems will be 
_ presented in a subsequent paper. 


‘dite 


g 


SQUARE ROOT OF ENTHALPY 


Summary 


This analysis is applicable only to the dissociation of dia- 
tomic systems wherein essentially free radicals are formed as a 
result of the molecular breakdown. It would, therefore, 
apply to any solid or liquid substance with a low dissociation 


DISSOCIATION ENERGY (ev.) 
Fig. 3 Variation of enthalpy with dissociation energy and 


temperature. A:—2A 
| 
dissociation. ‘The enthalpy may be written in the form ais P/P, =10 
cal CAT) 
gm VP 
where the first term contains the enthalpy contribution of the Ho (P/ 10) 
molecules, and the second term the contribution of the atoms. 4 
In expansion through a nozzle, the change in enthalpy be- as 2 
tween the inlet and exit temperatures (i and f) can be 2 1200 F 
written approximately as 
4 
an ~ ery - CT's) 5 1000+ 
= 
and consequently, the impulse in a dissociating system will " 
vary as 800F 
4 
x 100 fe) 
90 
E 80 200 F 
70 
«zo 
1 
8 ‘ol 102 DISSOCIATION ENERGY ( ev.) 
iP Pp Fig. 5 Variation of optimum specific impulse values with dis- 
4 ‘Vastation of enthalpy with sociation 
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energy. The general calculation may be extended to more | 


3.0 
i complicated molecular systems provided that in the process 
: e of molecular breakdown, fragments are formed which are 
2.8F stable in the high temperature region. The energy of dissoci-— 
a 1H ation which will increase the enthalpy at a given temperature _ 
were Fa : will then be returned to the system upon expansion. 
2.6F fore If, in the breakdown of the complicated molecular system, | 
a niet stable molecular products are formed at the high temper-— 
ature which will persist in the expansion process at low tem-_ 
peratures, performance gains may not be achievable. How- 
j4 it is expected that in the general breakdown of molecular 
lite systems free radicals will be formed. The extent of the free 
oT ve iicieuabas radical formation will determine to a large extent the ad- 
Zz are vantages or disadvantages of working fluids at high pressure 
levels. 
<> Probably the most interesting conclusion from these in- 
5 vestigations is that it may be possible to achieve a storable 
4 =10 atm dill fluid in the pressure range 10 to 1000 psi which will 
i *100 atm closely approach the performance of hydrogen at high pressure 
s ‘7a ~ ae Such systems would be composed of the lightweight 
rel atoms, such as carbon, hydrogen, oxygen, lithium, boron and 
a T23444.4°K 
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e ‘Use of Drag Modulation to Reduce 
Loads Dur 
Sin Atmospheric Entry CLARENCE B. COHEN? 


Space Technology Laboratories, Inc. 
Los Angeles, Calif. 


oe ‘a Paarup bitin The design of a space vehicle which must be capable of entering the Earth’s atmosphere involves, 
: it among other things, the consideration of the accuracy with which its landing point can be pre- 
a or dicted, the maximum deceleration loads which the vehicle will experience, and the aerodynamic 
a a a heating to which the vehicle is exposed. Since minimization of the landing point dispersion may 
require the use of entry angles sufficiently large to cause subsequent excessive deceleration loads to 
a the vehicle and its occupants, a scheme is examined whereby these loads may be reduced to a tolera- 
ble level. This is accomplished by variation of the drag of the vehicle (either discretely or con- 
tinuously) in a properly programmed fashion. It is found that the use of continuous drag modula- 
tion can reduce the deceleration loads by as much as 50 per cent, without significantly affecting the 
total aerodynamic heating of the vehicle. Particular numerical solutions of the complete equa- 
tions of motion are also included. These substantiate the results of the approximate analysis. In 
addition, several specific applications of drag modulation have been considered, such as to a ve- 
hicle returning from the moon. 


MONG the problems of recovering a vehicle from a geo- dynamic heating and — point prediction problems will 
centric orbit are: Minimization of the total heat input be considered briefly. 
to the body, which in turn minimizes the weight of the heat j ' 
shield necessary to protect the payload; increasing the ac- E ; f Moti 
curacy with which the landing point can be predicted; and 
control of the vehicle trajectory in such a way that the peak The problem at hand will be aia by using - the two- 
deceleration load experienced by the vehicle is held to a tol- dimensional equations of motion of a point mass being acted 
erable maximum. Although these and other factors are not upon by gravity, inertia and aerodynamic forces. These 
independent in their effect on the design of the vehicle or its equations are Meutviasl iin (1)3 thik etinen aA round, non- 
eae , » the latter subject herein treated in detail as a rotating Earth, a variable gravitational acceleration g, and 
necessary step in a system analysis. . that the drag force acts in a direction parallel and opposite to 
A method is suggested whereby the largest deceleration the inertial velocity vector.‘ 
load occurring during entry into the atmosphere can be held 
to a selected value tolerable by human beings or other deli- ‘7 . 
cate payloads. Essentially, this method involves a varia-  -V db L\ pV? i 
tion of the vehicle parameter W/C>pA in such a way that, for ag at ie (5) 2A + cos 9] 1 — g(R +h) (1 
a given dynamic pressure, there results a reduction in the de- _ ae) Se 
celeration load. In order to achieve this result, W/CpA _ _ ev? 4+ 2 
must be increased during that portion of the trajectory where 
peak dynamic pressure has not yet been realized, and de- ae — yw: 
creased when the peak has been passed. The physical [3] 
counterpart of W/CpA variation may involve a shape or Pa peat dt : 
among ee — eee a change in the drag co- where L/D is the lift-drag ratio, and the remaining notation 
ae p) or # change of the frontal area (A) of the used is illustrated in Fig. 1. Since it is not clear that Equa- 
vehicle. Jettisoning of depleted components which dimin- tions [1, 2 and 3] can be solved analytically unless some 
a ve — wl wide of course, also need to be con- simplifying assumptions are made, one may begin by defining 
= rg * © modulation may occur in a stepwise manner an area of interest in any given trajectory as follows: Each 
(sue as dropping a flared skirt at ed predetermined point in trajectory will be examined only to the point where the de- 
the trajectory) or <0 Aga continuously ’ such as by te celeration loads have become small enough to be of no con- 
pwhich prevail cern to the vehicle designer (or the pilot of the vehicle if it be 
which method of W/C>A variation is required. In addition manned). Subsequent to this point, the continued decelera- 
to treating the deceleration problem, the associated aero- ’ Numbers in parentheses indicate References at end of paper. 
EG ATS EA ‘In the analysis that follows, it has been implicitly assumed 
Received July 3, 1958. that the dynamic oscillations of the vehicle in no way affect the 
‘Member of Technical Staff, Astrosciences Laboratory. - motion of its center of gravity. For a discussion of the effects 
Member ARS. 7 of drag modulation on the entry dynamics of a rigid body see 
Associate Manager, Aerodynamics Dept. Member ARS. 
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tion will relegate the problem of successful recovery of the 
vehicle to the area of responsibility of the parachute and 
landing-system designer. With this philosophy in mind, one 
may examine more closely the equations of motion, partic- 
ularly the second term in the right member of Equation 
[1]. Notice first that this term will be identically zero at 
satellite velocities for which V? = g(R +h). The maximum 
satellite velocity (corresponding to minimum altitude satel- 
lites) is about 26,000 ft per sec. For velocities well above 
and below this value, the term [1 — V2/g(R + h)] is different 
‘rom zero and must be retained in Equation [1] if accurate 
solutions are to be expected. However, for Vz ~ 26,000 ft 
per sec, the inaccuracies introduced by neglecting this term 
lo not prove important, because the salient features of the 
trajectory under analysis will have already been determined 
luring the high speed portion of the trajectory, that is, be- 
‘ore the errors become significant. On the other hand, if one 
is dealing with entry velocities of 30,000 to 35,000 ft per sec, 
to neglect the aforementioned term in an approximate 
inalysis would yield unreliable results. This conclusion is 
slearly demonstrated in Fig. 2 which shows the results of 
“‘exact’’ machine solutions® of Equations [1, 2 and 3]. The 
peak g’s encountered during entry for a constant W/CpA 
trajectory are plotted as a function of entry velocity. Also 
indicated by dashed curves are the approximate solutions of 
the present analysis for the regions in which agreement 


with the exact solutions to within 10 per cent or better is Vg - Satellite Velocity 

achieved. It is apparent that as the entry angle is decreased Pr Vp- Retrorocket Velocity Increment 
the area of agreement is diminished. A fact not discernible ree ' V, - Resultant Vehicle Velocity 
from this figure but which is significant to the present analy- " oo - Fig. 1 Illustration of notation 


sis (and which has been demonstrated by the exact calcula- 
tions) is that the relative performance of different trajectories 42 he , 
with the same entry conditions is much more accurately * 


predicted by the approximate treatment than are the absolute 40 7 "oo 


values of the performance, particularly for low entry angles. L 
This will be further illustrated subsequently. /) 
With the above limitations in mind, the assumption is VIA 


1~ ——__|20 
g(R + h) 


32 
: 7 7 
Then, since this analysis is concerned with nonlifting trajec- 47 ea 


tories for which L/D = 0, Equation [1] reduces to LY 7 


Under these circumstances, the flight path angle @ remains 


\ 


essentially constant (and equal to its value at entry) during 
the portion of the trajectory which commands one’s interest. 
The utility of this approach was first demonstrated by Allen 


and Eggers (8). For entry velocities reasonably near 26,000 
ft per sec, the machine solutions of the complete equations 
of motion have shown @ to be constant within 5 per cent for 
the significant portions of the trajectories. 


Maximum Deceleration (g's) 


An additional simplification is introduced by the assump- 
tion that sin @ in Equation [2] is sufficiently small compared 
with the remaining terms so that it may be neglected. This 


12 


means that the component of the weight of the craft along EES — 

the flight path is small compared with the deceleration loads i 

the craft will experience; this is a good approximation, es- 
— 


pecially for small entry angles. Finally, the variation of 8 

gravitational acceleration with altitude is neglected. The sili vd 
foregoing assumptions and simplifications leave the following 
set of equations to be solved 


“Exect’ Calculations 


Approximate Theory 


| | | 


~ 64 2 2 30 32 34 


[4] 


Entry Velocity (feet/second x 10° 3) 


36 


g at ~ 2A 
’The computer programming for these solutions was per- Fig. 2 Maximum deceleration during a pure ballistic trajectory 
formed by J. Bundsen, whose assistance is gratefully acknowl- as a function of entry velocity. Entry angle arbitrarily selected 
edged. at 300,000-ft altitude 
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In effecting the solution of the above equations, it is mathe- 
matically convenient to change the independent variable 
from time (t) to the density ratio p/p» or a, and to introduce 
the parameter n, the deceleration expressed in terms of g. 


For the range of interest to the present analysis, the atmos- 
pheric omy v ariation is well approximated by 


Po 


If the constants are evaluated from the data of (4), the re- 
sulting values are 


2.6 X 10-3 (slug/ft?) = 1/23,350 (ft-) 


The deceleration can also be expressed as 


1 dV 1 dV do 
g do dt 


Equation [4a] becomes 


28 Asin 


For constant A this expression may be integrated to yield _ 


[6] 


and assuming gg = 0, one obtains the following expression 
for the velocity (V) as a function of o when A is constant 


ao 
Discrete Variation of W/CpA 


The simplest application of drag modulation a ypical 


entry trajectory will be considered first, that being a single 
stepwise change of W/CpA. Given an initial value of 
W/C DA, Vz and 6z, it is desired to find the point in the tra- 
jectory at which to increase instantaneously W/CpA to a new 
value such that the maximum deceleration during the entire 
flight will be the smallest of all possible values. One begins 
by considering a pure ballistic trajectory wherein no drag 
variation is introduced; i.e., W/CpA is considered to be con- 
stant at a value of Ay during the entire flight, as was done by 
Eggers and Allen (3). Then Equation [4a] is 


The deceleration n’ will be a maximum when 


y2 + 2a0V dV 
do 


d 
V2) = 
(oV?) 


Using Equation [7] for the velocity, one obtains for the values 
of o and V at maximum deceleration 


Ay sin 6¢ 


o = Om = ——— 


[8] 


[9] 


The maximum deceleration is 


gj 
poV x? sin [10] 


nN max = 


4ae 


Some of the inadequacies of this result, which are related to 
the simplification of the equations of motion, were discussed 
earlier with the aid of Fig. 2. In addition one may note 
that Equation [10] indicates that the maximum deceleration 
which occurs during a constant drag trajectory is independ- 
ent of W/CpA. Although this is true for a large range of 
values of W/CpA, 0g and V, pertinent to long range ballis- 
tic missiles, the complete machine solutions previously re- 
ferred to have shown n’max to depend upon W/C>p4A to the ex- 
tent that absolute values may be in error by 2 to 5g for extreme 
variations in W/CpA. Since the machine computations uti- 
lized the “real” variation of density as given in the ARDC 
model atmosphere (4), rather than the simple exponential 
variation (Eq. [5]) used in the approximate analysis, this may 
partially account for the trends just described. 
_ Proceeding, one may ask in what way the deceleration his- 
tory of a vehicle will be altered if at some point in the tra- 
jectory W/C>A is changed from Ap to some higher value, A; 
in a stepwise fashion (single step). The density ratio o = o* 
will characterize the point in the trajectory where W/CpA is 
to be modulated. Prior to this point, Equation [7] applies. 
To find an expression for the velocity when o> o *, one must 
return to Equation [4] 


Integrating, there results _ 
a(a — a*) 
i | A; sin 6x ] [ ] 


Now, one W ishes to determine the mz uxximum value of 


which the vehicle will experience subsequent to the point of 
drag modulation. 


Defining 
n= (- ) oV? 


one proceeds as before. (The subscript 1 denotes values after 
modulation, and the subscript 0 denotes values before.) 
Then, the density ratio and velocity at peak deceleration are 


2A, 


¢ - _ Ae sin 
2a 


where the peak deceleration is given by : 


Sketch A indicates the deceleration histories of a pure ballistic 
trajectory and of two drag-modulated trajectories with arbi- 
trary choices of density ratio, oa*, o»* at the modulation 
point. 

One may note from Equation [12] and with the aid of the 
sketch that, as o* is increased (thereby decreasing V*), 
Nim is decreased. This occurs because an increase in o* 
results in a sufficiently larger decrease in V* (which value is in 
turn squared), so that the net result is a decrease in Mijas 
Upon examination of Equations [10 and 12] and Sketch A, it 
becomes clear that the optimum location for o* is the one 


Po y= sin On 


4ae 
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which satisfies the equality 


[13] 


That is, with given values of Ap and A, the capability to re- 
duce the maximum deceleration is used to best advantage 
when W/CpA is increased from Ay to A; at the point in the 
trajectory such that the peak deceleration subsequently 
realized has exactly the same value as that which was exper- 
ienced at the last instant before W/CpA was modulated. 
Evaluating Equation [4a] at that instant yields 


= Nimax = 


= 

2A 
Combining this expression with Equations [12 and 13] one 


obtains 


sin | ] 
2A, 4ae A, sin 


Further, if one defines 


14 
A= A, sin On 
there results 
(3 [15] 
Ao 


Equation [15] may be solved for a range of (A;/Ao), thereby 
determining the optimum o* for a given set of entry condi- 
tions. Having determined the optimum value of o*, one can 
now proceed to find the reduction in maximum deceleration 
which can be realized if W/CyA is modulated at this o*. 
One may define the reduction in peak g’s as being character- 
ized by the fraction no*/n’,,,., Where n’,,, is the peak decel- 
eration load that would have been experienced had one chosen 
not to modulate W/CpA but to keep it constant at A for the 
entire duration of the trajectory. 
Then from Equations [10 and 4a] 


no* poo* 4ae 


n 2Ao po sin 
But recall 
= exp [ - Ap sin 
Therefore 
no* 2aa* 


| 


2aa* 
| 
Ao sin 6¢ 


ieee ver, by virtue of the definition of A; and Equation [15], 


= eX] 
Ao sin 


n max 


4=4, (constant) 


Sketch A  Deceleration history as a function of density eel 
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this may be arranged to yield 


no* Ai 


where the relation between A, and A,/Ao (Eq. [15]) is implied. 
Thus, as one might have expected, the maximum realizable 
reduction in peak g’s is a function solely of A;/Ao. These 
optimum solutions of A, and no*/n’,,,, as functions of (A,/ Ag) 
have been plotted in Fig. 3. Here it should be noted that 
Equation [15] has two solutions for each A,/Ao. However, 
the applicable solution can be selected simply by noting that 
_ the optimum o* must be less than om. That is, the modula- 
= tion must take place before n’,,,, is reached, or there will be 
no reduction in peak g’s whatsoever. 
Once the correct ¢* has been selected, it can be seen from 

Fig. 3 that there is a rather easily discernible point of dimin- 
ishing returns after which an increase in the W/C>)A ratio is 
not very profitable. It can further be seen that since no*/n’,,,, 
reaches a minimum value of only about 0.70, the application 

of single step drag modulation as a g-reduction scheme is 
limited to cases where peak deceleration loads need only to 
be reduced by 20 to 30 per cent. However, as will be seen in 
the next section, the use of continuous drag modulation can 
result in peak g reductions of as much as 50 per cent for rea- 
sonable values of Ao). 


Continuous Variation of W/CpA 


The opposite extreme of single-step drag modulation is the 
case which corresponds to an infinite number of infinitesimal 
steps in W/C>.A, that is, a continuous variation. These two 
extremes define an envelope which encloses all other possible 
ways of modulating W/CpA so as to reduce maximum decel- 
eration loads. The trajectory which results from the use of 
continuous drag modulation can be divided into three phases. 
Phase 1 will be the same as in the first part of single-step mod- 
ulation, namely, a portion of a constant W/CpA (pure ballis- 
tic) trajectory. In Phase 1 the velocity is given by Equation 
[7]. Phase 2 will begin at a point in the trajectory to be sub- 
sequently determined and at which o = o*. At this point, 
the deceleration will have reached a value which is to be main- 
tained constant thereafter by properly varyingW/C pA, specif- 
ically in direct proportion to the variation of the dynamic 
pressure, 3 p V%. The aforementioned value of decelera- 
tion load will be designated by no* since the maximum g force 
which will be experienced during the entire flight is the same 
as that value at which the modulation started. Finally, a 
point in the trajectory will be reached such that if one stops 
modulating W/CpA the deceleration load will drop off from 
that same value (n*). This point will initiate Phase 3 wherein 
W/C/_A is again held at a constant value, and, as in Phase 1 
Equation [11] with V* and o* replaced by V; and o describes 
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. 
0.6 
: 0. 
° 
a 


Ballistic Parameter Ratio, 3 | 


Fig. 3 Parameters for discrete modulation (single step) 
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the motion of the vehicle 


= exp 


alo — 
A, sin 


therefore 
B sin 0g 


Integration results in an expression for the velocity during 
Phase 2 


_2gno* 


or eliminating V** by means of Equation [4a] 
2no*A 2aa* 
V2 = (1 — _ hog, =) [18] 
poo* Ao sin 


The trajectory during Phase 2 must also satisfy the regular 
equation of motion [4a], which may be written 


Vy? = 2no* (A 
Po o 


where A is no longer considered to be a constant but is a func- 
tion of ¢. Eliminating V? between this expression and Equa- 
tion [18] yields 

A Ao 2a 


— — 
sin o* 


This equation describes the required variation of the ballistic 
parameter A with density ratio o during Phase 2 

Next it is necessary to determine the limits of Phase 2. It 
is desired to find the point in the trajectory, characterized by 
0, V; and A, where, if cessation of W/C'pA modulation occurs, 
the maximum deceleration subsequently experienced by the 
vehicle during Phase 3 of the trajectory will equal the value 
which was maintained during Phase 2, namely mo* 

The maximum value of deceleration that would occur after 
termination of modulation at any arbitrary density ratio o; is 
found by the same procedure as was used in the analysis of the 
single-step variation. Corresponding to Equation [8], one 
finds that the density ratio at that maximum deceleration is 
given by 


gy 
n/n’ 
max 
0.2 
10 


Fig. 4 Parameters for continuous modulation 


Substituting this expression into Equation [17] determines the 


velocity at this value of o 


ado, 


One is now in a position to determine from Equation [4a] the 
value of that maximum deceleration as 


26 9 

po V;? sin | 

Nimax da E Oe 


af = Vi exp 


The above requirement that n,,,,. = mo* can be used in con- 
junction with Equations [20 and 21] to evaluate o; and Vj. 
It can be shown that if the deceleration initially decreases at 
the end of modulation it will continue to do so. Since this is 
the only accepted deceleration pattern consistent with the 
above requirement, it follows that om,; = o1, so that Equation 
[20] yields 


o sinOe | 
= —— 

A, 2a 


[22] 
Substituting Equation [22] into Equation [21] then de- 
termines the velocity at the termination of Phase 2 


1 
posin 


Havi ing dete emined the limits of Phase 2 for specified trajec . 
tory and vehicle parameters, it is now possible to determine 
what values of no* may be realized or, more important, what 


values of m*/n’,,, may be achieved: Define 
9 * 
do = [24] 
Ay sin 


Then, at the end of Phase 1, the equation immediately pre- 
ceding Equation [16] still applies, and substitution of Equa- 
tion [24] therein yields 


no* = 

= exp (1 — do) [25] 


Now Equation [19], describing the variation of (A/a), must 
be consistent with Equation [22] which defines this ratio at 
the termination of Phase 2. Equating these expressions 
yields 

o 
= No loge No [26] 


1 
o* 


Furthermore, Equations [22 and 24] may be combined to 


give 
Ao 
so that Equations [26 and 27] yield 7 
A, 
28 
Xo exp 1) [28] 


From Equations [25 and 28] it is seen that Xo is a parameter 
which allows one to determine 79*/7' max in terms of A,/Ap ex- 
clusively. Furthermore, the density ratio at the start of drag 
modulation is given in terms of \» and the known values of a, 
Ay and 6, by Equation [24]. The values of mo*/n’max and Xo 
are shown as functions of (A,/Ao) in Fig. 4. 


Di 
iscussion 


Drag Modulation 


The most significant results of the preceding analysis are 
presented in Figs. 3 and 4. Indicated therein is the magnitude 
of the possible reduction in maximum deceleration load which 
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Table 1 Comparison of discrete and continuous modulation, Vg = 26,000 ft per sec, 0g = 0.1 radian, A, = 50 Ib per ft?, 
A; = 200 lb per ft? 


A> h*, ft 
Discrete 0.101 0.736 1.036 160,000 
modulation 
Continuous 0.270 0.560 0.692 x 10°% 170,000 
modulation 


tt ft 
per sec hy, ft per sec 
21,200 Not applicable to discrete modula- 12.2 16.6 

tion case 
1.026 Kk 10°? 107,000 


22,700 


11,800 9.3 16.6 


vehicle must experience, compared with that it would 
indergo in a pure ballistic (W/CpA constant) trajectory, if the 

‘ehicle possessed the capability of varying its ballistic parame- 
‘er W/C,A between the limits Ap and A;. The reduction de- 

vends on the ratio A,/ Ay and the program of modulation of 

W/CyA. As stated earlier, the two cases which have been 

‘onsidered (single-step and continuous drag modulation) con- 

stitute the boundaries of an envelope which encloses all other 

yossible manners of modulating to reduce peak g’s. 

In order to more fully illustrate the potentialities of drag 

‘ontrol as a g-reduction device, two representative examples 
have been caleulated. For an entry velocity of 26,000 ft per 
sec and an entry angle of 0.100 radian at an altitude of 300,000 
ft, a ratio (A,;/Ao) of 4 and a value of Ap = 50 Ib per ft? has 
been selected to show what values of peak deceleration may be 
expected when: 

1 W/C,A is increased in a stepwise fashion (single step). 

2 W/C_A is continuously increased to its maximum value. 
The velocities and altitudes which characterize the critical 
points in the trajectories, i.e., the points at which modulation 
begins and ends, are shown in Table 1. 

Finally, for the case of continuous modulation, Fig. 5 shows 
the various values of Vy, and @, which, when combined, will 
result in trajectories all of which have the same percentage 
reduction in peak deceleration. That is, for a given value of 
A\/ Ao there is a corresponding value of */n’max and hence 
Xo. Thus, for each constant value of n’max (Vg? sin O0¢ = 
constant) the value no* will also be the same. 

The dashed portions of the curves in Fig. 5 indicate the re- 
gions in which the present analysis is unreliable, the limit of 
reliability being the same as that defined previously, i.e., 
maximum g’s predictable to within 10 per cent of the exact 
values shown in Fig. 2. 

As was mentioned earlier, the use of discrete drag modula- 
tion must be restricted to somewhat limited applications, 
since the reduction in peak g’s which can be effected is ap- 
preciably less than that associated with continuous modula- 
tion. However, an attractive feature of the discrete modula- 
tion scheme is the probable simplicity of construction of such 
a device and subsequent adaptation to a given vehicle. One 
may visualize such a device as being a lightweight flare which 
would be fitted to the after portion of the vehicle and dis- 
carded by explosive bolts at the appropriate point in the 
trajectory. 

On the other hand, the use of continuous drag modulation 
may require a somewhat more complex mechanism to vary 
W/C)A in the necessary manner. This additional com- 
plexity, however, allows atmospheric entry under conditions 
which high g would otherwise make untenable. Continuous 
drag modulation could find application to manned vehicles 
returning from any satellite orbit, or even (as exact machine 
solutions of the complete equations of motion have shown) 
could allow a manned vehicle to safely enter the Earth’s at- 
mosphere at velocities as great as 35,000 ft per sec and entry | 
angles as large as 8 deg. Since the escape velocity of the Earth | 


is roughly 35,000 ft per sec, there would probably be few in- 7 


stances in which it would be necessary to enter the sensible 
atmosphere at velocities in excess of this figure. 


Of the many possible schemes which have been suggested — Fig. 6 Possible methods for varying drag 
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— Retractable 
t 


for varying drag continuously, three are illustrated in Fig. 6. 
The first scheme, (a), is simply a cruciform arrangement of 
flared panels which move in the indicated path. They might 
be actuated by a semipassive servo system which has a de- 
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celeration-initiated, predetermined time variation of W/CpA 
programmed into it; or it could be operated by a closed loop 
system having an “on-off” micro-switch which controls the 
position of the flaps so as to maintain a constant value of de- 
celeration. The second arrangement, (b), replaces the flares 
with paddles which retract into the interior of the vehicle, 
and which may be actuated by the same methods suggested 
above. Finally, (c) illustrates the use of a retractable spike 
which it might be possible to position so as to alternately in- 
crease or decrease the hypersonic drag coefficient in the proper 
manner, thereby achieving a variation in W/CpA. The suc- 
cess of such a scheme would depend on a sufficient knowledge 
of the flow separation induced by the spike at hypersonic 
speeds. The best method for modulating drag, of course, 


could only be determined from a system study with the point 


of view of minimizing the weight of the vehicle while retaining 
a maximum capacity of varying drag over a wide range. 

Here it should be observed that if a vehicle has the capabil- 
ity of increasing W/C>pA continuously from a low Ap to a 
higher value Aj, then it obviously possesses the capability of 
subsequently decreasing W/CpA from A, to Ap. In view of 
this fact, it should be evident that Phase 2 of a trajectory, 
wherein the vehicle drag is being continuously modulated, can 
be extended. This can be accomplished simply by continu- 
ously decreasing W/C pA in the proper manner, beginning this 
“reverse”? modulation at the point in the trajectory character- 
ized by o:, Vi, ete. This point is where peak dynamic pres- 
sure occurs. One possible advantage in extending Phase 2 of 
the trajectory would be to decrease further the vehicle velocity 
(at the prevailing high rate) an amount sufficient to lessen 
the demands on the recovery system. Ae 


Dispersion and Heating 


The discussion of any hypervelocity entry problem would, 
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Fig. 8 Relative heating due to continuous drag modulation 


of course, not be complete without including some mention of 
the attendant but no less important dispersion and heating 
problems. The studies of (5) show that when a vehicle leaves 
a satellite orbit by means of a retro-thrust, thereby establishing 
a transfer trajectory leading to the Earth’s surface, the great- 
est contributions to the dispersion of the landing point result 
from errors in the orientation and magnitude of the retrorocket 
velocity increment during the “deboost’’ period. Fig. 7 
shows these errors in range S per unit error in retrorocket 
orientation angle 6 and deboost velocity magnitude AV as a 
function of the entry angle 6, for an orbit at an altitude of 300 
nautical miles. From these curves it is evident that from the 
point of view of minimizing dispersion, it is desirable to have 
as large an entry angle as possible. However, Equation [10] 
has indicated that the maximum deceleration loads expected 
of a ballistic type of entry increase as the sine of the entry 
angle. Consequently, the reduction of dispersion to an ac- 
ceptable value may require unacceptable g-loads, unless the 
vehicle is provided with a deceleration control mechanism. 
Equipping the vehicle with such a device would permit it to 
enter the atmosphere at a relatively high angle, thereby keep- 
ing the dispersion low and at the same time holding the peak 
g’s to a tolerable value. 

Still within the systems frame of view, one must examine 
the effects of entry angle and W/C A on the total aerodynamic 
heating to which the vehicle is exposed. Several investiga- 
tors have shown (see (3)) that the laminar heating rate at the 
stagnation point of a blunt body of revolution is proportional 
to the square root of the free stream density and the cube of the 


free stream velocity thus 


[29] 
An expression for the ratio of the total heating during a drag 
modulated flight to the heat pulse which a vehicle would ex- 


perience during a constant W/C>)A trajectory with the same 
peak deceleration has been derived in (6). The result is 


(=) + 0.42 x 
nN max nN max 
A ( 1 ‘ 
38 do + 2) | 


This expression for Q has been evaluated and is shown in 
Fig. 8. It is seen that for modest g-reduction, the heating 
actually is reduced slightly (4 per cent for 18 per cent g-re- 
duction), although for appreciable modulation some increase 
in heating occurs (163 per cent for a 50 per cent g-reduction). 


Vp 


[30] 


6 This effect is magnified as the orbit altitude is increased. 
For a 500-nautical mile orbit, 08/05 decreases from 2460 nau- 
tical miles per deg at @¢ = 1 deg to —25.8 nautical miles per 


degree at 6 deg. 
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Concluding Remarks 


A scheme has been examined for limiting to a selected value 
the maximum deceleration load experienced by a returning 
space vehicle entering the atmosphere at a given flight path 
angle. Although in some cases reduction of the entry angle 
might be acceptable, in general, rather large entry angles (5 
deg and higher) are necessary if the landing point dispersion 
of these vehicles is to be held within reasonable limits. There- 
fore, drag modulation schemes may provide a solution to the 
problem of limiting the g-loads to vaiues tolerable by a human 
heing. The aerodynamic heating is not significantly affected. 

An interesting application of drag modulation which has 
heen considered is that of a returning moon vehicle entering 
the Earth’s atmosphere at about 35,000 ft per sec. As previ- 
ously pointed out, the limitations of the present analysis, 
when dealing with low entry angles, preclude accurate pre- 
dictions for the case of entry from the moon when absolute 
esults are desired. However, reliable values of 19*/7’max 
ean be obtained by the methods of this paper. This is evi- 
enced by the comparison of analytic and machine results in 
Cable 2. 

The results of machine computations for the case of entry 
rom the moon are shown in Fig. 9 where the absolute values of 
weak g’s during entry are plotted as a function of 6¢. If the 

ehicle attempts to enter at angles less than about 4 deg it 
vill not remain in the atmosphere but will ‘‘pierce’’ it and 
either return to the moon, assume a geocentric orbit, or per- 
orm a series of skipping maneuvers before finally descending 
to the Earth’s surface. For values of 0¢ greater than about 6 
deg, the deceleration loads experienced on a pure ballistic 
trajectory will be intolerable to a human passenger.’ As 
shown in Fig. 9, the use of drag modulation on such a vehicle 
might extend the ‘‘safe’’ range of entry angles from 5 + 1 deg 
to 6 + 2 deg. It will be appreciated that the protection 
afforded by widening this tolerance could easily make the dif- 


Nomenclature 

A = frontal or drag reference area of the vehicle, ft? 

Cp = drag coefficient of the vehicle, dimensionless 

D = drag force, lb 

‘ = base of the natural logarithm, 2.718. . . 

q = acceleration due to gravity, 32.2 ft per sec? at the 
Earth’s surface 

h = altitude above the Earth’s surface, ft 

L = lift force, lb 

n = defined as the quotient of the deceleration of the 


vehicle and the acceleration due to gravity g. 
Also referred to as g’s, dimensionless 
Qr = total aerodynamic heat pulse per unit area, Btu/ft? 
Q = ratio of the Qr for a drag modulated trajectory 
to the Qr for a constant W/CpA trajectory, 
dimensionless 


q = rate of aerodynamic heating to a vehicle, Btu/ft?- 
sec 

R = mean radius of the Karth, 2.09 107 ft 

S = distance traveled by the vehicle along the surface 


of a nonrotating Earth during the time elapsed 
from entry to impact—referred to as range, nau- 
tical miles 


0S/05 = error coefficient or the error in range S per unit 
error in 6, nautical miles/deg 

t = time of flight of the vehicle, usually taken as zero 
at entry, sec 

Vp = retrorocket velocity increment, ft per sec 

V = inertial velocity of the vehicle, ft per sec 

Vo = resultant inertial velocity of the vehicle immedi- 


ately after firing retrorocket, ft per sec 


7A tolerable upper limit on the maximum deceleration load 


has been arbitrarily selected as 20 g. 


2 


Table 2 Comparison of exact and analytical results 


(no*; n max) Q Q 
exact approx. exact approx. 
0.579 0.536 1.12 £23 
0.513 0.475 ae 1.22 
0.489 0.455 1.19 1.26 
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Fig.9 Maximum deceleration experienced by a vehicle returning 


Vs 


(W/CpA) 
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from the moon 


inertial velocity of the vehicle in a satellite orbit 
ft per sec 

vehicle ballistic parameter; weight divided by the 
drag coefficient times the drag reference area, 
Ib per ft? 

weight. of the vehicle, lb 

parameter defined to facilitate the solution of Equa- 
tion [21] 

atmospheric constant defined in Equation [6] as 
pog/28, lb per ft? 

constant used in Equation [5], ft~! 

mathematical symbol for the vehicle parameter 
W/C DA, |b per ft? 

angle, measured from the orbital flight path direc- 
tion, at which the retrorocket velocity increment 
is applied, deg 

flight path angle of the vehicle (see Fig. 1), radian 
or deg 

defined in Equation [27] 

defined in Equation [14] 

free stream or atmospheric density, a function of 
altitude, slug per ft’ 

reference density used in Equation [5], slug per ft* 

density ratio, p/po, dimensionless 


refers to conditions at entry (arbitrarily selected at 
h = 300,000 ft) into the Earth’s atmosphere 

refers to conditions in a drag modulated trajectory 
when A = Ao, the lower limit value of W/CpA. 
Also accompanies parameters which contain Ap 
in their definitions, e.g., 45. An exception occurs 
in the case of po, defined above. 
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refers to conditions in a drag modulated trajectory 
when A = Aj, the upper limit value of W/CpA. 
Also accompanies parameters which contain A, 
in their definitions, e.g., \1. 

= refers to a parameter evaluated at the point where 

maximum g’s occur 


* = refers to conditions in a drag modulated trajectory at which 
drag variation begins 

= refers to conditions in a trajectory wherein the drag is held 
constant at Ao (pure ballistic trajectory ) 


@ 


LTHOUGH the existence of radiation pressure has been 

known for many years, its application to spaceship 
propulsion has been mentioned only infrequently. A solar 
sail derives its propulsive force from the pressure due to the 
sun’s light falling on a sail, which may be a thin sheet of 
aluminum foil, or a thin plastic sheet silvered or aluminized 
on the sunny side to reflect the light. Some obvious ad- 
vantages of the solar sail are: Mass ratio of unity (mass of 
ship remaining constant), availability of sail force throughout 
ship’s entire journey, no need of fuel or propellant, no propul- 
sive powerplant aboard ship and its associated waste-heat 
disposal problem. Although the available force is small 
compared to that of chemical rockets, it can be applied for 
as long as it is needed. A trip to Mars or Venus can perhaps 
be made in less time by solar sail than by chemical rocket. 
A solar sail is equivalent to a rocket of mass ratio unity with 
an infinite propellant reserve. As far as propulsion is con- 
cerned, the ship can always return to Earth or make in-flight 
navigational corrections. Also, since the ship moves at a 
nonuniform speed under a noncentral force system, the 
ship’s contents are not absolutely weightless, although the 
weight of any object would be much less than that at Earth’s 
surface. In a recent paper (1),? Garwin states that the solar 
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Although the existence of radiation pressure has long been known, its application to spaceship 
propulsion has not been explored thoroughly. A solar sailing sip has a mass ratio of unity, needs 
no fuel or propellant, has no powerplant aboard ship and no waste-heat disposal problem. The 
propulsive force, although small, can be applied continuously in a controlled direction. The ship 
has an infinite “‘propellant”’ reserve to correct errors of navigation, or to return it to Earth at any 
time. A sailing project can be realized relatively soon. In this paper, the characteristics of a solar 
sail are examined in some detail, the equations of motion solved, the sail tilt angle optimized, the 
time of travel computed, and the performance compared with that of chemical-rocket and ionic 
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sail is of negligible cost, and is perhaps more powerful and less 
difficult than many often-cited competing schemes. We 
might add that a sailing project can be realized fairly soon 
without extensive research and development. In this paper 
the solar sail will be examined in some detail; the equations 
of motion will be solved; the time of travel will be computed, 
and the performance wil] be compared with that of chemical- 
rocket and ionic propulsion systems. Problems regarding 
departure from one planet and arrival at another will not be 
The Sail 
Perhaps the most important requirement of the sail is that 
it must have a low mass per unit area. The presence of 
holes, tears and other defects is unimportant. The radiation 
pressure on the sail can be calculated from Maxwell’s electro- 
magnetic theory of light or the quantum theory. On a nor- 
mally oriented reflecting surface at Earth’s orbit, the pressure 
is 


= 28o/c {1] 


If the sail is at distance r from the sun and oriented at an 


angle of incidence 9, the pressure is (2) = A 


[2] 


= po cos? 
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The numerical value of So is 1.94 cal/(cm?*) (min) or 1.35 X 
10° erg/(em?) (sec); c is 3 X 10” cm/sec; so that po = 0.9 X 
10-4 dyne/cm?. 

It is of some importance to compare the magnitude of radi- 
ation pressure with that of space drag, the latter being pro- 
portional to material density in interplanetary space. Space 
density is not known with certainty. There may be local 
variations. However, Alfvén states that (3) gas density in 
interplanetary space is probably intermediate between the 
density in interstellar space (1 atom/cm*) and that in the 
outer corona (105 atoms/cm#), and that a fair guessavould be 
10? or 10° atoms/em*. If we assume a mean density of 10° 
hydrogen atoms/cm‘, and a ship speed of 3 X 108 cm/sec, 
then space drag is on the order of 107 dyne/cm?. Since 
Alfvén’s estimate appears to be a reasonable one, we can say 
with some confidence that space drag is negligibly small com- 
jared to solar radiation pressure. 

There is a possibility that while the ship is in motion its 
sail may be punctured by meteoroids. Small punctures are 
inconsequential and need to cause no concern. Large punc- 
tures may need repair during the voyage. 

From the definitions given under Nomenclature, and re- 
inembering that @ is the acceleration due to radiation pressure 
joalone, we see that 


Hquation [3] is plotted in Fig. 1. The top line represents 
the hypothetical, limiting case where the sail material has no 
mass. The bottom line corresponds approximately to exist- 
ing technology, e.g., a sail material 0.1-mil thick with a spe- 
cifie gravity of 1.18. As an example, if a value of 0.2 cm/sec? 
is desired for a, and p = 2 X 1074 gm/cm?, then m/A = 
2.5 X 107-4 gm/em*. If m) = 500 kg, the required sail 
area A is 2 X 10° em?, which is equivalent to a circle 500 m 
(1650 ft) in diameter. With p = 3 >< 10-4 gm/cm?, the maxi- 
mum attainable value for @ is 0.3 em/sec?. This value in- 
creases to 0.45 cm/sec? for p = 2 X 1074 gm/cm?, and to 
0.9 cm/sec? for p = 10-4 gm/cm?. It is clearly important 
that the sail should have low mass per unit area. 


Equations of Motion 


Let us now consider a ship propelled by solar sail in inter- 
planetary space. The forces acting on it are shown in Fig. 2. 
F, is the sun’s gravitational force. F, is the force due 
to radiation pressure on the sail and is equal to pA. 
Space drag, being negligibly small, is not considered. 
We also neglect gravitational forces due to Earth and other 
planets. The equations to follow, therefore, are inaccurate 
when the ship is very close to a planet (say, within a million 
miles). The problems connected with departure from one 
planet and arrival at another are not treated in this paper. 
However, these problems are no more difficult with solar 
sail than with other means of propulsion. Under the speci- 
fied conditions, the equations of motion are 


—F, + Fs cos 0 v? 


F, A A ro\? ro\? 
Po costo = a cost ( ) [6] 
m m m r r 


The quantity F,/m can be written 


F F 2 2 
rs Fe (2) - [7] 
m m r 


JUNE 1959 


B 


T T 
12 
= 10 
8 
6 
e | 
4 
2 
ie) 2 4 6 8 10 12 rs 16 
sec2/em 
- 4 
Fig. 1 Characteristics of solar sail 
r 
Shi 
Light 
{ Fig. 2 Forces acting on spaceship : 


where do is the sun’s gravitational acceleration at Earth’s 
orbit. It has the numerical value 0.592 em/sec*. Substi- 
tuting Equations [6 and 7] into [4 and 5], we obtain, respec- 
tively 


ro \? v2 
+ 6) { -} =a—- {8 ] 
r r 
1 
—a sin cos? (") =v+ = [9] 
r r 


These equations are not difficult to solve in detail. For our 
purpose here, however, we shall make one simplification by 
neglecting the term u in Equation [8]. We shall show later 
that this approximation is a reasonable one for the type of 
motion here considered. The simplified equation reads 


y- 


v = — cos? [10] 


[8a] 


To 


(ay — cos? 8) ( 


which can be solved directly to give 


Eliminating v between Equations [9 and 10] and solving the 
resultant equation for vu, we obtain 


2a sin 8 cos? 6 


4 


423 


| 
84 - 
rat- 
igh 
rec 
8. 
te 
» 
ip 
ds 
he 
ip 
ay 
ar 
he 
‘SS 
Ve 
mn 
er 
ns 
d DA 
’ 
1g 
be 
[4] 
sin A uv 
of : =v+— [5] 
n r es 
. . . 
By using Equation [2], the quantity F,/m can be expressed 
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Fig. 3 Optimum sail setting and spiral angle of ship’s path 


Equation [11] can be integrated directly to give the time of 
travel between radii 79 and r 


1 — (ao/a — cos* 8)'/* 


From Equations [10 and 11], it is seen that Ah a Wait 
u| 2 sin 6 cos? 0 See 
[13] 
v a/a — 


The path of the ship is therefore a logarithmic spiral; the 
spiral angle Y = tan! u/v depends only on the acceleration 
ratio ao/a and the sail setting @. 

Suppose a planet is in orbit at a distance r from the sun. 
If the orbit is assumed to be circular, the orbital speed v, 0! 
that planet is given by 


> 
r 


Combining Equations [10, 11 and 14], we find 


w + 


Up Up 


2 sin2 [15 
8) 4 4(a/ao)* sin? 6 cos 


1 — cos? 0 


The ratio w/v, is independent of r. It depends only on a/ap 
and @. For a given sail design and sail tilt, the ratio of ship 
speed to planet speed, when the ship arrives at that planet’s 
orbit, is the same for all planets in the solar system. 

The power applied to the ship by the sail is 


P = (tangential component of sail foree)-v + 
(radial component of sail foree)-u 


which can be readily shown to be 


ay + a cos* 


P = —ma sin 6 cos? 6 Ze 16] 
; It follows that the work done by the sail in going from 79 to r is 
mro? (ao + cos? 8) /1 1 
Pap, 
va Optimum Sail Setting 
8 In order to go from re to r in the shortest time, Equation [12] 
shows that we must minimize the quantity 
This has been done. The result is shown in Fig. 3, which 
also shows the spiral angle Y computed from Equation [13] as 
/ a function of a for optimum sail settings. 


4 y 
.2 
0 
a,cm/sec® 
Fig.4 w/v, and u/(v?/r) at optimum 6 = 
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By differentiating Equation [11] we obtain u, which, when 
combined with Equation [8a], enables us to compute 


sin? 8-cos* 6 


[18] 
(ao — 6)? 

Equation [18] is plotted in Fig. 4 for optimum sail settings. 
It is seen that for a up to 0.25 cm/sec?, a is numerically less 
than 10 per cent of v?/r. Hence our original approximation 
of neglecting u in Equation [8] is a reasonable one. If very 


accurate results are desired, the term u should be retained, 
and the equations of motion solved more exactly. 

Fig. 4 also shows the ratio w/v, plotted from Equation [15], 
the sail setting @ being the optimum value corresponding to 
the value of @ used. 
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It is seen that for a up to 0.25 cm/sec’, 
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the difference between ship speed and planet speed is less than 
6 per cent, the ship being slightly slower than the planet. 


lrip to Venus 
As an example, let us consider a trip from Earth to Venus. 
The sail setting is shown in Fig. 5. The angle @ is taken to be 
positive when the setting is as shown in Fig. 2 or Fig. 5. By 
differentiating Equation [10] and using [11], we obtain 


v = (ro/r)? aw sin 8 cos? 6 [19] 


which is positive when @ is positive. Equation [11] shows 
that, for positive 6, u = 7 is negative, meaning that the ship 
will move toward the sun. Equation [17] shows that in 
moving from ro to a smaller radius r, the work done by the 
siil, W, is negative, meaning that energy is removed from the 


Fig. 5 Sail setting for Earth to Venus trip 
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Fig. 6 Trip time with optimum sail setting 
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Ll Fig. 7 Time of travel as a function of 6 and a 


ship. Therefore, as the ship moves in toward Venus, it gains 
speed and loses energy. It may be recalled that to bring an 
object from a larger orbit to a smaller one, it is necessary to 
increase the object’s speed and decrease its energy. The 
solar sail does this naturally. 

Substituting the numerical values ro = 1.495  10' em, 
r=r, = 1.081 X 10" cm, a = 0.592 em/sec?, in Equa- 
tion [12], and by using the optimum value of 6 for each value 
of a (Fig. 3), the time of travel t can be calculated as a func- 
tion of a. The results are shown in Fig. 6. A value for a in 
the neighborhood of 0.2 cm/sec? should be attainable in the 
reasonable future. Based on that value for a, the trip from 

arth to Venus takes 52 days. 

Fig. 7 shows trip time ¢ as a function of a@ and 6. The 
curves are calculated from Equation [12]. When a is small, 
the angle of sail setting is rather critical. 

It should be mentioned that the Venus orbit is inclined to 
the ecliptic at an angle of 3.394 deg. Therefore, unless the 
ship arrives at Venus’ orbit just when the planet is at one 
of its nodes (where the orbital planes intersect), the two will 
not meet, although the ship is at the correct radial distance. 
To insure contact, we must have at our command a force 
normal in direction to the ecliptic. We will not go into de- 
tailed calculations of such a force or the motion it causes. It 
suffices to say that such a force is readily obtained by prop- 
erly orienting the main sail or perhaps by using a smaller, 
auxiliary sail. Should the ship by miscalculation miss Venus 
on the first try, the sail will always provide the necessary 
propulsive force to turn it back for another try or return it to 
Earth. Since the ship is endowed with infinite ‘‘propellant”’ 
reserve, this advantage is not easily matched by other pro- 
pulsion schemes. 
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Fig. 8 Sail setting for] Earth to Mars trip 
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Fig. 9 Trip time along logarithmic spiral (¢) compared with that 
along transfer ellipse 


Trip to Mars 


Table 1 
for trip to Mars 


Solar 


sail 


mass of ship leaving circum- 
tellurian orbit/mass of ship 
entering circum-Martian 
orbit 

mass of ship leaving circum- 
tellurian orbit/mass of ship 
re-entering circum-tellurian 
orbit l 

mass of ship leaving circum- 
tellurian orbit/mass of useful 


load 1.8 
one-way voyaging time, days 118 
time to escape Earth several 

weeks 


several 
weeks 


time to adapt to circum- 
Martian orbit 


Comparison of chemical rocket and solar sail 


Chemical 
rocket 


9 (one way) 


73 (round trip) 


260 
negligible 


negligible 


For a trip from Earth to Mars, the sail setting is shown in 
Fig. 8. Here @ is negative. Hence, according to Equation 
[11], r is positive, meaning that the ship will move awa) 
from the sun. Equation [19] shows that v is negative 
Equation [17] indicates that W is positive, meaning that work 
is done upon the ship to increase its energy. As the ship 
moves out toward Mars, it loses speed and gains energy, again 
a desirable combination. 

By substituting ro = 1.495 X 10'% cm, r = r. = 2.28 X 
10'3 em, and a = 0.592 cm/sec’, the optimum trip time can 
be calculated from Equation [12]. The results are shown in 
Fig. 6. If we again take a value for @ of 0.2 em/sec?, the 
trip takes 118 days. By using chemical-rocket propulsion, 
it has been estimated (4) that the trip to Mars would tak 
260 days. It should be added, however, that the 118-day 
sailing time does not include the time required by the sail t 
escape Earth’s gravitation. If we assume the ship is first 
put into an orbit circling the Earth at a height of, say, 1000 
miles, it may take several weeks to escape from Earth. 

The Martian orbit is inclined to the ecliptic at an angle of 
1.850 deg. The remarks made regarding orbital inclination 
for the Venus trip are equally applicable here. 

For the return trip from Mars, the sail should be set in the 
same manner as for the trip from Earth to Venus (Fig. 5). 


Comparison With Chemical-Rocket Propulsion 


It has been shown that a solar sailing ship describes a 
logarithmic spiral. The time of travel between radii ro and r 
is given by Equation [12]. Another method of space pro- 
pulsion is to apply large thrust to the ship for short durations. 
By executing suitable thrust maneuvers (4), it is possible to 
make the ship follow a transfer ellipse whose major axis ex- 
tends from the Earth’s orbit (mean radius = ro) to the desti- 
nation planet’s orbit (mean radius = r). The ellipse contains 
the sun at one of its foci. If t, is the time required for a one- 
way journey along this ellipse, then, according to Kepler’s 
third law 

te = (20 


where g is the semi-major axis of the ellipse. Since q = 


(1/2)(ro + r), Equation [20] can be written as 


om 
(2)*/* ao To 


Combining Equations [12 and 21], and simplifying, we obtain 


t _ (2)%* {(ao/ax) [(do/a) — cos — 122] 


This is the ratio of time of travel along logarithmic spiral 
(i.e., by solar sail) to that along transfer ellipse. Equation 
[22] is plotted in Fig. 9, the optimum value of 6 being used 
for each value of a. To travel to very distant planets 
(Saturn and beyond), the solar sail would take longer, es- 
pecially when a@ is small. This is because the sail would have 
to spiral around the sun several times in order to reach a dis- 
tant planet. To go to nearby planets, however, sailing is per- 
haps the faster mode of transportation. The sailing ship 
appears attractive for a trip to Mars or Venus. 

A comparison of mass ratios brings out another important 
advantage of the sailing ship. Table 1 compares a trip to 
Mars by solar sail and by chemical rocket. 

A few words of explanation of Table 1 may be in order. 
The figures under chemical rocket are taken from (4). The 
term “useful load” means payload plus the necessary struc- 
tures and equipment to protect and preserve the payload. 


3 This is travel time to mean Martian orbital distance (2.28 X 
em). 
is 83 days. 


Travel time to perihelion distance (2.06 X 10! cm) 
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It is essentially the ship less the propulsion system. The ratio 
‘mass of ship leaving circum-tellurian orbit to mass of useful 
load” has the value 1.8 for the solar sail. This is based on 
the assumptions that a is 0.2 cm/sec? and the sail material has 
a value for p of 2 X 10-4 gm/cem?. For these conditions, 
Fig. 1 gives mo/A = 2.5 X 10-4gm/em*. Hence © 4 


A 


= 1. 
2.5 X 


mo Mo 


The author has not computed the time required by the sail 
to escape Earth, or that required to change the ship’s course 
from a spiral around the sun to an orbit circling Mars. A 
rough estimate is several weeks for each maneuver. When 
the time required for these maneuvers is included, it still 
appears possible that sailing might make better time to 
Mars. There is also a possibility that by using a combined 
solar sail and solar rocket, these maneuvers could be speeded 
considerably with relatively little increase of mass. 

The implications of the first three items in the table are 
clear. To send a 1-ton useful load to Mars and return it to 
ltarth, we would, if we chose chemical propulsion, have to 
erry a 73-ton mass to a space station circling Earth. If we 
chose to sail, the mass to be ferried would be only 1.8 tons. 
If the 1-ton load were abandoned to circle Mars forever, 
the ferried mass for chemical propulsion would be reduced to 
tons, which is still five times the 1.8-ton mass required for 
sailing. 


Comparison With Ionic Propulsion 


One scheme of spaceship propulsion that has received con- 
siderable attention is the ionic rocket. In this scheme an 
easily ionized material, such as cesium or rubidium, is first 
ionized and then electrically ejected behind the ship to 
create a thrust. The thrust is small but can be maintained 
for a relatively long period of time. This propulsion method 
has been studied by Stuhlinger (5, 6) and others. An ex- 
ample given in (5) shows that a ship with an initial acceler- 
ation of 0.0657 cm/sec? requires 124.5 days to escape Earth, 
222.5 days for the main journey from Earth to Mars, 54.3 
days for approach maneuvers when arriving at Mars, or a 
total of 401.3 days for a one-way trip. Thrust is applied 
during 390 of these 401 days. In (6), a performance chart is 
given for an ionic rocket with a propulsion time of 347 days. 
Using that chart, and for an initial acceleration of 0.0657 
cm/sec”, we estimate that the take-off mass is 2.5 times that 
of the useful load. 

If we chose to sail, the corresponding ratio would be 1.8; 
and the main journey would require only 118 days, as com- 
pared to 222 days by ionic propulsion. 

The comparison is a bit sharper for an Earth-Mars-Earth 
round trip. Excluding the time required for waiting and for 
escape and approach maneuvers, the main journey takes 
427 days by ionic propulsion vs. 236 days by sailing. The 
(take-off mass) /(useful load) ratio is 4.86 for the ionic rocket 
(5) vs. 1.8 for the sailing ship. 

These comparisons should be regarded as indicative but not 
conclusive, since they are based on examples which, although 
not unrealistic, are taken somewhat off-handedly. Never- 
theless, it seems clear that the solar sail possesses perform- 
ance characteristics at least on a par with the ionic rocket. 
It should also be remembered that the ionic engine is not a 


to carry a propellant, which is exhaustible. Surely, the 
solar sail also needs development work: Material and fabrir 
cation of the sail, its packaging, handling and control. But 
these appear to be relatively less difficult compared to what 
is required for the ionic engine. The author believes that a 
sailing ship is easier, less costly, and can be made operational ; 
sooner than an ionic rocket. - 


Concluding Remarks 


The solar sail seems to have important performance ad- 
vantages over either chemical or ionic propulsion for trips to 
Mars or Venus. It should be easier, cheaper and quicker to 
develop than an ionic engine. Some serious consideration 
toward the development of a sailing vehicle for interplanetary 


travel would be warranted. 

Nomenclature 
A = area of sail 

a = sun’s gravitational acceleration at Earth’s orbit; ao = 


Fgo/m 
c = velocity of light 
Fy, = sun’s gravitational force on ship at radius r 
Fg) = sun’s gravitational force on ship at radius ro ; 
"; = force on sail due to radiation pressure : 
m = total mass of spaceship (including sail) 


simple one. 


Its major components include: Nuclear reactor, 


= time derivative of @ 
= spiral angle = tan™! (w/v) 


my = mass of spaceship proper (excluding sail) 

P = power applied to ship by sail 

p = solar radiation pressure 

po = solar radiation pressure on a normally oriented reflecting 
surface at Earth’s orbit 

q = semi-major axis of transfer ellipse 

r = plane polar coordinate; radial distance from sun a 

* = time derivative of r 

ro = mean distance of Earth from sun 7 

ry = mean distance of Venus from sun A 

r2 = mean distance of Mars from sun 

So = solar constant at Earth’s orbit (normally incident energy 
flux outside Earth’s atmosphere) es 

l = time; travel time by solar sail - 

te = travel time by transfer ellipse 7 

u radial velocity = 

u = time derivative of u =» 

v = tangential velocity = r¢é 

v = time derivative of v 

vp = orbital speed of planet at distance r from sun 

W = work done by sail 

w = speed of ship = (u? + v?)'/2 

a = acceleration of ship due to radiation pressure po acting on 
sail area A; a = poA/m 

6 = angle of incidence : 

p = mass per unit area of sail 

@ = plane polar coordinate 

¢ 
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This paper presents a minimal set of equations applicable to radio detection of an ICBM or recon- 
_ maissance satellite and preliminary prediction of its path and impact point (if any). The input 
data, obtainable by modern radio tracking instrumentation systems, are assumed to consist of the 
position and velocity vectors of a “flying object” relative to a geodetic station. 
sufficient condition is derived for deciding whether or not the object is moving in a gravitational 
orbit. Equations are then given for each case; for a satellite, the elements of the orbit and the 
time of visibility are found; for a missile, the geodetic coordinates of the impact point and the time 
of impact are computed; for an aircraft, geodetic coordinates of position are found. The method 
may serve as the basis of a real-time computing system for identifying flying objects and predicting 


A necessary and 


their paths. 


INCE observations of an artificial Earth satellite or bal- 

listic missile obtained by modern radio tracking systems 
differ in character from those of classical astronomy, it ap- 
pears likely that the direct vectorial method (4)? leads to 
somewhat simpler determination of the orbit. The classical 
methods (5), e.g., those of Laplace and Gauss, are based on 
observation of angular positions relative to “fixed stars,” 
but most radio tracking systems give directly the position 
vector of the vehicle and its rate of change relative to the 
observing site. 

We shall present a minimal set of equations applicable to 
detection of an ICBM or reconnaissance satellite and deter- 
mination of its path and impact point (if any). The equa- 
tions, based on central-force theory, are “minimal”’ in that 
they suffice to determine a preliminary orbit or an impact 
point, but they are not intended for very accurate long- 
time predictions. 


General Theory 


Let the (U, V, W) orthogonal Cartesian axes have constant 
directions relative to the fixed stars, with origin at the cen- 
troid O of the Earth, and the W-axis positive toward the geo- 
graphic north pole. The Earth rotates with constant angular 
velocity Q about the W-axis. (Effects of precession are 
negligible in our problem.) The observing station S on the 
Earth has coordinates (U, V, W), which, as functions of time 
t, are 


U pcos 
V = psin Qt | {1] 
Sy 


where p and S,, are constants determined by the figure of the 
Earth. 

At station § there is an orthogonal coordinate system (S,, 
S2, S3) with the S;-axis positive upward along the normal to 
the Earth, and the S,- and S.-axes directed east and north, re- 
spectively. It is in this (S;, Se, S3) system that the radio 
instrumentation system measures the position vector P of the 
vehicle M and its velocity P. 

If S; is the unit vector along the S;-axis, we assume that M 
is observed by the instrumentation system during a set of 
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time values (¢;) which satisfy the “visibility” condition (1) 


P.S5;> 0 [2] 

Within an approximation sufficient for our purpose, we 

assume that the gravitational field of the Earth is spherically 

symmetrical, in which case the path of M is a fixed plane 

curve (conic section), relative to the (U, V, W) system. Let 

R and R be the position and the velocity vectors, respectively, 

of M relative to the centroid O of the Earth in the (U, V, W) 
system. 

According to Milne (4), the orbital motion of a particle in 

a central gravitational field is described by the vector equa- 


tions 
[3] 
[4] 


where H and B denote vector constants of integration, and 
M = 1.407735 X ft?/sec? 


denotes a constant of the central field. Generally, » = 
G(M, + m), G = universal gravitational constant, My = 
mass of central body, m = the (negligible) mass of satellite. 
Physically, 7 represents the angular momentum (3) per unit 
mass. 

The orbital curve is a conic section with one focus at the 
centroid O of the Earth. The vector B points from O along the 
major axis of the conic in a direction away from the nearer 
vertex (perigee) of the conic. The semi-latus rectum | and 
the eccentricity e are given by ae 

e= [6] 


Furthermore, the conic is an ellipse or a hyperbola when 
e< 1 ore> 1, respectively, and the length a of the semi- 


major axis is 
[7] 


Substituting Equation [3] in Equation [4] and eas tities 
the vector triple product, we find 


B = ((u/|R|) 


= {I/(1 — e?) | 


Evaluation of B 


- 
> 
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Since R=P+S5 [9] 
it follows that in order to compute B by Equations [8 and 9], 
we should express P and S in the (U, V, W) rege Then by 
diffe [9], we find R in the (U, V, W) sys- 
tem as 


[10] 


so that by substituting Equations [9 “wih 10] i in Equation [8] 
we can compute B in terms of known quantities. We now 


consider the necessary transformation. 


Coordinate Transformation 


An orthogonal transformation exists between coordinates 
in the (S;, Se, S3) system and those in the (U, V, W) system. 
Let 


S = position vector of station S relative to O 
Q = angular velocity of Earth rotation 
p = distance from W-axis to station S 
S, = distance from the (U, V) plane to S 
@ = geodetic latitude (2) of station S 

t) = initial time of observation 

Ao = longitude of station S at time fo 


measured in the (U, V) plane clockwise from the U-axis. 
Let (S,, So, 83) and (€ V, W) be unit vectors along the re- 
spective axes. Then the angle 7 between U and the projec- 
tion of 5 in the (U, V) plane, measured counterclockwise 
from U, is given at time t by 


= Q(t t do 11 
n 0) [ ] 
Hence 


S = (pcos n)U + (psin n) 7 + SLW [12] 
and 
Si —sin 7 cos 0 
| = | —sin @ cos — sin @ sin n cos V [13] 
5; cos@cosyn cos @sinyn sin @ W 
Since the observed data consists of 
= + + P38; {14] 
and 
PS + PAS: + PS; [15] 


we can now express P and P in the (U, V, W) system by 
Equations [13, 14 and 15]. 


Hence by Equations [9, 12, 13, 14 and 15] we find 


R=R,O0+ RV + [16] 
where 
R, = —P, sin — P2sin cos n + P; cos cos + pcosn 
R, = P,cosn — Pz sin sin n + P3 cos sin 7 + psin7n 
R, = P:cos@ + P;sin + Sy 
Differentiating Equations [16 and 17] we find 
R=R.O+ + RW [18] 
where 
R, = (—P, + QP2sin cos @ — p)] sin 7 + 
[—P,Q — Pz sin + P; cos cos 
R. = [-QP, — Pzsin + P; cos $] sin + 


— sin @ — P3 cos @ — cos 
Rw = + P;sin 
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. . . . . 
tively, and let B; and be their respective time derivatives, 


For each observation (P,, P,) and corresponding time & 
where, depending upon the computing facilities available, k 
may or may not include all points covered by the index a 
occurring in the “General Theory,” we evaluate R and & by 
Equations [16 to 19], and then compute ff and B by Equa- 


tions [3 and 4]. 


By the method just given, we thus generate a set of values 
To decide w or not the observations 
(P,, P,) are descriptive of a gravitational orbit, we may com- 4 
pute the derivative of a vector Z;(t,) by a finite-difference — 
formula, e.g., if the ¢,’s are equi-spaced at interval At, then 


Test for Orbit 


a=-+n 
3 aZa 
a=-—n 


Z = — (2s) ln= 
(At)n(n + 1) (2n + 1) 


(20) 
Let B, and R; denote unit vectors along B and R, resepec- 


as computed by Equation [20]. From Equation [8] we de- | 
rive a necessary and sufficient condition for existence of a — 
gravitational orbit, namely, that the position and velocity 
vectors must satisfy the inequality 
[Br | < [21] 

To see that condition [21] is necessary, we note that if the — 
gravitaticnal field of the Earth were spherically symmetric, — 
then the position R and velocity F of an orbiting vehicle must _ 
be such that B, by Equation [8], remains constant, and since _ 
R is never zero in the orbit, then evidently inequality [21] 
holds. The actual gravitational field of the Earth is not 
spherically symmetric, but even in this case we note that |B; | 
and |R;| represent the angular rate of rotation of B and R, 
respectively. The nonspherical part of the gravitational 
field causes |B, | to differ from zero, i.e., the line of apsides, or 
B, rotates slowly with time, but from the theory of perturba- 
tions ((5) pp. 333-334) we know that the change in B is quite 
small per revolution of the orbiting body, and so again in- 
equality [21] holds. 

To see that inequality [21] is sufficient, we consider a 
vehicle in nonorbital motion, say an aircraft or spaceship 
moving at constant speed and altitude along any great-circle 
route, and we use the position and velocity vectors R and R 
of the vehicle to compute a vector B by Equation [8]. Since 
R-R = 0, and since the vehicular speed differs from that of an 
orbiting body at the same altitude, Equation [8] shows that 
B is not zero and that B rotates at the same absolute angular 
rate as R, ie, |Br| = |R,|. If now the vehicle deviates 
from its course (as by starting to dive, climb or turn), then 
B will start to rotate faster than R, because Equation [8] 
shows that, in this case, B has nonzero components along 
both R and R. The nonspherical part of the gravitational 
field will not sensibly alter the preceding statement. Thus, 
when a vehicle is in nonorbital motion in the vicinity of the ' 

Earth, the vehicular position R and velocity R satisfy the 
condition 


|B, | > [22] 


Hence, a vehicle is in orbital motion about the Earth if, and 
only if, the vehicular position and velocity vectors satisfy 
inequality [21]. 


Determination of Orbit 


After deciding that a gravitational orbit exists, we are 
naturally interested in computing the elements of the orbit. 
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From the previously computed set (Ai, Br, ts), we estimate From the definition of J and H, it follows oe 


K pointing along the minor axis of the orbit is 
B. given by the vector product 
k K=H@B 


and a vector J pointing along the line of intersection of the 


Then using Equations [23] in Equations [5, 6 and 7], we find orbital and equatorial planes is given by > 
the semi-latus rectum /, the eccentricity e and the length a J=few - on 
of the semi-major axis. Then an equation of the orbit is ae 
s } Hence, since Q is the angle between vectors K and J, it is given 
= |R| (1 — e cos 8) we [24] by 


and since L is the angle between U, it is given by 
As a check on the calculation, we have 


L = are cos [29] 


using the minus or the plus sign when e < 1 or e> 1, respec- If it has been determined that the vehicle is moving in a 
tively. 5 _ gravitational orbit, we must decide whether or not this orbit 
it ice 7 intersects the Earth. Since the eccentricity e and the semi- 
Angular Orientation of Orbit (7) cod ae major axis a of the orbit have been found by Equations [6 and 
= = ; 7], the distance d from the center of the Earth to the perigee 
Although the vectors H and B in Equations [23] suffice to : siete ase 
point of the orbit is given by 
orient the orbit, usually the orientation is specified by three 
angular coordinates (J, L, Q) defined as follows: d = |a(1 — e)| [30] 
IT = inclination of the orbit. Wetake/ positive when meas- 
ured from the equatorial plane toward the orbital plane on 
the descending side of the orbit. 


Let ro be the equatorial radius of the Earth. Then if e< 1, 
the orbit is elliptical, and the vehicle will collide with the 


= longitude of a node (ascending or descending). In Earth if 7 “. 
Fig. 1 we show L as positive measured westward in the equa- d<ro o-., re 
torial plane from the U-axis to the line OF which points to whereas if ta - 
the descending node F. [32] 


Q = angle between the major axis of the orbit and a nodal ; f : , 
line. We take Q positive Ww hen measured from minor axis of the vehicle moves 1n a sate lite orbit about the Earth. 


the orbit to the line OF of the descending node F. Otherwise, if é > . the orbit is hyperbolic, and if the in- 
; ~ equality [31] holds, then the vehicle will collide with the 


~ Earth or recede into outer space according as |R | decreases or 

i increases in time, respectively. If e> 1 and the inequality 
_ [32] holds, the vehicle recedes indefinitely. 

: If impact conditions prevail, we determine the time and 
4 location of the impact. Otherwise, we find the time of 


i <n At any time ¢ of observation, we sii the polar angle @ by 

a - Equation [25], the eccentricity e by Equation [6], and the 

~i ~—-—C Semi-maajor axis by Equation [7]. If e< 1 (elliptic case), we 
compute 


E = 2are t (= 33 
= 2 arc tan [33 ] 


itr and the period of the orbit is 


Then the time 79 of perigee passage is found from Kepler’s 
equation as 


T =t — (T/2m) (E — esin E) 7 [35] 


eae) 1 (hyperbolie case), we compute 


2 are tanh “tan 5 136] 


Fig. 1 Curve GG, surface of the Earth. Curve EE, equator. 
Curve M,M, track of vehicle. Point S, observation station t+ (a*/ V/ uw) (F — esinh F) [37] 
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Visibility of Orbital Points 


After evaluation of 70, all six elements (a, e, 70, /, L, Q) of 
the orbit are known; i.e., the orbit is determined. Then the 
values of the times at which the vehicle is visible from any 
observing site may be computed as in (1). 


Prediction of Impact 


Assuming that impact conditions prevail, we may find the 
mpact point by constructing a vector equation for orbital 
points and then imposing the condition that such a point also 
coincide with a point on the spheroid representing the 
arth. 

By Equations [3 and 4], we have already evaluated the 
two vectors H and B which determine the orientation of the 
orbit. Another vector K, lying in the orbital plane, is con- 
veniently given by the vector product 

K=H@B [38 | 
Let Bo and Ko denote unit vectors along B and K, respectively. 
Then the position vector R of an orbital point may be ex- 
pressed as 
R = (By cos + Ky sin 4) [39 | 
can be found from Equations [16 and 17] as 


R\| = |R2+ R,2|'7 


where |R 


and the angle @ (the “true anomaly’) is defined by Equation 
[25]. Combining Equations [24 and 39] gives 


= {//(1 — e cos 9)] (By cos 8 + Ky sin 8) [41] 


Since H, B, K are vectors with known numerical com- 
ponents in the (U, V, W) system, we have 


bO+bV+6,0 By? =1 
Ro = + kV + k3 IV Ko? = | 


where b and & are known numerically. Hence, Equation 


[41] becomes = 
R=R0+ RV RW 
R; = (b; cos @ + k; sin 8)l/(1 — e cos 8) G = 1, 2, 3) 


If the surface of the Earth is represented by the spheroid 
+ — =} [44] 
the value of 6 corresponding to the impact point is, by Equa- 
tions [43 and 44], a solution of 
C? (R,? + + = AC? [45] 
We find that Equation [45] becomes 
q cos? 8 + q2 sin 20 + cos + =0 [46] 


where 
b;? e? 
qi = + be? — ky? — ke? + /2 
A? . 
q = bik, bok + bsks 
2A%e 
= 


Although Equation [46] can be rearranged in the form of a 
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+. Regardless of the type of vehicle, we are naturally inter- 


fourth degree polynomial equation, in cos @, it can probably | 

be more conveniently solved by iteration (e.g., Newton’s — 

method) in its present form, especially since we are inter- 
ested in just one root. 

After 0, 0 = 6, has been found by Equation [46], the © 

radius vector R, of the impact point is given by Equation [43]. 

_ The geocentric latitude y, of the impact point is determined 


by 
R,-W 
which, by Equation [43] with 8 = 6,, becomes = 2 
Wp = are sin [48] 


Hence the geodetic latitude of the impact point is aisle \, 


as 
= are tan [(A2/C?) tan [49 


where A and C are the semi-major and semi-minor axes of | 
the spheroid [44] representing the surface of the Earth. 
In order to find the longitude of the impact point, we must | 
first determine the time t, of impact, which can be found =~ 
from Equations [33 to 37]. 
Thus, knowing the true anomaly 6, of the impact point, — 
and having already found the time 7 of the perigee, we sub- 
stitute 6, for @ in either Equation [33] or [36] (according as 
é< 1 ore> 1, respectively) and solve either Equation [35] or 
[37] (as required), to find t = t,, the time of impact. 
At the time of impact, the azimuth 7, of the station S in 
the (U, V) plane is by Equation [11] 
Np = — to) — Ao [50] 
Hence, by Equation [43] with 6 = 6,, the longitude A, of 
the impact point measured westward from S (with S as 
origin) is 


A» = are tan (R2/R1) — np [51] 


Geodetic Coordinates of the Vehicle 


ested in the geodetic coordinates (6) of the projection Gy, of 
the vehicular position on the surface of the Earth. 

Since the position vector R of the vehicle at any time is | 


M 


Fig. 2 Vehicle’s geocentric and geodetic latitudes projected on 
the surface 
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given in the (U, V, W) system by Equations [16 and 17], 
then the corresponding geodetic position vector is 

G=R-A’N [52] 
where the vector N is the unit vector normal to the surface of 
the Earth at point Gy and pointing to the vehicle, and h is the 
altitude of the vehicle above the Earth. 

Let @ and wy be the angles between the equatorial (U, V) 
plane and the vectors hN and R, respectively. Let A and C 
be the semi-major and semi-minor axes, etenaiinds of the 
spheroid representing the Earth. wri 

The geocentric latitude is given by ae 


y = are sin |R}) [53] 
The geodetic latitude is given by ; 
= are sin{ + (At — C)R,2)'7} [54] 
The altitude h of the vehicle is given by 
|R| cos (@ — ¥) — — | sin? y) 


where 
|G |? = (At cos? @ + C* sin? @)/(A? cos? @ + C? sin? 


The longitude A of the point Gy as measured westward from 
station S with S as the origin is 


arc tan (R,/R.) — 7 
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Effect of Thrust Termination 


~ Process Upon Range Dispersion 


of a Ballistic Missile 


ARNOLD J. KELLY' 
Rocketdyne Division 

North American Aviation, Ine. 
Canoga Park, Calif. 


The contribution to the range dispersion of a ballistic missile caused by the thrust termination 


characteristics of its powerplant is analyzed and discussed. 


The range dispersion attributable to 


the nominal cutoff impulse acting on a variation in missile burnout mass is found to have a marked 


mee a effect as compared to a variation in cutoff impulse from nominal. 
_ reasonable estimates of engine and missile parameters. A generalized graph is presented which 


These results are based upon 


_ permits the estimation of the combined contribution of cutoff impulse deviation and burnout 
mass variation to range dispersion to be made once the nominal burnout conditions of the missile 
at engine are known. A numerical example of a hypothetical 500-nautical mile ballistic missile is 


presented. 


N GENERAL, the ultimate miss distance of a ballistic 
missile from its intended impact point can be attributed 
to errors associated with guidance, target mapping or power- 
plant. These errors are reflected as a variation in the burnout 
velocity vector or position from the intended nominal. This 
supposition assumes, as will be the case in the analysis, that 
the trajectory of the center of mass of the missile is not in- 
fluenced by any external influence other than gravity after 
burnout (no terminal guidance and vacuum re-entry). If 
mapping errors are neglected, the range dispersion of the 
missile is a function of irregularities in the burnout velocity 
vector and position caused by the guidance and powerplant. 
This study is concerned with the influence of the power- 
plant’s thrust termination characteristics upon the inaccu- 
Received Nov. 12, 1958. 
— Engineer, “Advanced Design Section. 
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racies in position and velocity at the end of powered flight on 
the dispersion of the missile. The relationship between the 
cutoff characteristics of the engine and the performance of 
the missile is then of particular importance, for exam pe 

1 It permits the calculation (for a given missile) of the 
thrust termination characteristics of the powerplant for a 
given contribution to range dispersion. 

2 Conversely, the range dispersion for a missile powered 
by an engine having specified cutoff characteristics can be 
estimated. 

The former is of importance, for the thrust termina- 
tion characteristics can be specified and the engine designed 
accordingly during the early stages of the powerplant de- 
sign. 


Discussion 
During the thrust termination process of a ballistic missile’s 
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HOW THE MILITARY SERVICES 
use Air Products CAPACITY 7 


Large-scale low-temperature systems, ultra pure gases and 
liquids, and a broad range of specialized cryogenic ‘‘hardware”’ 
are supplied by Air Products to the military. When large quantities 
of liquefied gases were needed for rocket engine development 


and missile testing, Air Products quickly designed, manufactured 
and put on stream complete production facilities. Typical facili- 
ties paid for themselves in less than a year’s time. Air Products 
also provides a broad line of portable air separators for field and 
shipboard use... and has advanced the development of exotic 


fuels. And, Air Products produces advanced design liquefied-gas _ 


pumps, cryogenic storage and transfer systems, electronic cool- | 
ing devices and refrigeration and distillation equipment for 


Air Products 

ER feedstock 
“preparation and 
oduct pliriiicatron 


Liquefied gases from Air 
Products equipment supply 
all major U.S. missiles. 


OU will find tangible of 
a growing technology. Applying “Cryo- 
genics” (the science of low temperatures ) 
and engineering broad new routes to low-cost, 


_ high-purity industrial gases is the main busi- 
~ ness of Air Products. 


Air Products combines original research 


_ knowledge with engineering and manufactur- 


ing capabilities and substantial operating 
These activities have 


major facilities provided by Air Products: 
Bethlehem Steel Braz fational Stee! - Celanese: 
Bmical « du Pont * Grace Chemical - Great Lakes Stee 
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HOW THE STEEL INDUSTRY 
uses Air Products CAPACITY 


In the blast furnace, the open hearth and the new 
converter processes — Air Products oxygen effi- 
cently increases steel mill capacity. Annealing 
rn trogen and other gases are also provided on a 
i w-cost tonnage basis. 

Air Products’ complete gas supply systems are 
stalled at steel mills without capital investment 

operating worries on the part of the users. 
yntinuity and reliability of supply are assured. 
n-site facilities pioneered by Air Products 
r duced the cost of oxygen 80% in 12 years — 
t ansforming oxygen from a costly chemical to 
practical working utility. 

Further progress marks on-the-job development 
vork now continuing around the clock at major 
s‘eelmaking facilities. Entirely new metallurgical 
techniques... and new profits... are available 
t rough Air Products. 


A Ao 


of .elped provide many Air Products customers 

o- vith distinct competitive advantages. 

3) Air Products is the world’s leader in 

st, APPLIED CRYOGENICS — the practical 

3i- ind profitable use of low-temperature science 

‘or industry. 

*h Perhaps this CAPACITY can help solve 

r- your problems — in cryogenics, in industrial 

ig gas supply systems, or in some new area where 

je “ground rules” are yet to be established. 
Air Products pipeline oxygen 
serves the Basic Oxygen Fur- 
nace Process. 


Laughlin Steel Osaka Oxygen - Spencer Chemical 
WS. Government-All major research and defense 


United States Steel + Venezuelan Ministry of 
irton Steel. 


HOW THE CHEMICAL INDUSTRY 


§ uses Air Products CAPACITY HIGH PURITY GASES AND LIQUIDS 


Air Products low-temperature systems permit — available like any other utility 
many modern chemical plants to improve operat- with the help of Air Products — 
ing efficiency and end-product quality — and to 

develop new processes and products. This results CAPACIT 
from the ready availability of low-cost tonnage 
quantities of oxygen, nitrogen, hydrogen, ammonia 
and methanol syn-gas, carbon monoxide and 
hydrocarbons such as purified methane, acety- 
lene and ethylene. Low-temperature separations 
of gaseous mixtures now make it practical to 
recover valuable components from natural gas, 
refinery off-gases, coke-oven gas and other 


“‘waste”’ gases. The versatility of cryogenics—as 
applied by Air Products—works profitably for the a. s 

chemical industry today... offers unparalleled 

future opportunity in this fast-growing industry. oan a P 
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... The supply is dependable... 
the price guaranteed .. . with ee 
Air Products on the job. 
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To find the new Applied may 
in the manufacture, separation and purification of Industrial 
Gases. 


To provide integrated research, engineering, manufacture and 
operation . ..complete services under a single responsibility. 


To forecast accurately and guarantee cost, 


performance and reliability. 
To put entire gas supply work for you 
capital investment on your part. : 

Air Products CAPACITY has helped our customers to step 
out ahead of competition in familiar fields...to open up 
entirely new areas of opportunity through new products or 
processes. A letter or telephone call will put Air Products 
CAPACITY to work for you. 


INCORPORATED 


ALLENTOWN, PENNSYLVANIA 


Information is 
available on: 

uperc ductivity 
frozen free radicals 
biplogical preservation 


environmental « control 
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owerplant, two distinct phenomena occur which contribute Pree 
to the range dispersion. These two effects are: — <i 

A deviation in the cutoff impulse from nominal. ; SEAS. THRUST CREAT 

A variation in the burnout mass from nominal acted on by 
the total cutoff impulse. 

The latter effect is, in reality, attributed to an inadequacy 
present in virtually all guidance systems which, lacking pro- 
‘ision to sense net thrust, cannot compute the instantaneous 
mass of the missile. 

The following analyses will present a method for assessing ete © INITIATION 
the change in the velocity vector because of these two effects . OF CUTOFF 
and will relate this to range dispersion. ; 


Effect of Cutoff Impulse Deviation on Range Dispersion 


During the thrust termination process of any rocket engine, 
liquid or solid, variations in sequencing or actuation cause 
the engine to exhibit an uncertainty in its thrust decay char- 
acteristics. The cross-hatched region in Fig. 1 illustrates this 
uncertainty. In order to relate this uncertainty to a vari- 
ation in burnout velocity, consider the following development. 

Writing the fundamental rocket equation (1)? 


Fig. 1 Typical thrust decay characteristics 


Mo 
AUs = Ue — Ue = Teco g In (*) —gsinf(te —t) [1] 
the change in velocity along the thrust vector during the cut- 
off process can be determined. Noting that 


and AM < Mo then 


AU; = —gsin — ty) 
Mo 
It should be noted that the time (t,) when the thrust decays 
to zero is usually poorly defined, because of asymptotic tail-off 
characteristics, or oscillations near zero thrust. Therefore, 
t. should be defined as the time when some percentage of 
nominal thrust is reached (say 5 per cent) where the remain- 
ing impulse due to tail-off or oscillation would be small. LOCAL HORIZONTAL 
The change in the velocity increment because of a devi- 
ation in cutoff impulse from nominal (A/,,) will be denoted Fig. 2. Velocity vector diagram 
by 6(AU,). It is assumed that t. — ¢, is invariant; is 
constant at the nominal value. This assumption is allow- 7 
able, as the normal variation to be expected in ¢. produces a : 7 : 
variation in 6(AU,) which is, for all practical cases, one order 
of magnitude smaller than that due to the variation in the * pts: 
first term of Equation [3]. Furthermore, a variation in ¢, in t > 0 for vehicles flying with a constant pitch attitude at 
essence prolongs or shortens the burning time of the missile 
and shifts the position of the burnout point in space from its ; The ran ‘hange because of the velocity magnitude vari- — 
nominal position. This effect is in opposition to the change a6 
in velocity because of At., so that the change in range is ation is then, combining Equations [5 and 6] 
0 = COS So ‘ 
KAU, Tey + Meo ro) OF My, \OUo hued 
Mo M The exchange factor O0R/OU> is a function of the missile — 
Defining the parameter a = Al«/Ieo, then Equation [4] re- parameters at burnout. _ The neunine burnout conditions will : 
: ie yes be used because the variations in these parameters are small. 
Guces to : ns The derivation of this exchange factor is shown later. 
6(AU;s) = aleo/Mo [5] Not only the magnitude, but the direction (6) of the ve- 
locity vector changes because of 6(AU,). Since tan 6(@) = 
As noted, this change in velocity is in the direction of the 5(O), and (Us) < Us, then 
thrust vector. It is reflected as a change in range because of 
two effects: Variation in magnitude, and variation in direc- 5(6:) = 
tion of the velocity vector (U). aia 
In Fig. 2, the change in the magnitude of the velocity vector 
(U) because of 6(AU,) can be derived, to a very close approxi- 
mation as shown below OR ale. sin & [OR 
8(U.) = 6(AU,) cos & (6] (ca) 
The exchange factor 02/06) is evaluated in the same manner 


l 


[8] 


The range change corresponding to this effect is then 


2 Numbers in parentheses indicates References at end of paper. 
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As noted above, because 6(Uo) < Uo, the factor 


as OR/OU>. 
How- 


5(Uo)/Uo would be expected to be exceedingly small. 


ever, because the magnitude of 0R/0# is unknown, the rela- 


tive magnitude of the range dispersion caused by the direc- 
tion ch: ange Equation [9] compared with that due to the 
change in the magnitude Equation [7] will have to be eve alu- 
ated numerically. 

Thus the total range dispersion because of a deviation in 
the cutoff impulse of the engine denoted as AR,, can be writ- 


where 


= 1/(ro + ho) 
Zr = 1/r 


The change in angular displacement is related to the range 
in nautical miles (from the burnout point) by 


ten as follows cine 
AR = 3438 [15] 
AR; = | * da 20. [10] Taking the partial derivative of the range equation with 
respect to the burnout velocity Uo yields 
etting 
oR oF sin 2g ok 6876 + 2ac — + + Uo?) | 
aU, cos s + 20, 2U, {11] Uo(b? + 4ac) V +(a/c) (b/c)Zn Zk 
Zr(b? + 2ac — WUo%) + (a + 
ar, . V +(a/c) + — Zr? 
Mo : and the partial with respect to the angle 6) at burnout is | 
Effect of Total Cutoff Impulse on Range Dispersion 3[(b? + 2ac)Z, + ab] i 


The typical radio-inertial or all-inertial guidance system 
senses acceleration and not net thrust; therefore the instan- 
taneous mass of the missile cannot be calculated. 


for variations in the burnout mass. 
because of uncertainties in structure weight and propellant 
residue. Propellant residue is by far the most significant. 
For example, if it is possible to have +1 per cent of the gross 
propellant as residue (or sliver), then a missile with a mass 
ratio of 10 could have approximately +9 per cent variation 
in burnout mass. This variation in the propellant residue 
can be caused by an engine which is running within specifi- 
cation limits. 

It can be seen that an engine having perfectly reproducible 
cutoff impulse, i.e., A7.. = 0, can still introduce an error in the 
cutoff velocity because of the above mentioned variation in 
the burnout mass. 

Considering Equation [12], the nominal range contribu- 
~ tion during the cutoff process can be written, for a missile 
having nominal burnout mass, as 


The change in range from nominal because of a variation in 
burnout mass will be designated as AR. Considering that a 
missile which is heavier than nominal at burnout will fall short 
of the intended impact point, for a given total cutoff impulse, 
then 


ad 
Derivation of the Exchange Factors tay, : 


In this study for the trajectory after burnout, the following — 


assumptions will be used: 

The Earth is spherical. 

The Earth is nonrotating. 

The effects of the atmosphere are neglected. 

These assumptions simplify the derivation of the ex- 
change factors and but the errors intro- 
duced by their use should be reasonably small and in keeping 
with the intended use of this work, i.e., the rapid estimation 
of the effect of thrust termination characteristics upon range 
dispersion. 

The range equation from which the partial derivatives 
are taken as exchange factors is Kepler’s problem in classical 
mechanics, which can be written as (2) 


—2cZp + =) 


This fact 
prevents the guidance system from making a correction 
These variations occur 


06) tan Bi(b* + +(a/e) + — 


2ac)Ze + ab 
V +(a/c) + (b/c)Zr — Zr? 


Once the nominal conditions at burnout have been as- 
— certained, it is then possible to calculate the contributions to 
range dispersion because of the cutoff process (AR; and AR). 
Equations [16 and 17] have been left in the form as shown 
because they do not simplify when the functions for a, 6 and 
c are inserted. 


Total Contribution to Range Dispersion Due to the Cutoff 
Process 


The two factors contributing to range dispersion are inde- 
pendent, because they stem in part from two separate aspects 
of engine operation. The variation in cutoff impulse is a 
function of the sequencing and actuation times of the cutoff 
process. The burnout mass variation is largely dependent 
upon the engine operation (performance) during mainstage, 
which is totally independent of conditions during cutoff. 

These two effects will be assumed to have normal distribu- 
tions which permit a root mean square (rms) sum to be made 
to take both into account. 

«Letting denote the rms range dispersion, then 


B \? 
= V(AR,)? + (AR)? = (-£,) [18] 
M, + i+, 


Another parameter ) is designated here to permit the de- 
velopment of the generalized plot shown in Fig. 3, where 


B 
= 2 = 19 

_ An interesting feature of Equation [19], which can be seen 
in Fig. 3, is that range dispersion is not symmetric in B as it 
isin a. Range dispersion is the same for positive @ or nega- 
tive a, while it is greater for negative 8 than for positive B. 
As mentioned in connection with Equation [13], a missile 
having a positive B will fall short of the intended impact 
point. Because of this, the following convention will be 
used: Positive 8 will always be associated with negative o 
and vice versa. 

It should be mentioned that the total cutoff impulse /,. re- 
fers to altitude conditions as does AJ... Because virtually all 
engine testing is conducted at or near sea level, it is desirable 
to be able to predict the engine’s altitude cutoff character- 
istics on the basis of these static tests. One method of ac- 
complishing this would be to use the chamber pressure decay 


M 0 
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characteristics to estimate the altitude thrust termination 
characteristics as shown in Equation [20] : 


F = 


This method is approximate for two reasons: 


1 The chamber pressure decay characteristics under flight 


conditions may vary from those recorded during static testing 
because of acceleration, ete. 

2 The thrust coefficient during the decay process varies 
because of fluctuations in the chamber mixture ratio. All 
tables of C; are determined for steady-state conditions, whereas 
the cutoff process represents a transient condition. 

However, this method should be adequate for preliminary 
design purposes. 


Consider a hypothetical missile having the following nomi- 
nal burnout characteristics: 


Sample Problem 


M, = 100 slugs 


= 10,000 fps 
= * 

6 = 30 ae, = 10 deg 

= 40 deg 

ho = 200,000 ft 


Teo = 5000 Ib-sec 


gross propellant wt. = 400 slugs 


Determine the maximum allowable cutoff impulse deviation 
for a range dispersion of 3 nautical mile (n. mi) and a pro- 
pellant residue variation of +1 per cent of the gross propellant. 

Substitution of the above characteristics into the range 
equation indicates that the missile wil] impact approximately 
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B=AOMo/Mo 


Fig. 3 a vs. for various 


a 


0.05 0.10 


A= Mo /Qleo 


480 n. mi down range of the burnout point. Then 


= 0.117 n. mi/fps 
OR/08 = 57.3 n. mi/rad 


Q =O0116 
d =0.08 
Referring to Fig. 3 using the above parameters 
8 = +0.04 
. a = +0.0769 for B = +0.04 
= +0.0752 for B = —0.04 


The minimum value for a (0.0752) which will give the 
maximum allowable range dispersion must be used. If the 
other value a were to be used (+8), the maximum allowable 
range dispersion would be exceeded. Then a = 0.0752. which 
corresponds to AJ, = 375 lb-sec. 

For this particular case 


¢ 
AR, = +0.436 n. mi 
and 


= 247 — 0.22: 
AR, = +0.247 — 0.223 n. mi 

The asymmetrical effect of 8 upon range dispersion can be 
clearly seen from the above numerical example. Maximum 
range dispersion corresponds to a negative value of B (the 
missile is lighter than nominal at burnout). 

It is significant to note that even if the cutoff impulse 
deviation could be reduced to zero (a = AIl,, = 0), the 
missile’s range dispersion would only be approximately 
halved. The same effect could be secured by simply halving 
the total cutoff impulse keeping @ and 8 constant at the 
values given above. 
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Using Fig. 3, the required thrust 
termination characteristics of a powerplant can be deter- 
mined quite rapidly for a given missile and mission with 
sufficient accuracy to satisfy preliminary design requirements. 

The total cutoff impulse (J..) acting on a variable burnout 
mass (AM,) can have a marked effect upon the range dis- 


to any ballistic missile. 


persion of a ballistic missile. The most effective manner in 
which to reduce range dispersion because of the cutoff char- 
acteristics of the engine is to reduce the total cutoff impulse 
rather than attempting a precision cutoff where only AI,. 
is reduced. A rapid cutoff process is the most desirable in 
effectively reducing this contribution to range dispersion. 


Nomenclature 


thrust chamber throat area, in.? 


thrust, lb 


QQ 
il 


g = gravitational constant = 32.2 fps? 
= height above a spherical Earth, ft 
Teo = cutoff impulse, Ib-sec 


specific impulse, sec 
average specific impulse during cutoff, sec 
missile mass, slug 

thrust chamber, chamber pressure, psi 
defined in Equation [11] 

distance from the center of a spherical Earth, ft 
radius of a spherical Earth = 20.903 X 105, ft 
distance measured along Earth’s surface, n. mi 
time, sec 
velocity, fps 


=~ 
~ 
| 


AMo/Mo 

angle measured between local horizon and thrust vector 

angle measured between local horizon and velocity vector 

defined by Equation [19] 

angle measured between thrust and velocity vector 

rms range dispersion, n. mi 

angle around Earth from burnout, rad 


oil 
b 
~ 
3 
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Three possible modes of droplet behavior have been ex- 
amined to find the requirements for consistency with 
Godsave’s equation for liquids burning or decomposing in 
spherically nonconvective systems with 
evaporation rate controlling. The empirical diameter- 
squared rule for small droplets is invalid for surfaces of 
low curvature, while the constancy of mass consumption 
rate per unit area observed for flat surfaces is untenable for 
small droplets. However, the assumption of a constant 
distance from liquid surface to flame front leads to results 
which are consistent with both situations. Agreement 
with some experimental investigators is shown. It is 
suggested that if the aforementioned assumption is valid 
over a large enough range, consumption rate data obtained 
from liquid strand burners can be used to predict droplet 
behavior. 


symmetrical, 
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Godsave’s Equation and the Diameter-Squared Rule 
HE equation derived by Godsave (1)* for the evaporation 
of liquid from a spherical liquid droplet surrounded by a 
concentric hot surface is 


p 


tins = (1) 


le 


Equation [1] is rigorously true only for the ideal situation 
which it is intended to describe. In real systems, complica- 
tions and deviations must be expected because of convective 
effects and variations in gas properties caused by the tempera- 
ture gradient across the gas film. However, it appears to be 
a good approximation in a number of cases of practical interest. 
If mass consumption rate per unit area is considered rather 
than total mass consumption rate, and quantities which are 
constant or nearly so for a particular combination of liquid 
and atmosphere are lumped as B, Godsave’s equation can be 
written as 


tay 


A 
A number of investigators (1,2,3) have found that the 


behavior of burning or evaporating droplets can be repre- 
sented empirically by 


D? = D,? — k't [3] 
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4 Relation of Droplet Consumption Rates 
eae to Liquid Strand Consumption Rates 
: 
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which implies that the total mass rate of consumption is pro- 
portional to the radius of the droplet (the so-called ‘‘first- 
power law’’), or that m,/A is inversely proportional to the 
droplet radius. Heat transfer considerations indicate that 
behavior of this sort might be expected, since it seems quite 
reasonable that either the ratio of hot surface area to liquid 
surface area will increase, or that gas film thickness will de- 
crease as droplet radius decreases. Conversely, it does not 
seem reasonable that the distance between the hot surface and 
the droplet surface should increase as droplet radius de- 
creases, which would be required if m;/A were to remain 
constant. 

If Equation [3] is to be consistent with Godsave’s equation, 
r./r; must remain constant. This behavior does not appear 
implausible, since it means simply that the system would 
maintain geometric similarity during combustion. Penner 
(4) claims some theoretical justification for this. However, it 
is apparent that Equation [3], with its corollary of constant 
r./T1, cannot be applied to large droplets or, in the limit, to 
systems in which the hot surface area is equal to the liquid 
surface area. It is inconceivable that the distance from the 
hot surface to the liquid surface would increase indefinitely. 
Furthermore, since m;/A must be inversely proportional to 
the radius of the liquid surface, the impossibility of combus- 
tion in the limiting type of system is implied. This conse- 
quence is contnary to experience, since it is well-established 
that combustion can take place from flat surfaces, and de- 
composition can take place in liquid strands (i.e., columns of 
liquid confined, in tubes) where the hot surface area is very 
nearly equal to, or possibly even smaller than, the liquid sur- 
face area. It is known that in systems of this kind, a constant 
distance from hot surface to liquid surface and. a constant 
m,/A are established. It was shown above, though, that the 
constancy of m,;/A is very unlikely for small droplets. 

There might appear to be, therefore, two distinct types of 
droplet behavior, depending on whether the gas film between 
the liquid and the hot surface is thick or thin. Further, 
direct extrapolation from either type of behavior to the other 
is not satisfactory. The distinction between the two types 
appears to lie in the mathematical approximations which have 
been employed rather than in the actual behavior of the drop- 
lets. There is undoubtedly a gradual transition between the 
two types, and it is probable that neither approximation ap- 
plies very well in this range. (It is also possible, but less likely, 
that either approximation would be satisfactory in the transi- 
tion regions.) The transition range is probably important in 
practice. 


Consequences of Constant Distance to Flame Front 


A somewhat more consistent model can be developed by 
using the observed characteristic of burning and decomposi- 
tion from nearly flat liquid surfaces, that the distance from 
the hot surface to the droplet surface remains constant. God- 
save’s equation for total mass consumption rate can be 


written 
Br(C + 11) 
ms = 


== 


where C is the distance from the hot surface to the liquid sur- 


[4] 


face. Also 
av 

np = pp — = 4 5 
Combining Equations [4 and 5] and separating variables 7 4 

ridr, _ (6) 
C+r pC 
Integrating, and evaluating the integration constant by 


letting ¢ 0 


Bi 
(C +r) = pC +ro —C In (C + rw) [7] 
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which implies that, when r, is large compared to C, r; will be 
a nearly linear function of t, and m,;/A will be nearly constant. | 
In other words, the behavior of liquids evaporating from drop- 
lets would become increasingly similar to that of liquids 
evaporating from flat surfaces as droplet radius is increased, 
which is a desirable result. 

The derivative, Equation [6], can be written in the form 


dr; 


dt 


from which it can be seen that, if 1/C can be neglected in com- 
parison with 1/r;, an equation having the form of Equation 
[3] will be obtained upon integration. If this is done, it is 


7 


= 
p 


The equations based on the assumption that the distance from 
the flame front to the liquid surface remains constant have, 
therefore, the advantage that they are consistent with ob- 
served results for the limiting cases of both small droplets and 
flat surfaces, and may offer some possibility of correlating 
results for intermediate types of droplet burning or decom- 
position. 

The table summarizes the implications of constant r. — r; 
for the two extreme ranges of liquid radius. Oddly enough, 


k’ = [9] 


Table 1 Consequences of constant distance to flame 


front 
Behavior of 
m A 


Behavior of liquid 
diameter with time 
D? = D,? — k't 
D Dy, — k’t 


Range of r; 
rl < 
rr >C 


my/A 
m,/A = constant 


it appears that either the assumption of constant r./r,; or 
the assumption of constant r. — r; leads to the result that r/? 
should be a linear function of time for small droplets. 


Comparison With Various Experimental Results 


In the investigation of Hall and Diederichsen (3), the times 
required for complete burning of droplets of varying initial 
diameter, rather than variation in diameter of single burning 
droplets, were measured. This method has a decided ad- 
vantage over the suspended droplet methods because of the 
wider range of variation in droplet radius which it is possible 
to study. A somewhat more reliable empirical description 
of the droplet behavior is to be anticipated for this reason. 
For example, the data given by Godsave (1), obtained from 
droplets burning in still air at atmospheric pressure, appear to 
be almost as linear when the first power of droplet diameter 
is plotted against time as when the square is plotted. This 
is because the droplet radii did not vary by a very large ratio. 
Hall and Diederichsen, by plotting the log of droplet lifetime 
against the log of initial radius for droplets of widely varying 
initial radius, came to the conclusion that droplet lifetime was 
proportional to the 2.05’ + 0.15 power of initial radius, and, 
therefore, that the droplets were behaving according to Equa- 
tion [3]. The observation was also made that r. — 7; was 
constant. To get some idea of how the data might fit the 
equations derived from the constant r. — r; assumption, the 


droplet lifetime from Equation [7] is 


If values of ¢ calculated from Equation [10] in the range from 
C = 0.5 ry to C = 5 ry (very approximately the range indi- 
cated by Diederichsen and Hall’s measurements of flame 
front distance at the forward stagnation point) are plotted vs. 
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rio on log-log coordinates, the droplet lifetime is found to be 
proportional to about the 1.8 power of initial radius, which 
appears to be in fairly good agreement with the experimental 
value. 

On the other hand, the data given by Agoston, Wood and 
Wise (5) for several materials burning on fairly large (up to 
1 cm in diameter) porous spheres do not fit any of the equations 
based on the constancy of a single parameter (r./r,, m,/A, 
re — 7). To quote the authors, ‘At constant pressure the 
distance of flame at the forward stagnation point from the 
surface of a sphere or suspended drop decreases with de- 
creasing diameter, but the ratio of flame diameter to drop or 
sphere diameter increases.”’ 


Conclusions 

From this brief study, it appears that the lack of a reliable 
way to predict the variation of r./r; or r- — r; with r; is a serious 
obstacle to the development of a completely analytical de- 
scription of droplet burning. However, if the constancy of 
re — Tr, is acceptable as an approximation, droplet consump- 
tion rates can be estimated from liquid strand burner data. 


Nomenclature 

my = mass rate of consumption 

A = average conductivity of gas between the hot surface and 
the liquid surface 

r = radius 

Cy = average specific heat of gas between the hot surface and 


the liquid surface 
AL = heat of vaporization of the liquid 


T = temperature 
B= a lumped constant 

D = droplet diameter 

t = time by, 

k’ = evaporation constant 

V = liquid volume 

p = liquid density a 

Subscripts 


= hot surface condition 
= liquid surface condition 
0 = initial condition 
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Chemical Reactions in Supersonic Flow’ 


F. A. WILLIAMS? 

Harvard University, Cambridge, Mass. 


ECENT experimentation first reported by Nicholls (1)* 
has demonstrated that the establishment of a stationary 
normal shock in a combustible gas stream may initiate a 
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chemical reaction in the flow behind the shock. In this ar- 
rangement the pressure and temperature rise across the shock 
wave are sufficient to ignite the combustible mixture, forming 
a reaction front which cannot propagate upstream into the 
supersonic flow. The obvious similarity between this phe- 
nomenon and the detonation produced by igniting a reactant 
gas in a tube has prompted the reasonable definition 
“steady detonation.” In addition to the fact that these 
new detonations are stationary in the laboratory reference 
frame, they differ from previously observed detonations 
in that the final state is not determined by the Chapman- 
Jouguet condition (Mach number M = 1). The observation 
that these combustion waves are strong detonations (M < 1 
finally) (2) invalidates previous inconclusive arguments that 
non-Chapman-Jouguet detonations are unstable (3). The 
chemical reaction appears to have a relatively small effect 
upon the shock Mach number (1) which is largely determined 
by the aerodynamic configuration. Thus, unlike deflagra- 
tions, there does not exist a unique detonation velocity. It 
might therefore be expected that in a suitable arrangement, 
weak detonations (M > 1 finally) may be produced, in which 
case combustion would occur in supersonic flow. It is the 
purpose of this note to present tractable governing character- 
istic equations for chemical reactions in steady supersonic 
flow and to offer an incomplete hypothetical answer to Gross’ 
challenge to find an application for supersonic combustion 
(2). 

For supersonic flow, it is reasonable to neglect the effects of 
viscosity, diffusion and thermal conduction in comparison to 
those of convective transport. It can then be shown that the 
general equations of mass, momentum and energy con- 
servation, chemical reaction rate, and state can be written, re- 
spectively, as (3, 4) 


a? Dt Cp jut 


pDv/Dt = —Vp 
Ds 
om 
Dt pl j=1 
DY; _ w; 


p 
p= pp, 

p = density aie 

velocity 
p = pressure vat 
7 = temperature 
8s = specific entropy 
Y; = mass fraction of chemical species j of the n-component 
mixture 

‘w; = mass rate of production of component 7 per unit volume 
oy; = chemical potential of species 7 


Here a? = (Op/Op)..x;, B = and cp = 
(0h/OT),.y; represent, respectively, the square of the sound 
velocity, the volumetric thermal expansion coefficient, and the 
specific heat at constant pressure for the mixture, all evaluated 
for the chemically frozen system; /A is the total enthalpy of 
the mixture, and h; = (Oh/OY;)p.p,y«ixj) equals the absolute 
specific enthalpy of species j for the case in which the mixture 
is an ideal gas. 

For steady-state flow, by using thermodynamic relations 
and the momentum equation, the entropy equation (Eq. [3]) 
can be replaced by an equation for the enthalpy A which is 
often more useful in combustion applications. The steady- 
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Fig. 1 Supersonic flow with combustion about a flat plate 


a. 


state forms of Equations [1 through 5] then become - 
~ 
pVev = - Vp - hyw; 
a’ Cp j=1 
po: Vr = —Vp [7] 
ve V(h + V2/2) =0 [8] 
= w;/p = 1, ,n) [9] 


[10] 


where V = 

In the case of plane two-dimensional or axially symmetric 
flow, this set of equations is easily expressed in finite difference 
form. If & is the distance along a streamline, 7 is the coordi- 
nate normal to £, M = V/a represents the Mach number re- 
ferred to the frozen sound speed, @ is the angle between the 
local flow direction and the axis of symmetry, and r is the 
radial distance from the axis of symmetry, then it is found 
that 


AE 


= “Ap = (M2 — 1)'2 x [11] 

AV = —Ap/pV for An =0 [12] 
_ A(h + V2/2)=0 for An =0 [13] 
AY; = (w;/pV) A&G = 1,...,n) for An =O [14] 
= pe, ¥:,..., Ye) [15] 


where the factor N in Equation [11] is zero for the plane case 
and unity for axially symmetric flow. These n + 4 equations 
constitute a tractable governing set for chemical reactions in 
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4 which illustrates that exothermic reactions 


supersonic flow and may be solved by characteristic tech- 
niques quite similar to those used for rotational flow. Equa- 
tion [11] determines the manner in which pressure changes 
propagate along Mach lines, while Equations [13 and 14] il- 
lustrate that the total stagnation enthalpy remains constant 
along streamlines, and composition changes propagate only 
along streamlines. 

A further simplification arises for flow in the neighborhood 
of a flat surface, in which case Equations [11 through 15] re- 
duce approximately to the one-dimensional flow relations for 
a reacting gas. In particular, changes in @ are small, and 
Equation [11] becomes approximately 


= hyw;A 16] 
(iM 

. 
hyw; > 0 

j=l 

in supersonic flow tend to increase the pressure, a well-known 
result which may conceivably be used for Jift or thrust aug- 
mentation. In the classical example of supersonic flow about 
a flat plate, fuel might be injected into the supersonic stream 
from the leading edge of the plate. The pressure and tem- 
perature rise across the oblique shock may ignite the mixture, 
producing a flow pattern for which typical streamlines and 
Mach lines are represented by the lines AB and AC, respec- 
tively, in Fig. 1. The flow along the flat plate is completely 
determined by the approach flow crossing the shock between 
the nose and the end point D of the limiting Mach line DE. If 
L is the length of the flat plate, p2 is the pressure above the 
flat plate, and p, is the pressure directly behind the oblique 


shock, then according to Equation [16], an approximate ex- 
pression for the fractional increase in the lift coefficient is 


— 1l)ep] | jw; 
fm) 


1 
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On Turbulent Supersonic 
| Diabatic Wakes 
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N THE past, the calculation of a theoretical model for the 
diabatic base region in transonic and supersonic flow (1)? 
had been somewhat hampered by its numerical complexity. 
The recent release of excellent charts of auxiliary integrals (2) 
greatly simplifies the non-isoenergetic jet mixing portion of the 
‘alculations. The purpose of this note is to illustrate results 
which may now be obtained easily through the use of the new 
generalized jet mixing information (2). Using the new charts, 
sample calculations of the base pressure ratio for a diabatic 
Received March 2, 1959. 
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wake were completed, and typical results are presented in this: 
note. The computations were carried out for the restricted 
theory (i.e., the theory for thin approaching boundary layers) 


step. Integral methods (2) were used to evaluate the energy 
addition to the wake from the free stream. 


Although tractable, the theory for calculating the base pres- _ 


sure ratio of a back-step is not easily condensed to a single 
simple explicit equation. Hence, the theoretical result is in- 


dicated here functionally sed 


P,P 
P, P; 


where 
M, > 1.0 
k = the ratio of specific heats ek! - 
H = the dimensionless bleed number bain Eq. [2]) 
0.9 
too 
2 #4 
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Fig. 1 Influence of bleed gas rate and temperature on a wake 
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HEAT TRANSFER RATE ,cA 
Fig. 2 Influence of heat transfer on a wake with zero bleed 
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: 


as it applies to supersonic flow past a two-dimensional back- _ 


the ratio of bleed gas stagnation temperature to 
free stream temperature 

dimensionless heat transfer number (see Kq. [3]) 

the similarity parameter for the homogeneous co- 


As 4 ordinate of the fully developed turbulent mixing 


profiles 

In order to show the flexibility of the theory, two examples 
are presented, one with no heat transfer from the walls but 
having mass addition to the wake, and a second with heat 
transfer from the walls but no mass addition to the wake. 
Both calculations were performed for an approach Mach 
number (immediately before the base) of 2.00 and a specific 
heat ratio of 1.40. 

Fig. 1 shows the effect of bleed mass rate and temperature 
on the base pressure ratio. Bydefinition 


¢ Wr 
_ Gs WT. ( 
w here 


4 Gp the mass bleed rate per unit width at which gas is 


Top/T 0a 


HP. 
added externally to the w ake 


= the step height 
R = the gas constant 
= the gravitational constant 


Ta 


the free stream stagnation temperature 
the stagnation pressure of the free stream before the 
base 


=e 


Physical appreciation for the dimensionless bleed number can 
be obtained by realizing that the quantity 


k 1 (k + 1)/2(k- 1) 
( +1) 


is the actual mass rate per unit width at which gas is added 
externally to the wake, divided by the maximum mass rate 
(from free stream initial conditions) per unit width at which 
gas could be added externally to the base. 

Fig. 2 illustrates the effect of heat transfer on a closed 


wake. By definition 
vf R\'* (k + 
where nat. 
cy» = the specific heat at constant pressure 
Q: = the heat transfer rate per unit width from the adjacent 


walls to the wake 


The other symbols are identified above. The dimensionless 
heat transfer rate may be interpreted as the actual heat trans- 
fer rate (per unit width) through the walls adjacent to the 
wake, divided by the maximum free stream convective energy 
rate (per unit width), i.e., the energy rate determined from the 
initial free stream stagnation enthalpy and the maximum 
mass rate (from free stream initial conditions) at which gas 
could be added externally to the base. 

It is noteworthy that these results are obtained without the 
use of empirical constants. However, for further application 
of Figs. 1 and 2, the numerical values of ¢ must be known. 
Although complete o data for supersonic mixing are not availa- 
ble, the temporary relationship 


[4] 


(based on data with M,, < 2) suggested by Korst and Tripp 
(3) may tentatively be applied until more experimental in- 
formation on supersonic mixing becomes available. (Meq is 
the value of the Mach number immediately after the Prandtl- 
Meyer expansion to the base pressure.) 

The theoretical trends indicated in Figs. 1 and 2 are in agree- 


= 12 + 2.758M-, 


ARS JouRNAL 


~ 
| = = 
= 
> 


ment with those obtained experimentally by heating the fore- The choice of X is constrained by the velocity equation 
body and base (4). y ae 
¢g=0= -- I, 
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peratures,”’ Proc. Fifth Midwestern Conf. on Fluid Mech., Univ. of Michigan, are 
April 1957, pp. 187-199. 
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Optimum Staging Techniques 7 _ where 7 is the Lagrangian multiplier. 
The solution to Equations [5] is 


LEON WEISBORD! 
ave T, sin J, — 


Curtiss-Wright Corp., Wood-Ridge, N. J. 1 + — 1) ik ~ 


ETHODS for the determination of the optimum staging — 61 
ratios for multistaged rockets were presented by Schur- Equations [6 and 4] must be solved for dy, As, .. . Ay by an 
mann (1)? and Weisbord (2). Although dealing with similar iterative method. 
conditions, there was a basic difference in the variables used. 
Schurmann minimized the gross missile weight for the case of 


Case 2 Variable ¢, I and rm 


different structural factors, and final thrust to weight ratios In many cases, it is more convenient to specify the initial 
among the stages, but the same specific impulse in every thrust to weight ratio ro rather than the final r;. For this 
stage. Weisbord handled the problem for the case of different situation, a minimization can be performed if Equation [3] is 
structural factors and specific impulses. rewritten in terms of 7 
Arens (3) commenting on (2) did not realize this distinction ” aa 
and stated that the solution of (2) is implicitly contained in In (An = 
Ton a 


(1). However, the solution offered in (3) is simply a solution —_ 
to the problem of (2) through the use of Lagrangian multi- 


The velocity Equation [4] using [7] becomes 
pliers. Moreover, this solution has already been presented 


by Vertregt (4) in 1956 (who used slightly different definitions V N Xi 

a g r=1 + €n(An 

Nevertheless, it might be worthwhile to derive solutions for ; 
the following more general cases: Different structural factors 7, sin Pn — | 
e, different specific impulses 7, and either different final thrust 
to weight ratios r, or initial thrust to weight ratios ro. (The Birman 2A oe ‘ 
definition of the structural factor € is the ratio of the gross Re-solving Equations [5] using Equation |8] as the new con- 
weight less propellant and payload weight of the stage to the straining relation, obtain 
gross weight less payload of the stage.) 
Wo, Won n 4 re rer = i, = constant [9] 


= 

Wo. — Won, : 

™ ii) For Case 2, a numerical solution can be obtained by an itera- 

tion between Equations [9 and 8]. 

_ ' Pili It should be noted that Equations [6 and 9] are general ex- 

Define the payload ratio of the nth stage A, as the ratio of pressions that contain certain more special cases. Forinstance, 
the gross weight of the nth stage to the gross — of the if 


(n + 1) stage 


Case | Variable ¢, J and r; 


r= gn 1 =e 


The overall payload ratio, W,,/ 1, which is to be minimized, 


and neither ro, nor r;, are considered, then 


is 
= =... = An 
= Mids,.... hy [2] [10] 
a is the criterion for minimization. Equation [10] will be 
ee eee ee recognized as the optimization criterion of Malina and 
The burning time of a stage can be expressed as Summerfield (5) 
On For the case where ro, and r;, are not considered, Equations 
At, = rp, [3] [6 and 9] can be rewritten as 
11] 
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Lastly, a conversion to the optimization criterion of (1) can 
be made by using the conversion . 


= at 12 


and setting J, = Jz =... =Ty. 

Mention should be made of the excellent articles on opti- 
mum staging techniques by Subotowicz (6) and Hall and Zam- 
belli (7). Both of these papers present weight minimization 
procedures for the case of variable specific impulses and struc- 
tural factors. 

It is hoped that this note will clarify the comment of Arens 
(3), and will be of some use to persons interested in missile 
weight minimization. 


Tn 


References 


1 Schurmann, E. E. H., “Optimum Staging Techniaue for Multistaged 
Rocket Vehicles,’ Jer Propusion, vol. 27, Aug. 1957, pp. 863-865. 

2 Weisbord, L., ‘“‘A Generalized Optimization Procedure for N-Staged 
Missiles,” Jer Propuusion, vol. 28, March 1958, pp. 164-167. 

3 Arens, M., “On a Generalized Optimization Procedure for N-Staged 
Missiles,”’ Jet PRoPpuULsION, vol. 28, Nov. 1958, p. 766. 

4 Vertregt, M., ‘‘A Method for Calculating the Mass Ratios of Step 
Rockets,”’ J. British Interplanet. Soc., vol. 15, March-April 1956, p. 95. 

5 Malina, F. J. and Summerfield, M., ‘‘The Problem of Escape from the 
Earth by Rocket,”’ J. Aeron. Sci., vol. 14, Aug. 1947, pp. 471-480. 

6 Subotowicz, M., ‘“‘The Optimization of the N-Step Rocket With 
Different Construction Parameters and Propellant Specific Impulses in Each 
Stage,”’ Jer Propu sion, vol. 28, July 1958, pp. 460-463. 

7 Hall, H. H. and Zambelli, E. D., “On the Optimization of Multistage 
Rockets,"’ Jer Propunsion, vol. 28, July 1958, pp. 463-465. 


Suppression of Instability 


in the Plungerjet Rocket Engine 


HERBERT SHIEBER' 
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ECENT issues of the ARS Journat have indicated a 
marked difference of opinion among researchers regarding | 
a mechanism for explaining combustion instability in rocket. 
engines (1 to 5).? 
tion expressed in (5) that “intellectual controversy is a strong 
stimulus to scientific progress,” some combustion instability 
phenomena experienced during the development of the 
Hughes Plungerjet XLR107-HT-1 monopropellant rocket 
engine are reported herein. These experimental phenomena, 
while unique to the Plungerjet engine geometry, pertain to 
the subject dispute in such a manner as to suggest possible 
experiments for resolving the argument, or at least for ampli- 
fying the contending combustion instability theories. 
In order to appreciate the pertinence of the Plungerjet ex- 


perimental experience to the controversy, it will be necessary | 


to comprehend the operation of the Plungerjet engine. Fig. 
1 is a cutaway view of this engine, showing its essential fea- 
tures. 
the propellant is a free 
by virtue of the area in- 


bustion chamber. This piston, 


equality existing on its liquid and gas sides, provides the pres-_ 


sure differential essential for propellant feed. As the pro- 
pellant is expended, the combustion chamber volume grows 
in the manner of an end-burning solid grain. 

It is thus seen that the Plungerjet engine traces the effect 
of chamber length on combustion instability during each 
firing. This is particularly relevant because much of the 
subject controversy appears to concern the relative suf- 
ficiencies of the contending theories to predict the 
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In the interest of supporting the tong 7 


It can be noted that the device employed to pressurize — 
piston, which also serves as the in- 
jector head separating the propellant reservoir from the com-— 


existence of a maximum chamber length beyond which the 
fundamental longitudinal mode of oscillation does not oc- 
cur. Peculiarly enough, the original configuration of the 
piston assembly (Fig. 2, Configuration A) was such as to in- 
hibit any combustion instability, irrespective of chamber 
length during operation. The choice of this configuration 
(piston rod vented into combustion chamber) was due at 
the time solely to structural considerations, that is, the pis- 
ton rod weight was minimized by reacting most of the propel- 
lant pressure on the piston rod with internal gas pressure 
However, in a later design of the engine, it was considered 
expedient to use the piston Configuration B shown in Fig. 2 
inasmuch as the saving in piston rod weight did not appear to 
justify the cost of sealing and insulating the interior of the 
rod. When Configuration B was static fired, however, the 
engine operation became a striking example of high-frequency 
instability. Although the length of the Plungerjet piston 
stroke was evidently not sufficient to demonstrate the afore- 
mentioned stability boundary effect of critical maximum 
chamber length, the exhibited instability otherwise conformed 


PISTON ROO 
‘SPACE 


Fig. 1 Cutaway view of Plungerjet rocket engine 
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Configuration A 
Original and Final Piston Assembly Design 


— 


Configuration B 
Intermediate Piston Assembly Design 


Fig. 2 Plungerjet piston assembly designs 
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very closely to that predicted by the theory and experiments 
reported in (1 and 2). That is, the frequency of the pressure 
oscillations showed dependence on instantaneous piston posi- 
tion (chamber length), and the factors affecting the onset 
and amplitude of the oscillations were likewise typical of 
longitudinal high-frequency combustion instability. 

The onset of instability appeared sooner and the amplitude 
of the oscillations were greater when the larger of two replace- 
able exhaust nozzle throats was installed in the engine. With 
the smaller of the two throats installed, there were practically 
no pressure oscillations, even at maximum chamber length. 
However, with the larger throat installed, the engine oper- 


Theoretical and Experimental Analysis" 
of a Cowling as a Means of a Drag 
Reduction for an Axisymmetric Center 
Body’ 

MARIAN VISICH Jr.2 and ANTHONY MARTELLUCCE 
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NTERFERENCE principles have been employed quite 
extensively for the reduction of the wave drag of various 


supersonic configurations. Linearized flow techniques have © 
been the major tool for the analyses conducted. However, in _ 
order to contain a sizable volume these techniques represent — 


ation was invariably unstable. 
A second factor that had a strong influence on the severity 


of the instability was the presence of inert matter in the pro- 
pellant. For considerations that are outside the scope of this 
Note, the propellant loading procedure included a thorough 
flushing of the propellant reservoir with Freon 12 gas. This 
gas is very soluble in the propellant, and on occasions when 
the engine was defueled and the propellant subsequently re- 
used after another exposure to the Freon flushing procedure, 
the ensuing instability was appreciably increased. 

The magnitude of the chamber pressure (which can be 
preadjusted in the Plungerjet) also seemed to affect the onset 
and amplitude of the combustion instability, with the higher 
chamber pressures appearing to create somewhat more 
severe pressure oscillations than did the lower chamber pres- 
sures. This effect, however, was not so pronounced as the 
aforementioned effects. 

Needless to say, this experience with the destabilizing in- 
fluence of piston Configuration B prompted a design return 
to the original piston Configuration A. From a developmen- 
tal point of view, this engine’s combustion instability tend- 
ency has now been effectively suppressed. From a more 
fundamental viewpoint, however, several interesting ques- 
tions remain, and these questions may have some bearing on 
the subject controversy concerning combustion instability 
explanatory mechanisms. 

For instance, is the stabilizing influence of the open piston 
rod (Configuration A) explained by the fact that the cavity 
represents enough additional chamber length to exceed the 
critical chamber length? There are objections to this ex- 
planation. In the first place, the open piston rod is forward 
of the injectors and, as such, should provide a different phase 
relationship between the pressure waves and the local burning 
rates than would be the case in the conventional combustion 
model with an equivalent increase in chamber length applied 
downstream of the injectors. Secondly, the open piston rod 
represents only a small portion of the cross-sectional area of 
the chamber, further departing from geometric equivalence. 

It would be most interesting if the research teams investi- 
gating combustion instability could correlate this cavity effect 
into their combustion theories. Specifically, is there a cor- 
relation between critical combustion chamber length and 
cavity length? Is the placement of the cavity forward of the 
injectors significant, or could it be placed at the exhaust 
nozzle end, or even radially disposed? Is there a cavity 
length that would reinforce rather than attenuate the pres- 
sure oscillations? Answers to these questions would not 
only illuminate the subject dispute, but might also provide 
an effective design tool for suppressing combustion instability. 
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slender bodies with relatively large fineness ratios. 

The present investigation employs nonlinear techniques for — 
the drag reduction of a system that contains a large volume _ 
per unit length. The configuration consists of a boat-tailed 
body of revolution surrounded by a large diameter ring which 
effectively reduces the wave or thickness drag of the center 
body. This reduction is accomplished by reflecting the coni- 
cal shock generated from the apex of the center body on to the 
boat-tailed rear part. Thus, the static pressure on the rear 
part of the body is increased, and the pressure drag considera- 
bly reduced in comparison to that of the center body alone. 
This type of drag reduction has been investigated for various 
configurations of body and ring by both the method of charac- 
teristics (1) and by linearized techniques (2). In this report, 
however, a finite-strength conical shock from both the nose of 
the center body and the lip of the cowling ring is considered. 
The result is only verified for a particular configuration at a 
design point Mach number of 3.09. It should be noted that a 
vehicle of this type is efficient in a limited range of Mach 
number close to the design value. 

A theoretical investigation has been conducted in order to 
determine a ring-body combination that will demonstrate 
the use of this method as a means of drag reduction. The 
body was designed for a Mach number of 3.09. A body with 
a large volume per unit length was desired; thus nonlinear 
techniques were required in the design. The center body 
length, leading edge semi-apex angle and maximum diameter 
were selected somewhat arbitrarily. In addition, a cuspated 
forebody was selected so that a greater part of the compression 
would be isentropic. The method of characteristics for 
axisymmetric irrotational flow was utilized to design the body. 
The generalized equations employed are given in (3). From 
the apex to the position of maximum thickness, the body 
shape was selected in a manner that would insure a gradual 
isentropic compression and an ensuing expansion which 
would avoid separation. With the forebody selected, the 
remaining body shape was established. 

The minimum diameter ring and its approximate axial 
position were then selected so that the area of the annulus, 
between the center body and ring, would permit supersonic 
flow to be established. An abbreviated characteristics analy- 
sis was then conducted to place the ring axially relative to 
the center body. 

The rear section of the center body was then calculated 
with the cowling ring in place. This calculation was more 
complicated due to the intersection of the tip shock with the 


Theoretical Investigation 
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Characteristics of the center body with the ring 
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shock emanating from the ring. The shape of the reflected 
shock was determined point-by-point, since conditions up- 
stream were known. From this shock, once again the method 
of characteristics for irrotational flow was used to deter- 
mine the body shape. This is permissible, since the varia- 
tions in entropy are small, and in a first approximation may 
be neglected while sufficient accuracy is maintained. The 
characteristic net for the center body with the cowling ring is 
shown in Fig. 1. An assembly view of the cowling model is 
shown in Fig. 2. 


Spark shadowgraph of 


7% 


the cowling model 


Fig. 2 Assembly drawing of cowling model 


900 


Experimental Apparatus and Testing Technique 

The experimental investigation of the body of revolution, 
surrounded by a cowling or ring, was conducted in the 10 X 
10-in. supersonic blowdown wind tunnel of the Polytechnic 
Institute of Brooklyn Aerodynamics Laboratory. Tests were 
conducted at a stagnation pressure of 200 psia with a corres- 
ponding test section Reynolds number per foot of 3.37 X 10’. 
A typical spark shadowgraph of the cowling model is shown 
in Fig. 3. 

The pressure drag was obtained by measuring the pressure 
distribution over the body and integrating it over the drag 
surface on which it acts. Tests were made of the center body 
without a ring and of the center body with the cowling. The 
experimental Cp distribution for the center body without a 
ring is shown in Fig. 4. Also shown for comparison is the 
distribution obtained by the characteristics analysis. The ex- 
perimental Cp distribution for the ring configuration is shown 


in Fig. 5. 
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Fig. 4 Comparison of theoretical and experimental pressure 
coefficient distribution for the center body alone 
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Volume = 31.068 Length = 9.000 as 
Estimate of 
Pressure Base Total skin friction 
Wetted drag, drag, drag, drag, 
center body 7.07 206 550 2.573 (0.2% 2797 0.194 
cowling body 7.07 6.2.06 172.0 of 0.702 0.18 0.887 0.605 
semi-apex angle 10.36 52075 1.040 2.115 
Haack-Adams 6.73 6.61. 58.2 0.697 0.661 1.358 
| TR 1271 6.73 6.61 58.2 0.560 0.661 1.221 


Table 1 Comparison of drag components for bodies of equal volume and length en 


Method of Data Reduction 
7 

The pressure drag of the center body was determined from 
a numerical integration of the static pressures over the surface, 
utilizing the experimental values of the recorded pressures. 
The pressure drag on the ring and struts was obtained by 
measuring the pressures on the respective surfaces. 

The base pressure on the center body was also measured, and 
the value of the drag coefficient computed. These results are 
of the same order of magnitude as the values obtained by 
Chapman (4) and Love (5). Since the pressure distribution 
along the body was known, using the method outlined in (6), 
an estimate of the skin friction was made. For this calcula- 
tion, it was assumed that the boundary layer was entirely 
turbulent and the wall temperature constant. 


Results and Conclusions 


For the center body without the ring, the pressure drag co- 
efficient based on the maximum diameter is 0.364, while for 
the configuration tested, the value for the center body alone 
is 0.060. This corresponds to a reduction of 83.5 per cent of 
the pressure drag of the center body. Comparison of the 
total drag of the cowling configuration (pressure drag of the 
entire system and skin friction on the ring and struts) to the 
total drag of the center body without the ring indicates a re- 
duction of 50.0 per cent. 

A comparison of the drag of the present cowling system 
with various linear systems of given length and volume cannot 
be achieved because of the relatively large volume per unit 
ength. Therefore, the system is compared to various bodies 
of revolution to further evaluate the reduced drag. Since the 
bodies considered have different maximum cross-sectional 
areas, the drag coefficients based on this reference area cannot 
be compared, but values of the product of the drag coefficient 
and the body maximum cross-sectional area are considered. 


- | O EXPERIMENTAL POINT | \ 


Fig. 5 Experimental pressure coefficient distribution of a body 
of revolution surrounded by a cowling ring 


JUNE 1959 


The results are illustrated in Table 1. Comparison of the | 
cowling system to the optimized Haack-Adams (7) body tends — 
to indicate that the drag reduction is appreciable. One must — 
note that the cowling body is capable of producing a substan- 
tial L/D. 

Thus the cowling center body presents a useful and prac-— 
tical system for employment as a flight vehicle with a circular | 
wing. In addition, the annulus contains air possessing a high — 
pressure recovery; thus a system of this type could also be — 
used as an auxiliary ramjet engine or as a bypass for a ramjet | 
engine. 
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An Optimum Transfer Path ie an 
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Cireular Orbit 
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HE PROBLEM considered in this paper is an optimum 

transfer path from the perigee of an elliptic orbit to any 
point on a higher energy circular orbit. The optimum con- 
sidered is the minimization of characteristic velocity, or the 
amount of fuel consumed, using two impulses. The time of 
transfer is left open. Lawden (1)? has treated this problem, 
and similar more general problems, exceptionally well. 

A rigorous mathematical derivation is presented, proving 
an optimum path to follow is the Hohmann (2) transfer el- 
lipse. 

Introduce a hypothetical circular orbit (B) (see Fig. 1) of 
radius r,. Such an orbit has a smaller major axis than the 
Received March 2, 1959. 
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Fig. 1 Orbits in a central force field 


major axis of elliptic orbit (A) and is therefore a lower energy 
orbit. At the point of tangency P of orbits (A) and (B), a 
space vehicle in orbit (A) or in orbit (B) would have the 
same potential energy. Therefore a space vehicle at point P 
in orbit (B) would have less linear velocity than a space 
vehicle traveling in orbit (A) at point P. 

Consider a space vehicle traveling in elliptic orbit (A). 
Suddenly an impulse i is applied at P, giving rise to a velocity 


AV, pi 


increment, AV, » (see Fig. 2). The magnitude is 
such that energy is added to orbit (A). Another way to 
describe this energy addition is a 
V, + AV,| 2 \V> [1] 
For such a AV, there exists a AV,’ (see Fig. 2) such that 
V, + AV, = V. + AV.’ [2] 


A comparison is next made between the transfer from orbit 
(A) to orbit (D), with the transfer from orbit (B) to orbit (D). 
Consider the following three characteristic velocities: 

Ellipse-Circle 


av, + 


2 Circle-Circle 


3 Hohmann Optimum 
Circle-Circle 


+ 


[aan 


1, the final correcting velocities AF appearing in 
1 and 2 are the same. 


From Fig. 2 we observe oe 
— | |-—— 
Bi 


“abe! eri 


«Fig. 2. Diagram of vector velocities 


Since the Hohmann circle-to-circle transfer ellipse gives an 
optimum path, we have 


+ < lave! + | 


From agi [3 and 4] 


lar 5) 


From Equation [5] 


opt| Feo! 


j——>» | lar. 
AV, + iA lar > [Av AV en + AF opt [6] 
The Hohmann characteristic velocity AV + AF on 


provides a lower bound for the class of characteristic velocities 
in the circle-to-circle transfer. Likewise the right-hand side 
of Equation [6] provides a lower bound for the class of 
characteristic velocities in the ellipse-to-circle transfer, since 
the lower bound given by Equation [6] can actually be 
achieved by taking the Hohmann transfer ellipse. This 
Hohmann transfer ellipse would be identical to the optimum 
transfer path taken from circular orbit (B) to circular orbit 
(D). 

Since the greatest lower bound for the class of characteristic 
velocities can be achieved, it follows 


AY, | | + laren! — 


Nomenclature 


_ = distance of perigee from focus 
AV, = velocity increment at perigee P in elliptic orbit (A) 
Ky = velocity at perigee in elliptic orbit (A) 
AV.’ = velocity increment at P in circular orbit (B) 
AF = velocity increment correcting final velocity of transfer 
path to circular velocity 
AV»: = Hohmann optimum velocity increment at P in circular 


orbit (B) 


q 
P 
[4] 
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42 
] 
450 
= 


—* So far, the two-dimensional model has been programmed. 
In this, it is assumed that the magnetic and rotational equa- 
tors are the same. The effect on charged satellites of various 


path to circular velocity 


V. = velocity in circular orbit (B) masses and charge have been investigated for an initially cireu- 
. ae 7 7 lar orbit. A charge-discharge procedure has also been studied 
y = difference between v. p and V. tg determine the feasibility of using this to control and change 
; : At present the general three-dimensional equations are 
The author expresses his gratitude to E. Beltrami of the being programmed. These include the angle between the ro- 
Operations Research Group, for valuable discussions con- tational and magnetic axes, the Earth’s rotation, and 
cerning the problem. Also to W. Lindorfer of the Dynamics charged body moving with any initial velocity. 


Research Group, for his excellent mathematical criticisms. 
Two-Dimensional Equations of Motion for a Charged Body _ 
References in the Magnetic and Gravitational Field 
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a —-GM.M + hQ6 
Electrically Charged Bodies Moving in = — 


the Earth’s Magnetic Field hQi 


= (276 + rb)M, 
W. W. FAIN! and B. J. GREER? 


Chance Vought Aircraft, Inc., Dallas, Texas : where 


r = distance from the Earth’s center to the center of the- 
N INVESTIGATION has been initiated to determine charged body 
the effect on the path of electrically charged bodies M, = mass of the charged body 
moving in the Earth’s magnetic field. The study is applicable Mr = mass of the Earth 7 a 
to h = strength of the dipole a 
1 Perturbation in satellite orbits (charging the satellite Q = charge on the body Ae 7 
can be used to afford a small measure of control over the orbit). 6 = angular distance measured from an arbitrary reference 
2 The effect on the path of a rocket fired from the Earth. G = gravitational constant 
3 The effect on charged particles entering the Earth’s oA ‘ ; 
magnetic field. (In particular it is hoped that the shape of Results for Two-Dimensional Model 
the radiation belt around the Earth can be deduced from The body was given sufficient angular velocity fora circular | 
this. ) orbit at a distance 6.5 X 10° km from the center of the model - 


" Earth, and the charge was “turned on.’”’ In Fig. 1 the dif-— 
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2 Associate Aerodynamics Design Engineer. as a function of the angular distance from the ‘turn-on’ 
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Ist Revolution 


5 th Revolution 


\ 39th Revolution 


Fig. 1 Variations from initially circular orbit as function of angle 
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Rmax- Ro * 2.58 X 10° Cm. 


Fig.2 Example of path followed by small mass with large charge Fig.3 Example of path followed by satellite with mass 10 times 
that in Fig. 2 


point. It can be seen that the maximum variation in distance 
from the initial orbit is approximately 20 km for a body of this 
mass and charge. With continued revolution of the satellite 
the peak distance shifts away from its origina] position, and 
the peak reduces in amplitude. At the same time, the average 
distance from the Earth’s center is increasing, and the indica- 


tions are that a new circular orbit is approached asymptoti- * Ms = 500gm 
cally. Q = 1073 Abcoulombs 


In Fig. 2. an example is shown in which the acceleration af- 
forded by the magnetic field on the moving charged body is 
large compared to the gravitational acceleration. The looping 
is not surprising, since a charged body originally in straight 


Charge turned on and off 


MAX. CHANGE IN ORBIT OF 
SATELLITE AS A FUNCTION \ 


OF MASS AND CHARGE 
R, * 6,500 Km FROM 
EARTH CENTER 


0 


Ry -Ro IN Km. 


1o¢ 


MASS OF SATELLITE IN GRAMS 


Fig. 4 Maximum variation in orbit for the first revolution as a Fig. 5 Effect of charge-discharge procedure on the orbit | 
function of mass and charge 
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line motion takes up a circular path in a magnetic field (cyclo- 
tron, for example). It will be noticed that the maximum dis- 
tance approaches twice the radius of the original orbit. Twice 
the distance is the upper limit, since at this distance all of the 
kinetic energy is converted to potential energy. It will be 
noted that the orbit does not necessarily repeat itself upon 
subsequent revolution. 

In Fig. 3 the orbit is shown for a satellite of 10 times the 
mass of the satellite in Fig. 2. It does not loop, and the maxi- 
mum variation in orbit from the original circular orbit is 2580 
km. 

In Fig. 4 the maximum variation in orbit (for the first 
revolution) is plotted as a function of mass of the satellite and 
ihe charge in abcoulombs. The dotted line at the top is the 
upper limit at twice the initial orbit distance. 

In Fig. 5 the charge was repeatedly turned on at the 
ninimum distance from the original circular orbit and turned 
ff close to the maximum distance. This tends to ‘‘walk” the 
satellite away from its original orbit. (It should be noted that 
the scale is considerably distorted since zero is at 6500 km 
from Earth center—the shaded area is actually very large 
-ompared to the variation.) 


The practical use for this method of satellite control is 
probably very limited, since for a realizable charge (10~* to 
10-2 abeoulomb is probably optimistic for reasonable radii) 
only small masses can be moved significant distances. It is 
possible that it could be useful for fine control (several km) 
and possibly for transfer of small masses from one orbit to 
another (possibly small missiles). 

The three-dimensional case will give further information, 
particularly as concerns rotating the plane of the orbit. 


(pplication 


On the Use of Side-Jets as Control 
Devices 


H. P. LIEPMAN'! 


University of Michigan, Ann Arbor, Mich. 


Wind tunnel experiments with side-jets, issuing laterally 
near the base of slender bodies in a supersonic stream, have 
suggested the existence of a sizable and usable interaction. 
With this interaction force, the use of jet reaction controls 
may be as attractive for flight within the atmosphere as it 
obviously is for flight outside the atmosphere. This note 
indicates the altitude regime of interest and the order of 
magnitude of the interaction bonus for a lateral control 
jet located near the base of a body of revolution. 


VER since flight beyond the Earth’s atmosphere has been 
under consideration, the value of jet reactions for control 
of attitude and flight path has been obvious. Furthermore, 
jet reactions may also provide useful control forces during 
atmospheric flight. In this case, the interaction of the gas jet 
or its component with the external flow field must be taken 
into account. 
A modest and exploratory study? of such interaction effects 
has been under way at the Aeronautical Engineering Labora- 
tories of the University of Michigan since 1950 (1 to 3). A 
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typical Schlieren picture from this program is shown in Fig. 1 
showing an air jet from a sonic orifice issuing into a Mach 
3.9 stream. Other test parameters are (d/D) = 0.08 and 
(po;/p) = 414. 

Recent results (2, 3) of tests with side-jets near the base of 
a missile indicate that considerable advantage may be gained 
by utilizing the interaction effect between a side-jet and super- 
sonic main stream as an additional control force contribution. 
This note summarizes the magnitude of the interaction effect 
and the altitude regime over which it may be useful. 


Interaction Effect 


Experimental] data and theoretical studies (2, 3) show that 
the interaction of a lateral jet near the base of a body of revolu- 
tion with supersonic flow leads to a significant normal force 
increment (N,) which is additive to the normal jet reaction 
(N ,) experienced in a vacuum where aerodynamic interactions 
do not exist. Side-jets located forward of the base, such as 
near the nose or the midpoint, do not produce significant addi- 
tive interaction forces (1, 2). 

A simple theoretical model (3) of the interaction effects of 
a two-dimensional side-jet indicates that the ratio of interac- 
tion force to vacuum-jet reaction (N,/N,,) depends on the jet 
stagnation to ambient pressure ratio (po,;/p), and on the 
specific heat ratio of the jet gas. Furthermore, theory as well 
as experiments to date have shown no significant dependence 
of the interaction force on the type of boundary layer, i.e., 
laminar or turbulent; hence, Reynolds number effects seem 
negligible. 

It has been found convenient (2) to express the experimental 
and theoretical results in terms of a normal force magnifica- 
tion factor K = 1 + (N,/N,) which is simply the ratio of the 
vacuum-plus-interaction normal force to the normal force in 
avacuum. Test results, using room temperature air exhaust- 
ing through a sonic orifice at three free stream Mach numbers, 
are given in Fig. 2 in terms of the normal force magnification 
parameter K/(. + 1) vs. the pressure ratio (py;/p). The 
denominator (7 + 1) was chosen in an attempt to eliminate 
the Mach number effects from this discussion. It is seen in 
Fig. 2 that this Mach number correlation is only fair, but use- 
ful for the subsequent discussion, and probably acceptable for 
Mach numbers from 2 to 5. 

From the curve of Fig. 2, some impressive values for the 
magnification of the vacuum normal force due to interaction 
effects can be predicted. Representative normal force magni- 
fication factors K are 


(po,/p) 50 100 500 
M 
1.50 1.35 110 
5 2.25 2.02 1.66 


Fig. 1 Schlieren picture of sonic side-jet issuing from ogive- 
cylinder model into M = 3.9 stream. (po;/p) = 414, (d/D) = 
0.08, exposure time }4o sec 
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Fig. 2 Normal force magnification parameter vs. jet stagnation 


pressure ratio 


10 

—— ALTITUDE IN THOUSANDS OF FEET 

Fig.3 Normal force magnification parameter vs. altitude and jet 
stagnation pressure 


Side-Jet Control Effectiveness 


It is clear that magnifications of jet reactions due to interac- 
tion effects can be used to advantage either to increase the 
effectiveness of reaction controls or to decrease the amount of 
control thrust needed. For convenience in predicting side-jet 
effectiveness, the normal force magnification parameter oi 
Fig. 2 has been plotted against altitude for various jet stagna- 
tion pressures in Fig. 3. The ICAO atmosphere was used t 
relate altitude and ambient pressure. It is seen that the in- 
teraction magnification is most pronounced at low altitudes 
However, at high Mach numbers, say 5, the magnification is 
still quite significant at altitudes as high as 150,000 ft. 

Fig. 4 has been prepared to show the resultant control force 
of a nominal 100-lb control jet with a stagnation pressure of 
500 psia as a function of altitude at flight Mach numbers of 3 
4 and 5. It is obvious that the interaction bonus can be 
traded for a sizable reduction in fuel, or size and weight of the 
control jet. In specific applications, the use of control jets 
within the atmosphere may be entirely competitive witl 
aerodynamic control devices. 


Conclusions 


Some results of a continuing investigation into the interac- 
tions between a side-jet and the supersonic flow field of a 
slender vehicle show that there exists a helpful interaction 
force which considerably magnifies the vacuum reaction of a 
control jet near the base. The magnification decreases with 
decreasing free stream Mach number and with increasing ratio 
of jet stagnation to ambient static pressure. 

For practical applications, the helpful interaction force is 
most significant at altitudes up to about 100,000 ft depending 
on Mach number and jet stagnation pressure. Since reaction 
controls are needed for flight beyond the atmosphere, their 
practical use within the atmosphere may, therefore, be quite 
competitive with conventional aerodynamic means of control. 


d = jet exit diameter ene 

K = normal force magnification factor = 1 + (N;/N>») 

M = free stream Mach number 

p = ambient pressure 

Poy = jet stagnation pressure 

N; = interaction normal force 
N, = normal force due to jet exhausting into vacuum 
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Technical Comments 


Remark on Thermodynamics of AIO, 
JOHN S. GORDON! Aa 


cone 
Thiokol Chemical Corp., Denville, N. 


HE IDEALIZED treatment of inert liquid Al,O; recom- 

mended at temperatures below 3000 K in a recent contri- 
bution (1)? will, in most cases of current interest, lead to 
pproximate and optimistic specific impulse values. 

The work of Brewer and Searcy (2) has indicated that the 
ignificant gaseous species occurring in the volatilization of 
are Al, Oo, O, AlO, and possibly Al,O2, but not 
Ale and are usua]ly omitted from propellant per- 
ormance calculations since no molecular constant data have 
een deduced for them. It is important, however, to allow 
or the interaction of aluminum-oxide vapors and halogens, 
ince the Al—CI bond energy is about 105 keal, as compared 
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Table 1 Thermodynamic functions of AlO, cal/mole, K 
-(F- H- 2698 16, 
T,K Ey)/T E/T s° keal/mole 
400 47.213 7.182 54.395 7.767 0.77186 
800 52.365 7.715 60.081 8.568 4.0714 
1200 55.563 8.050 63.613 8.829 7.5589 
1600 57.910 8.261 66.172 8.951 
2000 59.770 8.407 68.178 9.028 14.714 
2400 =8.516 69.829 9.086 18.337 
2800 62.632 8.601 71.234 9.1388 21.982 
3200 8.671 72.458 9.191 25.648 
3600 64.811 8.732 73.544 9.250 29.336 
4000 65.734 8.787 74.522 9.317 33.049 
4400 66.574 8.839 75.413 9.394 36.792 
4800 67.345 8.889 76.234 9.479 40.566 
5200 68.059 8.938 76.997 9.571 44.376 
5600 68.723 8.986 77.710 9.668 48.224 
6000 69.344 9.035 78.380 9.767 52.112 


Tabie 2 Thermodynamic functions of AICI, cal/mole, K 
—(F H - H 298.16, 
T,K £E/T s° keal/mole 
400 49 221 7718 56.940 8.560 0.85917 
800 54.764 8.252 63.016 8.917 4.3734 
1200 58.161 8.494 66.655 9.029 7.9654 
1600 60.625 8.637 69.263 9.100 11.592 
2000 62.564 8.736 71.300 9.159 15.244 
2400 64.164 8.811 72.975 9.212 18.918 
2800 65.527 8.872 74.399 9.263 22.614 
3200 66.715 8.924 75.639 9.313 26.329 
3600 67.769 8.970 76.739 9.362 30.064 
4000 68.716 9.011 77.728 9.411 33.819 
4400 69.577 9.050 78.627 9.460 37.593 
4800 70.366 9.086 79.453 9.512 41.388 
5200 71.095 9.121 80.216 9.565 45.203 
5600 = 71.772 80.927 9.622 49.041 
6000 72.405 81.593 59.684 52.902 
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to 119 keal for Al = 0 (3). Nor must AJN be ruled out _ 
a possible combustion product, since its heat of formation in 
the crystalline state is —57.7 kcal/mole (4). Theoretical per- a 
formance calculations have shown large amounts of AICI, AICI. _ 
and AlCl; gases at temperatures of 2500 to 3500 K, usually 
exceeding the amounts of AlO, Al,O, Al and AJH present. Un- 
fortunately, nothing whatev er is presently known about 
OAIC! gas which may be the most important species of all. _ 

We have recently recomputed the thermodynamic functions _ 
of AlO and AICI gases, as well as many other molecules, — 
from Herzberg’s data (5). These data are necessary in pro-— 
pellant performance calculations, and are shown in Tables 1 
and 2. Data have also been compiled for AlCl; and AICl;. 

More recent values for H {29s and S2ss of Al,O; are available | 
(6,7). 
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Comment on “Performance Calculations 
for Hybrid Nuclear-Chemical Rocket 
Propulsion Systems” 


4, 
JOHN S. GORDON! = 


any Thiokol Chemical Corp., Denville, N. J. 


Thermodynamics of Al2O3,'’ Jer Propuusron, vol. 28, 


HE PERFORMANCE estimates by Green and Carter (1)? 

for metal-H,O systems are subject to revision when more 
reliable data on light-metal combustion products become 
available. Thermodynamic functions for LixO, LiOH, LiO 
and other lithium species were generated by the author during © 
January 1957 and have been made available to Aerojet-Gen- 
eral Corp. personnel on two separate occasions; there is no 
evidence in the paper by Green and Carter that these data 
were considered. The thermodynamic functions for Li,O (g) 
are included here (Table 1), since this appears to be the only 
Li species considered in their paper. 

The sublimation of LizO has been studied by Brewer and 
Margrave(2) and Blue, Margrave, Berkowitz and Chupka (8). 
It has been demonstrated that Li,O vaporizes mainly to the 
elements, with minor amounts of Li,O and LiO, and that 


Li,O (s) ~ LiO (g) + 104 + 5 kcal/mole 
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ay 
ae Table 1 Thermodynamic functions of Li,O (gas) 

Keal/mole —(F — E.)/T -&,)/T C,’, eal/mole, K 
0.9363 48 23 8.396 56.63 
ahs 54.46 9.709 64 17 
Le 10.37 58.59 0.66 69.26 

15.65 61.75 11.29 73.05 

12.27 80.64 
42.96 12.60 | — 84.10 
-48.48 

76.17 13.09 91.17 


LiO (g) > Li (g) + 0 (g) + 83 kcal/mole 
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R. W. FANSELAU! 


N A RECENT article? Rao maximizes exhaust nozzle exit 
_# momentum by formulating an extremal problem involving 
two constraints: Nozzle length and propellant mass flow. 
The extremal problem is formulated with a noncharacteristic 
line as the control surface. An investigation of the extremal 
problem together with the assumptions adhered to in its 
formulation indicate that the last left-running characteristic 
line is assumed for the control surface. Subsequent verifica- 
tion of this assumption is treated therein as an independent 
result which emerges from the first variation and exclusively 
characterizes the solution sought. Even though the control 
surface of the maximum solution is recognized by Rao as the 
last left-running characteristic line, there appears to be an 
error in the suggested use of the Euler-Lagrange constants 
to obtain solutions for any nozzle length. It is, therefore, 
predicted that for a given ambient pressure, the nozzle con- 
tours obtained therein yield integrated exit momentum values 
other than the maximum. Comments concerning the fore- 
going extremal problem and its solution are given in the fol- 
lowing paragraphs. 

The extremal problem is formulated for a given pressure 
ratio. Thus, the formulation can represent an investigation 


for either a singly or a doubly infinite family of nozzle con- 

tours. The exit momentum J, the propellant mass flow M 

and the nozzle length L are originally expressed by the follow- 
Received April 20, 1959. 
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Comment on “Exhaust Nozzle Contour 


ing line along some control surface 


J gi(P, W 6, y)dy {la 
E 
g2(P, 6, y)dy =(C, [2a 
E 
Lo f = [3a 


Equations [la, 2a and 3a] represent the formulation for a 
doubly infinite family of nozzles; both isentropic and nonisen- 
tropic nozzles are investigated. The control surfaces ad- 
missible in the extremal problem are left-running character- 
istic lines of two types 


dO = sin 6 sin @ 
dy dyM{1 + [(y — 1)/2]M*} ysin @ 
and 
dé dM (M2? — sin 6 sin 
dy {1 + [(y — 1)/2]M2} y sin (5 
ds sin? cos @ 
dn yR sin 


Thus far, the formulation is quite general. However, a 
considerable loss of generality results by introducing the isen- 
tropic relations 


y¥-1 — — 
P = Py E + =) M? (isentropic perfect) 
— 
p = po [1 + (isentropic perfect) 
y-1 
(adiabatic perfect) 


Through the preceding relations, only the control surfaces 
satisfying Equation [4] are admissible. For this reason, 
Equation [4] need not be entered as an additional constraint. 
Rao introduces the preceding isentropic relations into Equa- 
tions [la, 2a and 3a] to obtain for J, M and L, respectively 


E 
Ai(M, 9, y)dy [1b] 
= r= 6, y)dy C2 [2b ] 
L= f. fs(p)dy = Cs [3b 


Thus, the doubly infinite family of nozzles is reduced to a 
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‘ 
singly infinite family of isentropic types. Fortunately, this where 
s of generality does not invalidate the solutions of the ex- en Be 
tremal problem, as the condition for a maximum is satisfied G = Gy, M, 6, ¢) 
by nozzles from this isentropic family. The Lagrange con- F: = F2(y, M, 9, ) 
stants F; = F;(y, M, 0, 
Solutions of the foregoing extremal problem are obtained 
a in the following manner. For a particular ambient pressure’ 
- P,, one chooses the Mack number Me, at the nozzle exit. a 
—); = ypW? sin? 6 tan (@ — 8) [7] The flow angle 6z at the nozzle exit is then obtained from 
so, from the first variation = 
in 20 cot {10} 
afi fs _ afr Wet 
oM 00 OM 
i 6 The flow parameters M and @ along the last left-running char-_ 
acteristic line are obtained by solving Equation [4] to gether 
Ofi Of: Oe _ with 
oM aM 
= de 6+a=¢ 
= for values of and 6, using Y/YE as a parameter. 
the sufficiency condition is satisfied everywhere along the obtains the flow parameters M and @ along the same control =| 
ntropic left-running charae ter istic line by the J Jae obian surface directly from Equations [6 and rae Solutions of the — 


extremal problem obtained in this way will satisfy the Euler — 
equation, but do not necessarily satisfy the sufficiency condi-— 
tion given by Equation [9]. To insure the sufficiency condi- — 
tion, the solutions of the foregoing extremal problem must 
also satisfy Equation [4], which can be derived by differ- 
entiating the Lagrange constants with respect to the inde-. 
pendent variable. The location of the nozzle kernel design 
point and the construction of the nozzle contour follow as — 
discussed therein. 
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What Lockheed is doing today to develop 
Tomorrow’s missiles and spacecraft 


The world’s first polar-orbiting satellites...the 
world’s fastest ramjet target drone...a “fuel cell” that 
produces auxiliary power for spacecraft in a radically dif- 
ferent way...a TV camera-transmitter—the world’s 
smallest —that gives engineers and scientists on the ground 
a televised report of what takes place in test missiles dur- 
ing flight. These are just a few of the many activities and 
achievements of Lockheed’s Missiles and Space Division. 

Lockheed is System Manager and Prime Contractor 
of the ARPA DISCOVERER satellite series—and the U. S. 
Navy’s POLARIS missile. Both are highest priority pro- 
grams of the U. S. Government. And both are ahead of 
schedule. 

A Lockheed MSD work-force of over 16,000—includ- 
ing 3,500 scientists and engineers—is engaged in all 
phases of missile and space technology: satellite systems 
development; space communications; electronics; ionic, 
nuclear, and solar propulsion; magnetohydrodynamics; 
computer research and development; flight sciences; 
materials and processes; human engineering; electromag- 
netic wave propagation and radiation—and many other 
advanced fields. 

From these efforts at Lockheed will come many sig- 
nificant breakthroughs and scientific “firsts’—to speed 
the development of tomorrow’s missiles and spacecraft. 


Ramjet-powered supersonic target, Lockheed Q-5 kinc- First polar-orbiting satellites, the Lockheed-built 
FISHER electronically simulates enemy attackers—is being DISCOVERER I and II were orbited in two successive 
used by U. S. Army to evaluate and sharpen our nation’s launches. Lockheed is System Manager for this ARPA 
missile marksmanship, increase our defense capabilities. advanced scientific research program. 
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Transmitting performance data from missiles and test vehicles, 


direct to ground stations, Lockheed-developed TV camera- 
transmitter and PAM-FM telemetry systems are acknowledged to 
be the world’s smallest, world’s best. 


Revolutionary “fuel cell,”” under development at Lockheed, con- _Erupting from beneath the sea, a full-scale POLARIS 

verts chemical energy directly into electrical power—for auxiliary test vehicle demonstrates how operational missile will be 

power requirements of miniaturized spacecraft instrumentation — launched from submerged U.S. Navy subs. Ultimate range _ 
of the PoLaRIs will be 1500 nautical miles. 
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rate gyros (2,868,023). C. R. 
Bonnell, Columbia Heights, Minn., assig- 
— to Minneapolis-Honey well Regulator 


Liquid disposed within an enclosure so 
that the unfloated mass attached to either 
side of a flexible diaphragm is equal, per- 
mitting displacement of a hollow sealed 
housing about two axes at right angles to 
each other. 

Anti-backlash gear (2,868,028 ). 
J. C. Ziegler, Maywood, N. J., assignor to 
Bendix Aviation Corp. 

Gear arranged to drives in either direc- 
tion through a clutch and two associated 
gear units. A friction member associated 
with one unit applies a relatively low 
torque to the other unit, and a torque less 
than the torque to the first unit, through 
the driven gear. 

Rocket motor (2,868,127). H. M. Fox, 
Bartlesville, Okla., assignor to Phillips 
Petroleum Co. 

Piston disposed within a fuel tank to 
form a movable wall. The piston head 
within an oxidizer tank forms a movable 
wall of the tank and combustion chamber. 
A pressure sensitive means seals a hyper- 
golic fuel within the hollow piston shaft. 
A resilient member ruptures the sealing 
means, moving the piston into the oxidizer 
tank and into the fuel tank when released 
by a starting means. 

Propulsion of vehicles (2,869,479). P. 
Hutchinson, Ruislip, Middlesex, England. 


System for propelling a vehicle along a 
track, using gas in jet form directed into 


ducts, fixed to the track. Vehicle ducts 
register with successive track duct inlets 
and receive jet flow from track duct out- 
lets as the vehicle travels along the track. 
Convertiplane with tiltable cylindrical 
wing (2,868,476). IE. W. Schlieben, Con- 
cord, Mass. 

Means for propelling a wing comprising 
: circular body of toric section. A fuselage 
> mounted in the wing in hinged relation- 
ship. 
Rocket control (2,868,478). 
Cloughy, Morris Plains, N.J 

Steering device comprising a tubular 
housing in which a lining of refractory 
material is rotatably mounted. The hous- 
ing extends in a space-surrounding relation 
to the jet stream. Pivoting the housing 
causes part of the housing to engage the 
perimeter of the jet stream. 


T. Me 


Epitor’s Nore: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies of 
patents may be obtained from the Com- 
missioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 
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Friction type fuse (2,869,464). J. G. 
Villepigue and E. Schmued, Los Angeles, 
Calif., assignors to Northrop Aircraft, 
Inc. 

Rocket with centrifically releasable re- 
straining means holding a firing pin in 
cocked position to fire a percussion ele- 
ment when the rocket is spun prior to 


launching. A_ longitudinally movable 
cylinder is mounted to rotate on the spin 
axis in accordance with forces produced 
when the rocket is decelerated after 
launching. A friction member and friction 
surface generate heat to fire the explosive 
charge when the member is rotated around 
the surface. 


Control of air flow over an aircraft wing 
(2,868,479). M. Kadosch, J. LeFoll, 
F. M. L. Maunoury and J. Bertin, Nevilly- 
sur-Seine, France, assignors to SNECMA. 

Nozzles extending spanwise of a high 
camber wing and positioned at a distance 
above the upper surface, ahead of the trail- 
ing edge. Fluid pressure applied to the 
nozzles forms jets crosswise of the air flow, 
preventing separation of the flow from the 
wing surface. 

High lift supercirculation system using 
supersonic blowing (2,868,480). J. S. 
Attinello, Falls Church, Va. 

Method of increasing circulation while 
maintaining constant attack angle when 
the airfoil behaves in accordance with the 
Kutta-Joukowsky law. A stream of ex- 
panding and accelerating fluid is injected 
into the flow field above the airfoil at a 
velocity of at least Mach 1 to attract the 
flow in proximity to the injected stream. 
Ceramic bodies ¢ 2,868,658). W. W. Cof- 
fean, Westfield, N a assignor to Metal & 
Thermit Corp 

Body avian low electrical conductance 
consisting of 0.03 per cent (by weight) to 
0.3 per cent of at least one metal oxide 
(iron, nickel, zinc, cobalt, manganese or 
copper), 0.5 per cent to 50 per cent bis- 
muth oxide, and the remainder tin oxide. 
The analysis of the metal content is ex- 
pressed in corresponding metal oxides. 
Impact accelerometer (2,868,611). H- 
Carleton, Silver Spring, Md. 


Measuring device containing a magnet 
which locks a shaft and disk in place 
against the forces of a wound spring. A 
predetermined amount of acceleration 
causes the spring to rotate the disk while a 
stylus moves to mark a record in a soft 
metal indicating washer. 

Control system for dirigible craft (2,869,- 
063). J.J. Hess Jr., Garden City, N. Y., 
assignor to Sperry Rand Corp. 

Servomotor system for maintaining a 
craft at a predetermined reference attitude 
about an axis. A reversible servomotor 


modifies the normal reference provided by 


the system, providing a signal in accord 
ance with the output of the servomotor. 
bandpass filter receiver detects malfunc 
tion of the system, and means are provider 
for comparing the signals and controlling 
the craft. 
Hypergolic fuel for developing thrust 
(2,869,320). J. M. Burton, Waco, Texas 
assignor to Phillips Petroleum Co 
Bipropellant components injected sepa 
rately and simultaneously in contact witl 
each other in the combustion chamber of 
reaction motor to produce spontaneou 
combustion. Component consists of be 
tween 70 and 80 per cent (by volume) of 
pyrrole, and between 30 and 20 per cent of 
pyridine. 
Rocket fuse (2,869,465). J. G. Villepigue 
Inglewood, Calif., assignor to Northroy 
Aircraft, Inc. 


Impact switch connected with an energy 
source, an acceleration operable switch 
and thermally actuated means operative 
over a predetermined period after launch- 
ing a rocket. The impact switch is posi- 
tioned to close and detonate the charge 
after an additional time interval has 
elapsed 
Apparatus for nondestructive tests of 
initiators (2,869,364). I. Kabik and J. N. 
Ayres, Washington, D. C., assignors to the 
U.S. Navy. 

Potentiometer connected across a source 
of electrical energy, and a condenser in the 
circuit chargeable to a level below that 
required for firing. Means spatially dis- 
placed relative to the connected initiator 
for indicating the magnitude of the energy 
radiated by the conductive path. 

Vortex frequency airspeed indicator 
(2,869,366). A. Nitikman, Inglewood, 
Calif., assignor to Northrop Aircraft, Inc. 

Airfoil shaped member mounted near its 
leading edge to swing into the direction of 
the airstream. The streamlined leading 
edge separates the airstream and generates 
vortices. Sound produced by the vortices 
is transmitted to a microphone on the 
member. 
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Dyna-Soar (or dynamic soaring) is a joint 
project between the Air Force and the NASA, 
and is an attempt to solve the technical prob- 
lems of manned flight in the sub-orbital regions. 
Advance knowledge on the project indicates 
how a boost-glide vehicle can operate from the 
outer fringes of the atmosphere where it can 
maneuver and be recovered undamaged. Studies 
show that by varying the original rocket boost, 


Tmk(s) ® 


RADIO CORPORA TION of AMERICA 


DEFENSE ELECTRONIC PRODUCTS 
CAMDEN, N. J. 


and thus the velocity, and with the control 
available to the pilot, the Dyna-Soar aircraft 
can circumnavigate the earth, followed by a 
normal and controlled landing. Boeing Airplane 
Company, one of the competing companies for 
the development contract for the complete 
boost-glide system, has delegated to RCA the 
responsibility for the development of important 
electronic components of Dyna-Soar. 
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Force-ratio-responsive system (2,869,367 ). 
Moore, Pacific Palisades, Calif., 
assignor to Servomechanisms, Inc. 

with a pair of electrically 
spaced exciting windings and a rotor wind- 
ing. First and second signals 
_ vary with the forces, and are applied to 
the windings. A servo responsive to the 
signals produced adjusts the rotor to re- 

_ duce the signal to zero, the position of the 
_ rotor winding being representative of the 
ratio between two forces as indicated on a 


number scale. 

Optical (2,869,369). S. N. 
- Howell, Huntington, N. Y., assignor to 

Servo Corp. of America. 

Electrically responsive radiation py- 

--- rometer containing a heat sensitive ele- 
-- ment and collecting optics. The reading 
Of an indicator is optically displayed when 

looking into the sight. 
Balancing machine (2,869,372). 


K. E. 


Sihronen, Detroit, Mich., assignor to 
General Motors Corp. 

Means for measuring unbalance in a 

rotating workpiece. A pivoted reflector, 
directing a band of parallel light rays on a 
remote screen, moves in response to vibra- 
tions in the workpiece to produce a light 
pattern having a size proportional to the 
magnitude of unbalance. 
Means for detecting and suppressing ex- 
plosions (2,869,647). A. Mathisen, Lon- 
don, England, assignor of one-half to 
Graviner Mfg. Co. Ltd. and The Wilkinson 
Sword Co. Ltd. 

Photo-multiplier adapted to produce an 
electrical signal Siena exposed to the illumi- 
nation from an_ incipient explosion. 
Amplified current from the signal shatters 
a container of a liquid suppressant, pro- 
jecting the liquid at an initial velocity of 
150 fps. 


Vehicles (2,869,803). J. D. McGee, 


Ealing, London, England, assignor to 
Electrical & Musical Industries, Ltd. 

Remotely controlled space vehicle con- 
taining a television camera for transmit- 
ting signals indicating the position of a 
target relative to the course of the craft. 
Gyro controls maintain the craft on the 
course established from a remote control 
station. 


Apparatus for measuring accelerations 
and inclinations (2,869,851). H. B. 
Sedgfield, M. L. Jofeh and R. Albrecht, 
London, England, assignors to The 
Sperry Gyroscope Co. Ltd. 

Inertia element mounted in a casing for 
relative movement along the axis of an 
axial member. Signal generating means 
produce a signal dependent upon the dis- 
placement of the element from a_ null 
position along the axis exerting an electro- 
magnetic force restoring the member to its 
null position. 


Celestial Mechanics, edited by E. Finlay- 
Freundlich, Pergamon Press, New York, 
1959, 150 pp. $7.50. 
Reviewed by S. F. SINGER 
University of Maryland 


It is interesting to see the renewed 
appreciation of celestial mechanics which 
has been brought about by the rapid 
development of astronautics. That in 
itself would be sufficient reason to write 
new books on the subject. The present 
volume, however, was written for quite a 
different reason. The author, a dis- 
tinguished astrophysicist and astronomer, 
wants to stress the importance of celestial 
mechanics to astrophysical problems, such 
as the motion of close binary systems, 
and includes also revisions brought about 
by the theory of relativity. 

The result is a book which is very 
concise but, unfortunately, addressed 
primarily to the astronomer, or at least 
to someone who has had a considerable 
background in classical mechanics. The 
reader who is looking for guidance either 
for planning of space orbits or for reduc- 
tion of orbit observations will be dis- 
appointed. On the other hand, if he is 
looking for an elegant advanced presenta- 
tion of the principles of celestial mechanics, 


Ali Bulent Cambel, Northwestern University, Associate Editor 


social and political changes being wrought 
by the new source of energy, and to the 
general, nontechnical reader who desires 
to understand the implications in these 
revolutionary scientific developments now 
worldwide in their effects. 

The author has succeeded remarkably 
well in condensing within the short space 
of 200 pages a coverage—in many fields 
quite superficial—of the world energy re- 
sources of all types. The ever-increasing 
energy demands of industrialized countries 
are discussed in regard to future require- 
ments; the various potential sources 
beyond those of the heat forms of atomic 
energy, the various known methods of 
releasing atomic energy, are all discussed 
in a clear, completely nontechnical way. 

After treating many means of potential 
utilization of this form of energy, such as 
electric power generation, space heating, 
marine power, locomotives, airplane, road 
transport and medical applications, the 
author discusses the other side of the coin, 
namely, the risks involved in airborne 
contamination, isotope selective concen- 
trations, waste disposal problems and 
third party insurance liabilities inherent 
in such applications, and also the legal 
and sociological steps which must be 
taken if the maximum good is to be ex- 


the book can be heartily eres if pected from such applications. 


Finally, a short but illuminating dis- 


- eussion of the world impacts of the new 


The Atom and the Energy Revolution, hy 
Norman Lansdell, Philosophical Li- 
brary, New York, 1958, 200 pp. $6.00. 
Reviewed by Donatp H. LoucuripGE 
General Motors Corporation 


4 This work is addressed to the business- 
man who faces the necessity of accom- 
modating himself and his industry to the 
new factors introduced into commerce 
by the energy of the atom, to the citizen 
who must assist in making intelligent 
democratic decisions with respect to the 
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military, industrial and world trade effects 
being brought about by the atom and 
the alterations in state and private in- 
dustry’s relationship, provide impact 
which should prove most thought-provok- 
ing to the serious reader. 

If a negative criticism were to be men- 
tioned, it would be the discouraging effect 
on the average reader of so many tables 
and statistical data in the first few chap- 
ters. These should have been relegated 


to an appendix, with references thereto 


for smoother reading. 


An Introduction to Rocket Missile Pro- 
pulsion, by the Rocketdyne Training 
Department, published by Rocket- 
dyne, Canoga Park, Calif., 1958, 130 pp. 

Reviewed by T. SNYDER 
=n North Carolina State College 

This book is written for those who do 
not have a background in rocket propul- 
sion theory or in rocket engine operation. 
The book does a very fine job of present- 
ing the fundamental principles of rocket 
propulsion with sufficient clarity to be 
useful to a nontechnically trained reader. 
The scope of the book is illustrated by the 
chapter headings as follows: 1 and 2 
Rocket Fundamentals, 3 Rocket Missile 
Flight, 4 Liquid Propellants, 5 Liquid 
Propellant Feed Systems, 6 Thrust 
Chamber Assembly, 7 Rocket Engine 
Controls, 8 Solid Propellant Rockets. 

Each chapter includes a series of review 
questions at the end, and numerous ex- 
ample problems are included in the book. 
The terminology and vocabulary of mis- 
sile technology are introduced and defined 
in simple terms. Well-drawn_ illustra- 
tions add to the overall clarity of the 
presentation. The appendix contains a 
complete list of symbols and units used, a 
summary of the more important formulas 
used in the text, a condensed table of 
properties of the atmosphere, a rather 
detailed flow diagram of the V-2 engine 
including an engine operation sequence, 
and a flow diagram of the NAA 50-4.5 
rocket-sled engine. The book does not 
contain an extensive bibliography. 

The general level of the book, from a 
technical point of view, is quite low, but 
this statement is intended as one of evalu- 
ation rather than criticism. No mathe- 
matics more sophisticated than elemen- 
tary algebra is employed, and it is the 
broad, qualitative description of missile 
propulsion which gives the book signifi- 
cance. With one evening of rather rapid 
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As revolutionary in design 
and engineering as the technological transition 
from flying machine to space vehicle...Achieveda 
quarter of acentury ago, the onlyregulators which 
could be used for the first missile experiments 
in the Navy Experimental Station at Annapolis 
... Still today the standard for the critical control 
of high pressure fluids... Now on the threshold 
of space, Grove is uniquely prepared—through 
advance research and development—to meet 
the challenge of even higher pressures ahead. 


GROVE POWREACTOR DOME REGULATOR MODEL GH-408—5S0-6000 PSI INLE 3000 PSI REDUCED PRESSURE 


GROVE VALVE and REGULATOR COMPANY 


6529 Hollis St., Oakland 8, California » 2559 W. Olympic Bivd., Los Angeles 6, California 
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reading, a reader can cover the book and 
acquire a useful degree of understanding 
of rocket propulsion theory and the hard- 
ware components of rocket engines. 


The Gas Turbine Power Plant, The Turbo- 
prop, Turbojet, Ramjet, and Rocket 
Engines, vol. II of Aircraft and Missile 
Propulsion, by M. J. Zucrow, John Wiley 
and Sons, Inc., New York, 1958, 636 + 
xiv pp. $13.00. 

Reviewed by ALLEN E. Funs 
Northwestern University 


The title of volume II is a succinct book 
review in itself and is indicative of the 
thoroughness with which the topics have 
been covered. Gas turbines which are de- 
signed primarily to provide shaft horse- 
power are discussed in the first tenth of 
the book. This includes gas turbines for 
locomotives, electrical power generation 
and marine propulsion. It seems rather 
strange to find such a discussion in a book 
on aircraft and missile propulsion; how- 
ever, an understanding of these gas tur- 
bines is helpful for a mastery of the turbo- 
prop. All the cycles which can be ob- 
tained by the permutations and combina- 
tion of regenerative, nonregenerative, re- 
heat, nonreheat, intercooling, noncooling, 
open, semiclosed and closed cycles are 
analyzed for the ideal gas turbine. In 
the following chapter, the turboprop and 
stationary turbines are examined com- 
ponent-by-component, from the compres- 
sor, combustion chamber and turbine to 
the exhaust nozzle. Combining the per- 
formance of each component, the design 
performance of the complete gas turbine 
is discussed, along with the effects of 
changes in ambient conditions. The 
enthalpy of the products of combustion is 
tabulated to assist in numerical calcula- 
tions. The designer of a turboprop has 
the choice of how much of the available 


= power to proportion to the propeller and 


Technical Literature Digest__ 


ery M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
ae The James Forrestal Research Center, Princeton University 


Aerodynamics of Jet — 
Propelled Missiles 


Hypersonic Flow About a Thin Body 
of Revolution, by R. Timman, NATO, 
Advisory Group Aeron. Res. and Dev., 
Rep. no. 141, July 1957, 10 pp. 

Effects of Nose Angle and Mach Num- 
ber on Transition on Cones at Supersonic 
Speeds, by K. R. Czarnecki and Mary W. 
Jackson, NACA TN no. 4388, Sept. 1958, 
17 pp. 


Eprror’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 
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how much to the exhaust gases. The 
best way to make this choice is not dis- 
cussed in detail. 

Almost a half of the book is devoted 
to the turbojet. After presenting the 
thrust equation and the definitions of 
numerous turboject performance param- 
eters, static and in-flight performance 
at the design point is considered. In 
many books the discussion of similarity 
and its application to corrected mass 
flows, engine speeds, etc., is dismissed 
with a statement “It can be shown 
that ....’’ Not so in this book; there 
is an excellent section which deals with 
similarity starting with the fundamental 
concepts. Off-design performance, vari- 
able area nozzles and thrust augmenta- 
tion conclude the turbojet chapter. The 
dual-spool turbojet is not discussed; 
Prof. Zucrow probably saves that as a 
final exam problem for his students! 

A comparison of the turbojet and ram- 
jet leads off the brief chapter on ramjets. 
As in the previous chapters, only the 
component performance and not the de- 
sign details is presented. A considera- 
tion of component design is deferred to 
volume III. The chapter includes sec- 
tions on efficiencies, definition of per- 
formance parameters, the ideal ramjet, 
the deviations from ideal and a ramjet 
with variable geometry. 

From the front cover to the rear cover 
there is a spectrum of jet propulsion en- 
gines; the final chapter is concerned with 
the missile mainstay, the rocket motor. 
The aspects of rocket propulsion com- 
mon to both liquid and solid propellant 
rockets are considered first; this includes 
thrust equation, performance parameters 
and nozzle flow. Solid propellants and 
the internal ballistics of the solid propel- 
lant rocket are discussed. A familiarity 
with combustion and thermochemical 
equilibrium is essential background for 
understanding frozen or equilibrium ex- 


pansion. A section on these subjects is 
inserted between the solid and liquid 
propellant rocket discussions. Exterior 
ballistics of rocket propelled vehicles con- 
clude volume II. 

In summary, each chapter has near) 
identical format: An introductory sec- 
tion which gives the history and _per- 
spective, a section defining the pertinent 
parameters, and a discussion of com- 
ponent performance. With this back 
ground the ideal device is presented, fol 
lowed by the consideration of real effects 
Many figures and graphs enhance th 
value of the book, as do the numerou 
numerical examples solved in the text. 

Formulas for thrust, specific fuel con 
sumption, thrust coefficients are not de 
veloped as explicit functions of known 
quantities, i.e, Mach number, heat ad- 
dition and ambient conditions. Rathe: 
a stepwise calculation procedure is fol- 
lowed. For example, the nozzle exit 
velocity appears in the thrust equation: 
several supplementary calculations are re- 
quired to obtain this velocity. Of course, 
it can be argued that explicit functions 
become unwieldy and obscure rather 
than clarifying. In this book the conver- 
sion factors J and g are carried along 
through algebraic manipulations. By 
suitable choice of units, the need for / 
and g is eliminated easing the manipula- 
tion test—for example, using ft-lb in- 
stead of Btu for enthalpy dispenses with 
J. Ducted fans, turborockets and simi- 
lar jet propulsion engines might profitably 
have been included. 

The book is written for the senior or 
first-year graduate engineering student. 
The reader, whether he is a student or 
practicing engineer, will find here a book 
of continuing usefulness. As a textbook, 
it provides an understanding of the 
fundamentals and a store of information. 
As a reference, it provides numerous 
tables and graphs of important formulas. 


Solution of the Three Dimensional 
Laminar Boundary Layer Over the Surface 
of a Rotating Cone, by Ching-Sheng Wu, 
Princeton Univ., Dept. Aeron. Engng. Rep. 
no. 415 (AFOSR TN 58-301; ASTIA 
AD 154, 211), April 1958, 7 pp. 

Hypersonic Viscous Flow over Slender 
Cones, by Lawrence Talbot, Toyoki 
Koga and Pauline M. Sherman, NACA 
TN 4327, Sept. 1958, 34 pp. 

Inviscid Flow around a Blunt Body of a 
Reacting Mixture of Gases, Part A: 
General Analysis, by W. Lick, Rensselaer 
Polytech. Inst., Dept. Aeron. Engng., Tech. 
Rep. AE 5810 (AFOSR TN_ 58-522; 
ASTIA AD 158335), May 1958, 52 pp. 

Inviscid Hypersonic Flow Around 
Spheres and Circular Cylinders, by E. Dale 
Martin, Rensselaer Polytech. Inst., Dept. 


Tech. Rep. AE 5807 


Aeron. Engng., 
ASTIA AD 


(AFSOR TN _ 58-448; 
158254), April 1958, 70 pp. 

Preliminary Heat-transfer Studies on 
Two Bodies of Revolution at Angle of At- 
tack at a Mach Number of 3.12, by Nor- 
man Sands and John R. Jack, NACA TN 
4378, Sept. 1958, 29 pp. 

An Investigation of Supersonic Turbu- 
lent Boundary Layers on Slender Bodies 
of Revolution in Free Flight by Use of a 
Mach-Zehnder Interferometer and Shad- 
owgraphs, by Alvin Seiff and Barbara 
J. Short, NACA TN 4364, Sept. 1958, 57 
pp. 

Simultaneous Effects of Pressure Gra- 
dient and Transverse Curvature on the 
Boundary Layer along Slender Bodies of 
Revolution, by Michiru Yasuhara, Univ. 
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Development of Ultra-Refractories for Prolonged Service 
at 2000° F to 5000° F in Nuclear Flight Program 


Ceramic problems facing engineers and 
scientists at General Electric’s Aircraft 
Nuclear Propulsion Department equal — 
and transcend —the materials require- 
ments created by missile re-entry 
conditions. 

Consider this contrast: a nose cone 
plunging to earth at 15,000 mph reaches 
extreme temperatures for less than one 
minute... but ceramic components envis- 
aged for nuclear flight will encounter 
operating cycles in the 2000°F to 5000°F 
range for periods of hours or days. 

Current opportunities for ceramists in- 
clude both applied research and materials 
development: 


® LEAD TECHNICAL DIRECTION of major mate- 
rials development activity utilizing ceramics. PhD. 


& SUPERVISE APPLIED CERAMIC RESEARCH 
including fabrication processes, testing, develop- 
ment of novel ceramic materials. PhD. 


m@ CONDUCT RESEARCH in ceramic technology in 
particular fields of fabrication, chemical reactivity, 
high temperature testing, development of high 
strength ceramic materials. MS, BS. 


= TAKE PROGRAM PLANNING RESPONSIBILITY 
...including financial estimates... for major ce- 
ramic materials development activity. PhD, MS, BS. 


@ PROVIDE entire Applied Materials Research op- 
eration with scientific consulting advice in the 
fields of ceramics and physical chemistry related 
to aircraft nuclear propulsion. (Heavy experience 
plus graduate study essential.) 


CERAMISTS who value the opportunity to do original work 
with a company that fosters free inquiry and initiative, are 


invited to inquire about positions now open in the areas to the 


left. Please include salary requirements with resume. ‘vase 2 — <5. 


Write to Mr. P. W. Christos, Div. 35-MF : 


AIRCRAFT NUCLEAR PROPULSION DEPARTMENT B 


GENERAL@Q ELECTRIC 


P.O. Box 132 Cincinnati 15, Ohio 
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RESEARCH 
DEVELOPMENT 
PRODUCTION 


CESSNA 


MODEL 141 
MA-1A 


Continental Aviation & Engineering Corporation 
has more than 10 years’ experience with small 
gas turbine engines. ...Four versions of the J69 
= turbojet—delivering 920 to 1700 Ibs. thrust— 
are in operation today. . . . The Model 141 
turbine air compressor engine is widely used 
in jet aircraft ground support equipment. . 
Research and development programs are con- 
_ tinuously in process on uprated turbines, new 
turbine configurations, solid fuel propellents, 
and solid fuel ramjet 
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Tokyo, Aeron. Res. Inst., Rep. 335, Aug 
1958, 169 pp. 

The Effects of Blunt Leadi — on 
Delta Wings at Mach 5.8, by Kenneth F. 
Nicholson, Calif. Inst. Tech., Guggenhein 
Aeron. Lab., Hypersonic Res. Proj., Men 
45, June 1958, 67 pp. 

The Boundary-layer Equation for Axi- 
ally Symmetric Flow Past a Body of Re- 
volution—Motion of a sa by David 
Meksyn, J. Aero/Space vol. 25, Oct 
1958, pp. 631-634, 664. 


Heat Transfer and ca 
Fluid Flow 


Similar Solutions for the Compressib|- 
Boundary Layer on a Yawed Cylinder with 
Transpiration Cooling, by Ivan E. Bec! 
with, NACA T'N no. 4345, Sept. 1958, 7. 
pp- 

On the Mechanism of Thermal Ioni 
zation Behind Strong Shock Waves, b 
H. D. Weymann, Univ. Maryland, Ins 
Fluid Dynamics Appl. Math., TN 
BN-144 (AFOSR TN no. 58-788; AS 
TIA AD 2-3113,) July 1958, 36 pp. 

Small Perturbations in Hydromagnetics 
Free Oscillations in a Rectangular Box 
and the Perfectly Conducting Gas, b 
G. 8. 8. Ludford, Univ. Maryland, Insi 
Fluid Dynamics Appl. Math., TR no 
BN-143 (AFOSR TN no. 58-700; ASTI A 
AD 162, 234), July 1958, 24 pp. 

Theory and Experiments on Supersonic 
Air-to-Air Ejectors, by J. Fabri and J 
Paulon, NACA Tech. Mem. no. 1410 
Sept. 1958, 30 pp. (Translated fron 
France ONERA NT no. 36, 1956.) 

Determination of High Gas Tempera- 
tures at High Speeds, by E. Neher, Ar- 
chiv Tech. Messen no. 270, July 1958, pp. 
133-136. (In German.) 

Evaporation into a Boundary Layer, by 
Ernest Bauer and Martin Zlotnick, Phys. 
Fluids, vol. 1, July-Aug. 1958, pp. 355- 
356. 

Heat Transfer to a Sphere at the Transi- 
tion from Free Molecule Flow, by P. 
Hammerling and B. Kivel, Phys. Fluids, 
vol. 1, July-Aug. 1958, pp. 357-358. 

Higher Order Interactions in Homo- 
geneous Turbulence Theory, by Robert 
H. Kraichnan, Phys. Fluids, vol. 1, July- 
Aug. 1958, p. 358. 

Heat Transfer by Laminar Free Convec- 
tion in Enclosed Plane Gas Layers, by G. 
Poots, Quart. J. Mech. Appl. Math., vol. 
11, Aug. 1958, pp. 257-273. 

Three-dimensional Steady, Radial Flow 
of Viscous, Heat-conducting Compressible 
Fluid, by Akira Sakurai, Quart. J. Mech. 
Appl. Math., vol. 11, Aug. 1958, pp. 274 
289. 

Improved Method of Gamma-ray Calo- 
rimetry, by I. T. Myers, Rev. Sci. Inst. 
vol. 29, Sept. 1958, pp. 758-761. 

Calculation of the Laminar Boundary 
Layer on an Insulated Flat Plate by the 
Klunker-McLean Method, by Leslie M. 
Mack, Calif. Inst. Tech., Jet Propulsion 
Lab., Prog. Rep. no. 20-352, July 1952, 
120 pp. 

A Theoretical Study of Stagnation-point 
Ablation, by Leonard Roberts, NACA TN 
no. 4392, Sept. 1958, 29 pp. 

Aerodynamic Heating Characteristics 
of a Nose Cone at Angle of Attack, by 
William M. van Camp, McDonnell Air- 
craft Corp., Res. Dept., Rep. no. 5978, May 
1958, 175 pp. 

Interaction of a Shock Wave with a Tur- 
bulent Boundary Layer at M-3.85, by 
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July- depended on AMMONIUM PERCHLORATE from American 
Potash & Chemical Corporation since the very beginning. 
‘by G. First in the field with this essential oxidant, AP&CC was 
., vol. for many years the only domestic producer of ordnance- 
grade NH;ClO;. Today, with a growing network of 
l toe é strategically located plants and increased technical ; 
‘Mech ' knowledge, Trona still leads the industry. Supporting the © 
. 274 big tonnage production of AMMONIUM PERCHLORATE at 
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I. E. Vas and S. M. Bogdonoff, Princeton 
Univ., Dept. Aeron. Engng., Rep. no. 294 
(AFOSR TN 55-199), April 1955, 9 pp. 

The Performance of a Small Helical 
Tube Heat Exchanger, by J. C. Breeze, 
Atomic Energy Comm., CF-57-9-33, Dec. 
1957, 16 pp. 

Free Convection Heat Transfer to Hori- 
zontal Cylinders from an Ordinary Fluid 
Containing a Volume Heat Source, by 
Jose Pineda de Guzman, Atomic Energy 
Comm., AECU-3579, Jan. 1958, 95 pp. 
(M.S. Chem. Engng. ). 

Axisymmetric Stagnant Flow of a Vis- 
cous and Electrically Conducting Fluid 
Near the Blunt Nose of a Spinning Body 
with Presence of Magnetic Field Part I. 
Exact Solution of Incompressible and Con- 
stant Properties Model, by Ching-Sheng 
Wu and Wallace D. Hayes, Princeton 
Univ., Dept. Aeron. Engng., Rep. no. 431 
(AFOSR TN 58-405; ASTIA AD 
158208), April 1958, 44 pp. 

Steam Slip and Burnout in Bulk Boiling 
Systems, by A. E. Galson, Atomic Energy 
Comm., GEAP-1076, June 1957, 21 pp. 

Effect of Magnetic Field on Forced- 
convection Heat Transfer, by R. Siegel, 
J. Appl. Mech., vol. 25, Sept. 1958, pp. 
415-416. 

High-temperature Heat Conductivity of 
Some Metal Oxides, by John C. Jamieson 
and A. W. Lawson, J. Appl. Phys., vol. 28, 
Sept. 1958, pp. 1313-1314. 


On Reducing Aerodynamic Heat-trans- 
fer Rates by Magnetohydrodynamic Tech- 
niques, by Rudolph C. Meyer, J. Aero/- 
Sci., vol. 25, Sept. 1958, pp. 
572. 


The Effect of Extreme Cooling and 
Local Conditions on Boundary-layer 
Transition, by John R. Jack and Richard 
J. Wisniewski, J. Aero/Space Sci., vol. 25, 
Sept. 1958, p. 592. 

An Experimental Investigation of the 
Stability of the Hypersonic Laminar 
Boundary Layer, by A. Demetriades, J 
Aero/Space Sci., vol. 25, Sept. 1958, pp. 
599-600. 

Simple Test Method for Determining 
the Drop Size Distribution of a Hollow 
Cone Spray, by Hikmet Binark and W. E. 
Ranz, Project Squid TR no. PSU-2-P 
(ASTIA AD 200822), Aug. 1958, 13 pp. 
(Available only on microcards. ) 

Film Coefficients for Heat Transfer to 
Liquid Drops, by E. R. Elzinga Jr. and J. 
T. Banchero, AIChE Preprint no. 3, Aug. 
1958, 49 pp. 

Graphs of Reduced Variables for Com- 
puting Histories of Vaporizing Fuel Drops, 
and Drop Histories Under Pressure, rd 
G. L. Borman, M. M. El Wakil, O. A. 
Uyehara and P. S. Myers, NACA rN 
no. 4338, Sept. 1958, 55 pp. 

Void Volumes in Subcooled Boiling 
Systems, by P. Griffith, J. A. Clark and 
W. M. Rohsenow, AS ME Paper no. 58- 
HT-9, Aug. 1958, : 20 pp. 

Heat Transfer to a Boiling Liquid; 
Mechanism and Correlations, by K. E. 
Forster and R. Greif, ASME Paper no. 
58-HT-11, Aug. 1958, 33 pp. 

Bubble Growth Rates in Highly Sub- 
cooled Nucleate Boiling, by S. G. Bank- 
hoff and R. D. Mikesell, AJChE, Preprint 
no. 5, Aug. 1958, 12 pp. 

A Prediction of Surface Temperature at 
Incipient Boiling, by S. G. Bankoff, 
AIChE Preprint no. 4, Aug. 1958, 18 pp. 

Generalized Correlation of Boiling Heat 
Transfer, by S. Levy, ASME Paper no. 
58-HT-8, Aug. 1958, 6 pp. 

A Preliminary Study of Boiling Burnout 
Heat Fluxes for Water i in Vortex ee by 


W. R. Gambill and N. D. Greene, AIChE, 
Preprint no. 29, Aug. 1958, 29 pp. 

Selected Topics from the Theory of Gas 
Flow at High Temperatures. V. The 
Application of Transfer Equations to the 
Calculation of Diffusion, Heat Conduction, 
Viscosity and Electric Conductivity, by 
J. M. Burgers, Univ. Maryland, Inst. 
Fluid Dynamics Appl. Math., TN nos. 
BN-124a and b (AFOSR TN 58-427 and 
427a; ASTIA AD 158230 and 158230a), 
May 1958, 125 pp. 

The Second Approximation to the Stress 
Tensor and the Heat Flux Vector in a Gas 
for Krook’s Form of the Collision Equa- 
tion, by Hsun-Tiao Yang (supplement to 

“Selected Topics from the Theory of Gas 
Flow at Temperatures,’ Chapter IV), 
Univ. Maryland, Inst. Fluid Dynamics 
Appl. Math., TN no. BN-142 (AFOSR TN 
58-753; ASTIA AD 201612), July 1958, 
10 pp. 

Investigation of Boiling Burnout and 
Flow Stability for Water Flowing in Tubes, 
by Warren H, Lowdermilk, Chester D. 
Lanzo and Byron L. Siegel, NACA TN 
no. 4382, Sept. 1958, 51 pp. 

Analytical and Experimental Investi- 
gation of Temperature Recovery Factors 
for Fully Developed Flow of Air in a Tube, 
by R. G. Deissler, W. F. Weiland and 
W. H. Lowdermilk, NACA TN no. 4376, 
Sept. 1958, 35 pp. 

A Cooled Gas Pyrometer for Use in 
High-temperature Gas Streams, by Lloyd 
N. Krause, Robert C. Johnson and George 
E. Glawe, NACA TN no. 4383, Sept. 1958, 
32 pp. 

An Experimental Investigation of a 
Simplified Method for Measuring Heat 
Transfer at Supersonic Speeds, by Robert 
M. O’Donnell, Univ. Texas, Defense Res. 
Lab., CF 2709 (DRL-427), Aug. 1958, 20 
pp, 11 figs. 

Statistical Distribution of Solid Phase 
in Two Phase Turbulent Motion, by S. L. 
Soo and Richard L. Peskin, Proj. Squid, 
Tech. Rep. PR-80-R (ASTIA AD 158- 
714), May 1958, 55 pp. (Available only on 
microcards. ) 

Effects of Normal Surface Vibration on 
Laminar Forced Convective Heat Trans- 
fer, by J. A. Scanlan, /nd. Engng. Chem., 
vol. 50, Oct. 1958, pp. 1565-1568. 

Reduced Thermal Conductivity Cor- 
relation, Gaseous and Liquid Hydrogen, 
by C. A. Schaefer and George Thodos, Ind. 
Engng. Chem., vol. 50, Oct. 1958, pp. 
1585-15 

Electroconvectional Heat Transfer in 
Gases, by Sigurds Arajs and Sam Legvold, 
J. Chem. Phys., vol. 29, Sept. 1958, pp. 
531-536. 

Energy Transfer through a Dissociated 
Diatomic Gas in Couette Flow, by John 
F. Clarke, J. Fluid Mech., vol. 4, Sept. 
1958, pp. 441-451. 

On the Melting of a Semi-infinite Body 
of Ice Placed in a Hot Stream of Air, by 
Leonard Roberts, J. Fluid Mech., vol. 4, 
Sept. 1958, pp. 505-511. 

Nucleate Boiling—a Correlation, by 
Charles H. Gilmour, A/ChE Preprint 27 
Aug. 1958, 12 pp. 

Active Sites for Nucleate Boiling, by 
H. B. Clark, P. S. Strange and J. W. West- 
water, AIChE Preprint 13, Aug. 1958, 15 
pp. 

Forced Convection, Local Boiling Heat 
Transfer in Narrow Annuli, by Louis Ber- 
nath and William Begell, ATChE Preprint 
9, Aug. 1958, 13 pp. 

The Effects of Superimposed Forced 
and Free Convection in Horizontal and 
Vertical Rectangular Ducts, by M. Alt- 
man and F. W. Staub, AJChE Preprint 2, 
Aug. 1958, 39 pp. 


Condensing Heat Transfer within Hori- 
zontal Tubes, by W. W. Akers, H. A. 
Deans and QO. K. Crosser, AIChE Pre- 
print 1, Aug. 1958, 22 pp. ; 

Mass Transfer Cooling Near the Stag- 
nation Point, by Leonard Roberts, NACA 
TN 4391, Sept. 1958, 42 pp. 

A Literature Survey on the Applica- 
bility of Steady State Heat Transfer Co- 
efficients to Transients in Power or in Flow 
with Single Phase Flow and with Two 
Phase Flow, by R. H. Norris, AEC 
AECU-3611, Nov. 1957, 15 pp. 

Mechanically Aided Heat Transfer, b) 
D. Q. Kern and H. J. Karakas, AIChE 
Preprint 17, Aug. 1958, 24 pp. 

Internal Energy-Compressibility-Tem- 
perature Equation for Air in Dissociativ« 
and Ionic Equilibrium, by Perrin Win 
chell, J. Chem. Phys., vol. 29, Oct. 1958 
pp. 687-688. 

Statistical Mechanical Theory of Trans- 
port Processes in Liquids, by Frank C 
Collins and Helen Raffel, J. Chem. Phys. 
vol. 29, Oct. 1958, pp. 699-710. 


Effect of Combustion on Heat and Mas: 
Transfer in a Laminar Boundary Layer 
by E. R. G. Eckert and J. P. Hartnett, 
ZAMP., vol. 9B, March 25, 1958 (Jakob 
Ackeret anniv. vol. ), pp. 259-265. (In 
German. ) 

The Dissociation of a Pure Diatomic 
Gas —— a Strong Normal Shock Wave, 
by G. Jarre, ZAMP., vol. 9B, March 25, 
1958 (Jakob Ac keret anniv. vol.), pp. 

389-395. 

Magneto-gasdynamic Channel Flow, 
by E. L. Resler Jr. and W. R. Sears, 
ZAMP, vol. 9B, March 25, 1958 (Jakob 
Ackeret anniv. vol.), pp. 509-516. 

On Magneto-aerodynamic Boundary 
Layers, by V. J. Rossow, ZAMP, vol. 9B, 
March 25, 1958 (Jakob Ackeret anniv. 
vol.), pp. 519-527. 


Boundary Layer Stability Diagrams for 
Electrically Conducting Fluids in the Pres- 
ence of a Magnetic Field, by Vernon J. 
Rossow, NACA TN 4282, Aug. 1958, 32 
pp., diagrs., tabs. 

Interaction of a Turbulent Boundary 
Layer with a Step at M = 3.85, by I. E. 
Vas and S. M. Bogdonoff, Princeton Univ., 
Dept. Aeron. Engng., Rep. 295, April 
1955, 25 np. (Air Force, Office of Scien- 
tific Res. TN 55-200.) 

Influence of Solid-body Rotation on 
Screen-produced Turbulence, by Stephen 
3 Traugott, Johns Hopkins Univ., 
NACA TN 4135, Aug. 1958, 100 pp., 
diagrs., photos. 

Report on Boiling Heat Transfer (Rus- 
sian Work), by Novak Zuber, AEC, 
AECU-3569, Sept. 1957, 5 pp., 27 refs. 

A Study of Burnout Heat Fluxes As- 
sociated with Forced-convection, Sub- 
cooled and Bulk Nucleate Boiling of 
Water in Source-vortex Flow, by W. R. 
Gambill and N. D. Greene, AEC CF{57 
10-118, Oct. 29, 1957, 15 pp. 

Cryogenic Data Book, by Dudley B. 
Chelton and Douglas B. Mann, AEC 
UCRL-3421, May 15, 1956, 116 pp. 

On the Spectrum of Turbulence, by 
Donat G. Wentzel, Phys. Fluids, vol. 1, 
May-June 1958, pp. 213-214. 

Shock-front-thickness Measurements 
by an Electron Beam Technique, by H. N. 


Ballard and Douglas Venable, Phys. 
Fluids, vol. 1, May-June 1958, pp.’ 225- 
229. 


Relaxation Time for Reactions Behind 
Shock Waves and Shock Wave Profiles, by 
Russel E. Duff, Phys. Fluids, vol. 1, May— 
June 1958, pp. 242-244. 
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Variational Approach to Magnetohydro- 
dynamics, by Philip Rosen, Phys. Fluids, 
vol. 1, May—June 1958, p. 251. 

Magneto Hydrodynamic Distortion of a 
Magnetic Field Due to a Uniform Flow, 
by David S. Falk, AVCO Res. Lab., Res. 
Rep. 29, April 30, 1958, 5 pp. 

Ionization Phenomena of Shock Waves 
in Oxygen-nitrogen Mixtures, by S. 
Lin, AVCO Res. Lab., Res. Rep. 33, June 
1958, 23 pp., 15 refs. 

Measurements in a Shock Tube of 
Heat-transfer Rates at the Stagnation 
Point of a 1.0-inch-diameter Sphere for 
Temperatures to 7,900°F, 

Alexander Sabol, NACA TN 4354, 
n ug. 1958, 15 pp. 

Turbulence and Temperature Fluctua- 
tions Behind a Heated Grid, by R. R. 
Mills Jr., A. L. Kistler, V. O’Brien and S. 
Corrsin, NACA TN 4288, Aug. 1958, 67 
pp. 

Temperature and Thermal-stress Dis- 
tributions in Some Structural Elements 
Heated at a Constant Rate, by William 
A. Brooks Jr., NACA TN 4306, Aug. 1958, 
pp. 

Approximate Formulas for the gr 
and Thermal Conductivity of Gas Mix- 
tures, by Richard 8. Brokaw, J. Chem. 
Phys.. vol. 29, Aug. 1958, pp. 391-397. 

Turbulent Convection in a Pipe in the 
Case of a Fluid of High Prandtl Number, 
hy G. Ribaud, (Aerothermodynamics 
Seminar of the Paris Faculty of Science, 
1956-1957), France, Ministere de I Air, 
Publications Sci. et Tech., Note Tech., 73 
1958, 192 pp. (In French.) 

Effects of External Stream Flow and 
Afterbody Variations on the Performance 
of a Plug Nozzle at High Subsonic Speeds 
by R. J. Salmi and E. M. Cortright Jr., 
NACA Res. Mem. F56F lla, Oct. 1956, 
19 pp., diagrs. (Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 128, 7/28/58. ) 

Internal Characteristics and Perfor- 
mance of an Aerodynamically Controlled 
Variable-discharge Convergent Nozzle, 
by Jack G. MeArdle, NACA TN 4315, 
July 1958, 33 pp., diagrs., photo., tabs. 

An Example of Unsteady Laminar 
Boundary Layer Flow, by Itiro Tani, 
Univ. of Tokyo, Aeron. Res. Inst., Rep. 
331, (vol. 24, no. 2), June 1958, 42 pp. 

Extended Applications of Hot Wire 
Anemometry to High Speed Turbulent 
Boundary Layers, by Mark V. Morkovin 
and Ralph E. Phinney, Johns Hopkins 
Univ., Dept. Aeron., (. 1 FOSR-TN-58- 469; 
ASTIA AD 158279), June 1958. 

An Unsteady Turbulent Boundary 
Layer, ky Sture K. F. Karlsson, Johns 
Hopkins Univ., Dept. Aeron., 1958, 59 pp., 
20 figs. 

Boundaries of Supersonic Axisymmetric 
Free Jets, by Eugene 8S. Love, Mildred J. 
Woodling and Louise P. Lee, NACA Res. 
Mem. L56G18, Oct. 1956, 98 pp., diagrs. 
(Declassified from Confidential by author- 
ity of NACA Res. Abstr. 125, 2/26/58.) 

A Review of the Thermodynamic, 
Transport, and Chemical Reaction Rate 
Properties of High-temperature Air, by C. 
Frederick Hansen and Steve P. Heims, 
VACA TN 4359, July 1958, 33 pp., diagrs. 

Tables of Thermodynamic Properties of 
Argon Free Air to 15000°K, by Joseph 
Hilsenrath and Charles W. Beckett, 
Arnold Engng. Dev. Center, TN 56-12, 
(ASTIA AD 98974), Sept. 1956, 45 pp. 

Spreading of Exhaust Jet from 16-inch 
Ram Jet at Mach Number 2.0, by Fred 
Wilcox and Donald Pennington, NACA 
Res. Mem. E52F25, Aug. 1952, 14 pp., 
diagrs., photo., tab. (Declassified from 


ar \ Electronic systems engineering covers the full spectrum of 
- applied research and development at IBM. Currently, engi- 
_ neers and scientists concerned with broad systems problems 
in business, science, and government are working toward self- 
optimizing computers that may some day program themselves 
a and arrive at the one best solution to a problem. Progress is 
being made in advanced studies for radically different data 
systems for terrestrial and stellar navigational problems. For 
problems like these, and many others, IBM needs people who 
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You will enjoy unusual professional freedom, comprehensive 
education programs, the assistance of specialists of diverse 
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Confidential by authority of NACA Res. 
Abstr. 128, 7/28/58.) 

Cylindrical Shock Waves from Explod- 
ing Wires, by Frederick D. Bennett, 
Aberdeen Proving Ground, Ballistic Res. 
Labs., Rep. 1035, April. 1958, 19 pp. 

An Experimental Investigation of Vis- 
cous Effects on Static and Impact Pressure 
Probes in Hypersonic Flow, by Malcolm 
L. Matthews, Calif. Inst. Tech., Guggen- 
heim Aeron. Lab., Hypersonic Res. Proj., 
Mem. 44, June 1958, 38 pp. 

Free Field Measurements of Sound 
Radiated by Subsonic Air Jets, by Robert 
Lee, David Taylor Model Basin Rep. 868, 
Dec. 1953, 13 pp. 

Mass Transfer Cooling in a Laminar 
Boundary Layer in Steady Two Dimen- 
sional Stagnation Flow, by A. A. Hayday, 
Minn. Univ., Mech. Engng. Dept., Heat 
Transfer Lab., TN 19, (AFOSR-TN-58- 
337; ASTIA AD 154241), April 1958, 
17 pp., 19 figs. 

Bodies in an Expanded Electromagnetic 
High Frequency Field, by H. Stolz, An- 
nalen Physik, vol. 1, June 1958, pp. 334 
343. (In German.) 

Schumann-Runge Oxygen Absorption 
in Shock-heated Air, by W. H. Wurster, 
C. E. Treanor and H. 8. Glick, J. Chem. 
Phys., vol. 29, July 1958, pp. 250-256. 

Finite Amplitude Cellular Convection, 
by W. V. R. Malkus and G. Veronis, J. 
Fluid Mech., vol. 4, July 1958, pp. 225- 
232. 

Characteristics of the Equations of 
Motion of a Reacting Gas, by L. J. F. 
Broer, J. Fluid Mech., vol. 4, July 1958, 
pp. 276-283. 

The Thermal Boundary Layer on a Non- 
isothermal Surface with Non-uniform 
Free Stream Velocity, by E. M. Sparrow, 
J. Fluid Mech., vol. 4, July 1958, pp. 321- 
328. 

Detached Shock Waves Ahead of Gas- 
sampling Probes, by Robert Friedman 
and Donald R. Boldman, J. Aero/Space 
Sci., vol. 25, Aug. 1958, p. 526. 

On the Performance of a Double-dia- 
phragm Shock Tube Using the Reflected- 
shock Method and a Light-gas Buffer, by 
Charles J. Schexnayder Jr., J. Aero/- 
Space Sci., vol. 25, Aug. 1958, pp. 527- 
528. 

A New Function in the Theory of Fluids 
and an Equation of State for Liquids and 
Gases, by A. G. McLellan, Phil. Mag., 
vol. 3, July 1958, pp. 707-715. 

Classical Theory of Transport Phenom- 
ena in Dilute Poly-atomic Gases, by N. 
Taxman, Physical Rev., vol. 110, June 
15, 1958, pp. 1235-1238. 

Magnetic Field Effects on Bow Shock 
Stand-off Distance, by Richard W. Ziemer 
and William B. Bush, Physical Rev. 
Letters, vol. 1, no. 2, July 15, 1958, p. 58. 

Development of the Calorimeter Heat 
Transfer Gauge for Use in Shock Tubes, 
by Peter H. Rose, Rev. Sci. Instr., vol. 29, 
July 1958, pp. 557-560. 


Combustion, Fuels and 
Propellants 


Burning Velocities of Hydrocarbon 
Flames—Fuel Lean Mixtures, by S. A. 
Weil, J. R. Hasenberg and R. T. Ellington, 
Inst. of Gas Tech., June 1958, 83 pp. 

The Hyperscopicity of Lithium and Am- 
monium Nitrates and Perchlorates, by R. 
F. Muraca and L. L. Taylor, Calif. Fast 
Tech., Jet Prop. Lab., Progr. Rep. 20-347, 
Jan. 1958, 4 pp. 
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Effect of High Temperature Storage 
upon the Ballistic Stability of Small Arms 
Ammunition, by M. E. Levy, Frankford 
Arsenal, Pitman Dunn Labs. Group, Rep. 
R1433, Jan. 1958, 58 pp. 

Preparation and Physical Properties of 
Di- and Trinitrotoluene Isomers, by 
Cornelius Conklin and Frank Pristera, 
Picatinny Arsenal, Tech. Rep. 2525, July 
1958, 24 pp. 

The Underwater Shock-wave Initiation 
of Cast Pentolite, by C. H. Winning, Proc. 
Roy. Soc., London, vol. A 246, July 29, 
1958, pp. 288-298. 

Sensitivity of Explosives, VIII: The 
Effect of Compression of Occluded Gas 
Bubbles on the Initiation of Liquid Com- 
bustibles, by Julius W. Enig, Navord Rep. 
6093, June 2, 1958, 16 pp. 

Similitude in Aerothermochemistry, by 
Daniel E. Rosner, Aero. Chem. Res. Labs., 
Inc., Tech. Pub. 5, July 1958, 35 pp. 

The Analysis of Red Fuming Nitric 
Acid, I: A Critical Survey, by R. F. Mur- 
aca, Calif. Inst. Tech., Jet Prop. Lab., 
Progr. Rep. 20-343, Dec. 1957, 26 pp. 

Thermodynamic Properties of Some 
Boron Compounds, by William H. Evans, 
Donald D. Wagman and Edward J. Pro- 
sen, Nat. Bur. Standards, Rep. 4943, Aug. 
1956. 

The Vapor Pressures of Some Boron 
Compounds, by William H. Evans, Don- 
ald B. Wagman and Edward J. Prosen, 
Nat. Bur. Standards, Rep. 5663, Dec. 1957. 


Combustion of Elemental Boron, by 
Experiment, Inc., Quarterly Summary 
Rept. for Feb. through April 1958, 77M 
1038, May 1958, 7 pp. 

An Investigation of the Ignition Tem- 
peratures of Solid Metals, by W. C. Rey- 
nolds, Stanford Univ., Dept. Mech. Engngq., 
June 1957, 76 pp. 

Thermal Stability of Pentaborane in the 
Range 329° to 419°F, by Glen E. McDon- 
ald, NACA Res. Mem. £54G16, Dec. 1956, 
4 pp., diagrs. (Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 128, 7/28/58.) 

Cooperative Test Projects of the Joint 
Army Navy Air Force Panel on Analytical 
Chemistry of Solid Propellants, Part I: 
Introduction to the First Twelve Round 
Robins and Round Robin I, by R. H. 
Pierson, NOT'S 1937, Part I, Feb. 1958, 53 
pp. 

Dilution of Cryogenic Liquid Rocket 
Propellants During Pressurized Transfer, 
8. Greenfield, Aircr. Engng., vol. 30, July 
1958, pp. 210-215. 

Energy Requirements for the Ignition 
of Seven Solid Explosives, by G. J. Bryan 
and C. E. Noonan, Proc. Roy. Soc., Lon- 
don, vol. A 246, July 29, 1958, pp. 167-175. 


Initiation of Condensed Explosives by 
Compression of the Surrounding Gas, by 
J. I. Evans and A. M. Yuill, Proc. Roy. 
Soc., London, vol. A 246, July 29, 1958, 
pp. 176-179. 

Studies on the Oxidation Induced by 
Ultrasonic Radiation on a Solid System of 
Lauryl Aldehyde and Potassium Chlorate, 
by F. A. H. Rice and D. Levine, Proc. 
Roy. Soc., London, vol. A 246, July 29, 
1958, pp. 180-188. 

Nitromethane, Physical Properties, 
Thermodynamics, Kinetics of Decomposi- 
tion, and Utilization as Fuel, by A. Mak- 
ovky and L. Lenji, Chem. Rev., vol. 58, 
Aug. 1958, pp. 627-644. 

Determination of the 0—0 Bond Energy 
in Hydrogen Peroxide by Electron Impact, 
by L. P. Lindeman and J. C. Guffy, J. 
Chem. Phys., vol. 29, July 1958, pp. 247- 
250. 


Materials of Construction 


Proceedings. High Temperature Ma- 
terials, Their Strength Potentials and 
Limitations, Sagamore Ordnance Materials 
Research Conference, 4th, Aug. 21-23, 1957 
Co-sponsored by the Ordnance Materials 
Research Office and the Office of Ordnancc: 
Research, Arrangements by Syracuse Uni 
versity Research Institute (ASTJA Al 
157100), 348 pp. 

Elastomers for High Temperature Ap- 
plications, by E. R. Bartholomew, NATO 
AGARD, Rep. 178, March-April 1958, 2°; 
pp. 

Application of Structural Adhesives in 
Air Vehicles, by David L. Grimes, N ATO 
AGARD, Rep. 181, March-April 1958, 2s 
pp. 

A Review of the Development of Cer- 
mets, by G. C. Deutsch, A. J. Meyer Ji 
and G. M. Ault, NATO, AGARD, Re; 
185, March-April 1958, 20 pp. 

Fatigue and Ageing, by F. Bollenrat! 
NATO, AGARD, Rep. 157, Nov. 1957 
18 pp. 

Thermal Stresses in Thin Cylindrica 
Shells Stiffened by Plane Bulkheads fo 
Arbitrary Temperature Distributions, b: 
D. J. Johns, Cranfield, Coll. Aeron., Not 
83, July 1958, 22 pp. 

Thermal Insulation Ceramic Coatings 
by Alan V. Levy, Inst. Aeron. Sci., Pre 
print 857, Oct. 1958, 11 pp. 

Ceramic Materials—Properties fo 
Structural Applications, by W. J. Knap; 
and F. R. Shanley, Aero/Space Engng. 
vol. 17, Dee. 1958, pp. 34-38. 

Relief of Thermal Stresses Through 
Creep, by H. Poritsky and F. A. Fend, 
J. Appl. Mech. vol. 25, Dec. 1958, pp 
589-597. 


Instrumentation, Data 
Recording, Telemetering 


The influence of Modulation on the Re- 
cording of Resonance Phenomena, by F. 
Bruin and D. van Ladesteyn, Appl. 
Sci. Res., Section B., vol. 7, no. 4, 1958, 
pp. 270-274. 

Determination of the Dynamic Form of 
Pressure Indicators for Use in a Tube un- 
der High Pressure with Shock Waves, by 
J. Zucrow, Fusées et Récherche Aéron., vol. 
3, no. 2, July 1958, pp. 97-102. (In French. ) 

Current Distribution in Modulated Mag- 
netically Focussed Electron Beams, by 
M. Chodorow, H. J. Shaw and D. K. Wins- 
low, J. Appl. Phys., vol. 29, Nov. 1958, 
pp. 1525-1533. 

Dynamics of Electron Beams from Mag- 
netically Shielded Guns, by A. Ashkin, J 
Appl. Phys., vol. 29, Nov. 1958, pp. 
1594-1604. 

Measurement of Reflectance and Emis- 
sivity of Graphite at Arc Temperature with 
a Carbon Arc Image Furnace, by M. R. 
Null and W. W. Lozier, J. Appl. Phys., 
vol. 29, Nov. 1958, p. 1605. 

Properties of a Thermoelectric Cell, by 
G. M. Grover, D. J. Roehling, E. W. Salmi 
and R. W. Pidd, J. Appl. Phys., vol. 29, 
Nov. 1958, pp. 1611-1612. 

Disturbance Phenomena in Probe 
Measurement of Ionized Gases, by 
Takayoshi Okuda and Kenzo Yamamoto, 
J. Phys. Soc., Japan, vol. 13, Oct. 1958, 
pp. 1212-1223. 

Infrared Pyrometer for Temperature 
Measurements in Turbine Blades, by J. 
Crabol and J. Van Kote, La Récherche 
Aéron., no. 66, Sept.-Oct. 1958, pp. 3-12. 
(In French. ) 
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all, by Copper in production 
Beryllium in limited production 


e Reinforced plastics in development 


WYMAN-GORDON IS FORGING RE-ENTRY SHIELDS 


Probe 
by 
moto, 
is | UA SZ IVE A Gs ORD Oni 
rature FORGING 
by J. 
us ee ALUMINUM MAGNESIUM STEEL TITANIUM BERYLLIUM MOLYBDENUM COLUMBIUM 
3-12. AND OTHER UNCOMMON MATERIALS 
WORCESTER, MASSACHUSETTS 
GRAFTON, MASSACHUSETTS FRANKLIN PARK. ILLINOIS | 


DETROIT, MICHIGAN FORT WORTH, TEXAS ae ANGELES, CALIFORNIA 


foremost 

designers 
and 

manufacturers 


receivers : 


e CLARKE 


has been designing and manufacturing precision 


electronic equipment since 1909 when it first 
produced communications devices for the United 
States Government. This half-century of experi- 
ence is recognized by government and industry 
alike. Today more than 95% of the telemetry 


receivers in use at United States missile fest sta- 


tions and ranges were designed and built by 
Nems-Clarke. Among many installations now 


using this equipment are: 


WHITE SANDS MISSILE RANGE 
EGLIN AIR FORCE BASE 
ARMY BALLISTIC MISSILE AGENCY pane? 


— We welcome inquiries on 
problems in the telemetry field 


LAREKEH COMPANY 


DIVISION OF VITRO CORPORATION OF AMERICA 
SILVER SPRING. MARYLAND ° JUNIPER 5.1000 


| 
| 


Transistorized Relay Amplifier, by R. A 
Bruns and J. H. Wilcher, Calif. Inst. Tech., 
Jet Prop. Lab., Mem. 20-170, Sept. 1958, 
15 pp. 

A Note on the Magnetoresistance Ef- 
fect of Strain Gauge Wire, by Hideo 
Takaki and Toshiro Tsuji, J. Phys. Soc., 
Japan, vol. 13, Nov. 1958, p. 1406. 

Pressure Gauge for Corrosive Gases in 
the Micron and Submicron Region, by J 
R. Anderson, Rev. Sci. Instr., vol. 29, Dec. 
1958, pp. 1073-1078. 

Radiation Attenuation Method of Meas- 
uring Density of a Two-phase Fluid, b. 
Michael Petrick and Bernet S. Swanson, 
Rev. Sci. Instr., vol. 29, Dec. 1958, pp 
1079-1085. 

Small Oven for High-Temperature Op- 
tical Measurements, by D. A. Patterso:, 
and H. 8. Goulart, Rev. Sci. Instr., vol. 29 
Dee. 1958, pp. 1141-1142. 

On-off Temperature Control for Electri 
cally Heated Filaments, by Roger E. Litt] 
and J. D. McKinley Jr., Rev. Sct. Instr 
vol. 29, Dee. 1958, pp. 1143-1144. 

Improved Four-Pi Porportional Ga: 
Flow Counter, by Troy C. Martin and O 
E. Green, Rev. Sci. Insir., vol. 29, Dec 
1958, p. 1147. 

Simulation Techniques in Aeronautics 
by J. J. Foody and R. J. A. Paul., J. Roy 
Aeron. Soc., vol. 62, Dec. 1958, pp. 878 
889. 


Guidance Systems and 


Components 


Instrumentation Aspects of Inertial 
Guidance, by C. S. Draper, JSA J. 
(Instrum. Soc. Amer.) vol. 5, Nov. 1958, 
pp. 54-61. 

Radio-Astronomy and Navigation, by C. 
M. Cade, J. Roy. Aeron. Soc., vol. 62, Nov. 
1958, pp. 805-828. 

On the Representation of the Trajectory 
Band and Its Perturbation for Fire-Direct- 
ing Apparatus, by E. Roth-Desmeules, 
Zeitscrift fiir Angewandte Mathematik und 
Physik, vol. 9a, no. 3, Sept. 25, 1958, pp. 
235-250. (In German. ) 

Feasibility and Limits of Infrared View- 
ing, by T. Schneider, Zeitschrift fiir Ange- 
wandte Mathematik und Physik, vol. 9a, 
no. 3, Sept. 25, 1958, pp. 251-259. 

Inertial Navigation, by Ted Newman, 
Arma Engng., vol. 2, Oct.-Nov. 1958, pp. 
8-14. 

Aided Inertial Systems, by C. J. Mundo, 
Arma. Engng., vol. 2, Oct.-Nov. 1958, pp. 
15-19. 

Formula for Inertial Guidance Gyro 
Evaluation, ‘100 Tons of Concrete Plus 2 
Seconds of Arc,’’ by Aaron Becker, Arma 
Engng., vol. 2, Oct.-Nov. 1958, pp. 20-24. 

Real-time Airborne Digital Computers, 
by Gerhard L. Hollander, SA J. (Instrum. 
Soc. Amer.), vol. 5, Dec. 1958, pp. 26-31. 

Information Theory Plays Major Role in 
New Radar Designs, by James P. Hola- 
han, Space/Aeron., vol. 30, Dec. 1958, pp. 
138-141, 144-149. 

Lear Designs Guidance for Manned 
Satellite, Missiles and Rockets., vol. 5, Jan. 
5, 1959, p. 26. 

State of the Art: Inertial Guidance, by 
W. G. Wing, Space/Aeron., vol. 31, Jan 
1959, pp. 24, 82-87. 

State of the Art: Radar, by James 
Holahan, Space/Aeron., vol. 31, Jan. 
1959, pp. 25, 94-96, 98-104. 

State of the Art: Infrared, by Raymond 
H. McFee, Space/Aeron., vol. 31, Jan. 
1959, pp. 26, 88-93. 
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Spiral Antennas Meet Needs of New 
Vehicles, by Edwin M. Turner, Space/- 
Aeron., vol. 31, Jan. 1959, pp. 56-62. 

Far Field Scattering from Bodies of 
Revolution, by K. M. Siegel, Appl. Sci. 
Res., Sect. B., vol. 7, no. 4, 1958, pp. 293- 
328. 


Flight Vehicle Design 
and Testing 


Booster Case Design for Hypersonic 
Vehicles, by Peter E. Grafton and Eugene 
I. Styer, Inst. Aeron. Sci., Preprint 835, 
July 1958, 24 pp. 

Optimum Thrust Direction for Maxi- 
nium Range, by G. Leitmann, J. Brit. 
I nterplanet. Soc., vol. 16, Sept.-Oct. 1958, 
pp. 503-507. 

Missiles Emerging and Missles Re-en- 
tering, by P. Magot-Cuvru, Fusées et 
Récherche Aéron., vol. 3, no. 2, July 1958, 
pp. 53-58. 

Studies and Tests of a Liquid Propellant 
Ground-to-ground Missile and a Quasi- 
rectilinear Trajectory, by C. Belmundo, 
Pusées et Récherche Aéron., vol. 3, no. 3, 
July, 1958, pp. 65-68. 


e 
Servomechanisms and | 
Controls 


D-C Amplifiers for Control Systems, by 
Larry 8. Klivans, Electronics, vol. 31, Nov. 
21, 1958, pp. 96-98. 

Inertial Navigation, Part II, by C. F. 
0’ Donnell, J. Franklin Inst., vol. 266, Nov. 
1958, pp. 373-402. 

Servo Phase Control Shapes Antenna 
Pattern, by E. W. Markow, Electronics, 
vol. 32, Jan. 2, 1959, pp. 50-52. 

Pneumatic Actuating System for 1000° 
Operation, by John A. Osterman, /nst. 
Aeron. Sci., Preprint 846, Oct. 1958, 14 
pp. 

Short Account of the Theories for the 
Study of Servo Systems Defined by Linear 
Differential Equations with Constant Co- 
efficients and Influenced by Stochastic 
Disturbances, by M. Sundstrom and K. 
Flodin, Sweden, Roy. Inst. Techn., Dept. 
Aeron. Engng., TN 46, 1957, 11 pp. 

Relations between Spectral Densities of 
Inputs and Outputs in Linear Servo Sys- 
tems with Time Variable Coefficients, by 
M. Sundstrom, Sweden, Roy. Inst. Techn., 
Dept. Aeron. Engng., TN 45, 1957, 7. pp. 

Hypersonic Digital Control Systems, by 
A. E. Cooper, Inst. Aeron. Sci., Preprint 
839, July 1958, pp. 32. 

On the Use of Spectral Densities in the 
Study of Servo Circuits Defined by Linear 
Differential Equations with Time Variable 
Coefficients and Stochastic Right Member, 
by M. Sundstrom and K. Flodin, Sweden, 
Roy. Inst. Tech., Dept. Aeron. Engng., TN 


44, 1957, 9 pp. sles 


Space Law and Metalaw: Jurisdiction 
Defined, by Andrew G. Haley, J. Brit. 
Interplanet. Soc., vol. 16, no. 8, July-Aug. 
1958, pp. 472-479. 

The Rotation of the First Russian 
Earth-satellite, by J. H. Thomson, Phil. 
Mag., vol. 3, no. 32, Aug. 1958, pp. 912-916. 


JUNE 1959 


Hearings on H.R. 11881, 85th Congress, 
2nd Session, Select Committee on Astro- 
nautics and Space Exploration, April 
15—May 12, 1958, 1542 pp. 

Influence of the Self-rotation of the 
Central Body on the Motion of Satellites 
after the Einstein Gravitation Theory, 
by N. St. Kalitzin, Jl Nuovo Cimento, 
vol. 9, no. 3, Aug. 1, 1958, pp. 365-374. 
(In German.) 

Fundamentals of Inertial Guidance and 
Navigation, by William E. Frye, J. 
Astronaut. Sci., vol. 5, Spring 1958, pp. 
1-4, 10. 

Plastic Balloons for Planetary Research, 
by Malcolm D. Ross, J. Astronaut. Sci., 
vol. 5, Spring 1958, pp. 5-10. 

Space Cabin Design, by Alfred M 
Mayo, J. Astronaut. Sci., vol. 5, Spring 
1958, pp. 11-14. 

A Survey of Propulsion and Space 
Dynamics, by Douglas B. Cross and 
Jorgen Jensen, J. Astronaut. Sci., vol. 5, 
Spring 1958, pp. 15-19, 30. 

Basic Sciences Must be Part of Space 
Team, by L. Eugen Root, Missiles and 
Rockets, vol. 4, Oct. 6, 1958, pp. 31-33. 

Departure and Return in Interplane- 
tary Flight, by Karel J. Bossart, Aero/- 
Space Engng., vol. 17, Oct. 1958, pp. 44- 
52. 


The Effect of the Earth’s Oblateness on 


_ the Orbit of a Near Satellite, by D. G. 
_ King-Hele, Proc. Roy. Soc., vol. 247A, no. 


1248, Sept. 9, 1958, pp. 49-72. xy- 


Atomic Energy 


Fusion Chain Reaction—Chain Reac- 
tion with Charge Particles, by Michal 
Gryzinski, Physical Rev., vol. 111, no. 3, 
Aug. 1, 1958, pp. 900-906. 

The Containment of Plasma by the 
Pinch Discharge, by Sir George Thomson, 
Phil. Mag., vol. 3, no. 32, Aug. 1958, pp. 
886-896. 

The Stellarator Concept, by 
Spitzer Jr., Phys. Fluids, vol. 1, 
July-Aug. 1958, pp. 253-264. 

Equilibrium of a Magnetically Confined 
Plasma in a Toroid, by M. D. Kruskal and 
R. M. Kulsrud, Phys. Fluids, vol. 1, no. 
4, July-Aug. 1958, pp. 265-274. 

On the Stability of Plasma in Static 
Equilibrium, by M. D. Kruskal and C. R. 
Oberman, Phys. Fluids, vol. 1, no. 4, 
July-Aug. 1958, pp. 275-280. 

Some Stable Hydromagnetic Equilib- 
ria, by J. L. Johnson, C. R. Oberman 
R. M. Kulsrud and E. A. Frieman, Phys. 
Fluids, vol. 1, no. 4, July-Aug. 1958, pp. 
281-296. 

On the Ionization and Ohmic Heating 
of a Helium Plasma, by J. M. Berger, I. 
B. Bernstein, E. A. Frieman and R. M. 
Kulsrud, Phys. Fluids, vol. 1, no. 4, July- 
Aug. 1958, pp. 297-300. 

Heating of a Confined Plasma by Oscil- 
lating Electromagnetic Fields, by J. M. 
Berger, W. A. Newcomb, J. M. Dawson, 


Lyman 
no. 4, 


3 A. Frieman, R. M. Kulsrud and A. 


Lenard, Phys. Fluids, vol. 1, no. 4, July- 


Aug. 1958, pp. 301-307. 


Generation and Thermalization of Plas- 
ma Waves, by Thomas H. Stix, Phys. 
Fluids, vol. 1, no. 4, July-Aug. 1958, pp. 
308-317. 

On the Theory of Stochastic Processes 
in Nuclear Reactors, by L. Pal, Nuovo 
Cimento, Supplemento, no. 1, 1958, pp. 
25-42. (In English.) 

Plasma Loss from Magnetic Bottles, 
by E. Persico and J. G. Linhart, J! Nuovo 
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Help us solve } 

the problems of a, 

PROPULSION 


Our Advanced Design Group is 
seeking engineers and scientists to 
work on low-thrust propulsion sys- 
tems that will make possible the 
manned exploration of Outer Space. 
These are a few of the challenging 
positions : 


Theoretical Physicist 


Experienced in the field of particle 
physics, magneto-hydrodynamics and : 
plasma physics, comparative perform- + 
ance of electrical propulsion systems, 
new methods of developing thrust 
from electrical propulsion systems, 
and direct conversion processes. 


Chemical Physicists 
To make studies on propellants for 
ion propulsion. Experienced in elec- 
tronics, atomic physics, and physical 
chemistry. 
Physicists 
Primary activities will be in the 
field of ion sources and high vacuum 
techniques. Experienced in particle 
accelerators, ion sources, and elec- 
tronics. 
Electrical Engineer 
Analysis of electrical power gen- 
erators, preliminary stages designs of 
permanent magnet, AC induction and 
electrostatic generators. 
Please write: Mr. AF. Jamieson, 
Engineering Personnel Dept., 6633 
Canoga Avenue, Canoga Park. Calif. 


ROCKETDYNE 


A DIVISION OF NORTH AMERICAN AVIATION, INC. 
FIRST WITH POWER FOR OUTER SPACE 
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ENGINEERS 
SCIENTISTS 


Here is your opportunity to grow 
with a young, expanding subsidiary 
of the Ford Motor Company. Out- 
standing career opportunities are 
open in Aeronutronic’s new RE- 
SEARCH CENTER, overlooking the 
Pacific at Newport Beach, and the 
facility in Glendale, California. 
You will have all the advantages of 
a stimulating mental environment, 
working with advanced equipment 
in a new facility, located where you 
can enjoy California living at its 
finest. 


PhD and MS RESEARCH SPECIAL- 
ISTS with 5 to 7 years’ experience in heat 
transfer, fluid mechanics, thermodynamics, 
combustion and chemical kinetics, and 
thermoelasticity. , To work on theoretical 

and related to re- 
technology advanced rocket pro- 
pulsion. Specific assignments are open 
in re-entry body design, high temperature 
materials studies, boundary layer heat 
transfer with chemical reaction, thermal 
stress analysis, and high temperature 
thermodynamics. 


PROPULSION ENGINEERS with 5 years’ 
experience in liquid and solid rocket design 
and test. Familiarity with heat transfer 
problems in engines desirable. To work 
on program of wide scope in R & D of ad- 
vanced concepts in rocket engine com- 
ponents, and for missile project work. 


APPLIED MATHEMATICIANS, 3-5 
years’ recent experience required and A.B. 
or M.A. degree. Experience in numerical 
analysis and ter work in 

with rockets and rocket propellants. 


STRESS HIGH TEMPERATURE MA- 
TERIALS. Mechanical or Metallurgical 
Engineer, must be familiar with thermal 
stress and shock as well as elasticity and 
plasticity. Application to nose cones and 
thrust chambers. 


FLIGHT TEST & INSTRUMENTATION 
ENGINEERS with 5 to 10 years’ experi- 
ence in laboratory and flight test instru- 
ill develop tech- 
niques utilizing advanced instrumentation 
associated with space vehicles. 


THEORETICAL AEROTHERMODY- 
NAMICIST. Advanced degree and at 
least 5 years’ experience in high-speed 
aerodynamics. Knowledge of viscid and 
gas flows required. To work on program 
leading to advanced missile configurations. 
Work involves analysis of the re-entry of 
hypersonic missiles and space craft for 
determining optimum configuration. 


CERAMICIST. M.S. or Ph.D. required 
and 3-5 years’ recent experience with high 
temperature materials, structures and cer- 
mets. 


ENGINEER or PHYSICIST. With ex- 
perience in the use of scientific instruments 
for making physical measurement. Work 
related to flight test and facility instrumenta- 
tion. Advanced degree desired with mini- 
mum of 3 years of related experience. 


Qualified a; 
résumés and inquiries to Mr. R. W. Speich, 
Aeronutronic Systems, Inc. 


AERONUTRONIC 


a subsidiary of Ford Motor Company 


icants are invited to send 


2234 Air Way Bidg. 19, Glendale, Calif. 
CHapman 5-6651 


Newport Beach, Glendale, Sante Ana, 
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Cimento, vol. 8, no. 5, June 1, 1958, pp. 
740-753. (In English. ) 

Dynamic Stability of a Self-pinched 
Discharge, by H. W. Wyld Jr., J. Appl. 
Phys., vol. 29, Oct. 1958, pp. 1460-1465. 

Research on Power from Fusion and 
Other Major Activities in the Atomic 
Energy Programs, Atomic Energy Com- 
mission, Jan.-June 1958, 410 pp. 

Hydromagnetic Instability in a Stellar- 
ator, by M. D. Kruskal, J. L. Johnson, 
M. B. Gottlieb and L. M. Goldman, Phys. 
Fluids, vol. 1, no. 5, Sept.-Oct. 1958, pp. 
421-429. 

A Programme for Solving the Multi- 
group Neutron Diffusion Equations in 
Two Space Dimensions on the Ferranti 
Mercury Computer, by A. Hassitt, Gt. 
Brit., Atomic Energy Res. Estab., AERE 
T/R 2487, Feb. 1958, 39 pp. 

A High Current Proton Source, by 
P. C. Thonemann and E. R. Harrison, 
Gt. Brit., Atomic Energy Res. Estab., 
AERE GR/R 1190, Jan. 1955 (published 
1958) 15 pp. 

Special Fuel Distribution and Cooling 
Problems in the Reactors of Nuclear 
Rockets, by F. Winterberg, Astronautica 
Acta, vol. 4, no. 2, 1958, pp. 138-165. 
(In German. ) 

Hydromagnetic Shocks Used in Nuclear 
Fusion Engine, Aviation Age, vol. 30, 
Aug. 1958, pp. 30-34. 

Cleanup of Atomic Hydrogen, by Kath- 
erine B. Blodgett, J. Chem. Phys., vol. 29, 
July 1958, pp. 39-44. 

Energy Spectrum of Electrons Emitted 
from Gases Bombarded by Positive Ions, 
by David E. Moe and Otto H. Petsch, 
Phys. Rev., vol. 110, June 15, 1958, pp. 
1358-1362. 

Neutron Production in Linear Deute- 
rium Pinches, by Oscar A. Anderson, Wil- 
liam R. Baker, Stirling A. Colgate, "John 
Ise Jr. and Robert V. Pyle, Phys. Rev., 
vol. 110, June 15, 1958, pp. 1375-1378. 

Instability, Turbulence and Conduc- 
tivity in Current-carrying Plasma, by O. 
Buneman, Phys. Rev. Letters, vol. 1, no. 1, 
July 1, 1958, p. 8. 

Neutrons of Possible Thermonuclear 
_— by W. C. Elmore, E. M. Little and 

E. Quinn, Phys. Rev. Lelters, vol. 1, 
July 1, 1958, pp. 32. 

R-Matrix Theary of Nuclear Reactions, 
by A. M. Lane and R. G. Thomas, Rev. 
Modern Phys., vol. 30, no. 2, Part 1, April 
1958, pp. 257-277. 

Summary of High-energy Nucleon- 
nucleon Cross-section Data, by Wilmot 
N. Hess, Rev. Modern Phys., vol. 30, no. 2, 
Part 1, April 1958, pp. 368-471. 

Electromagnetic Structure of the Neu- 
tron, by L. I. Schiff, Rev. Modern Phys., 
vol. 30, no. 2, Part 1, April 1958, pp. 
462-472. 

Neutron-electron Interaction, by Leslie 
L. Folday, Rev. Modern Phys., vol. 30, no. 
2, Part 1, April 1958, pp. 471-481. 

Electromagnetic Structure of the Proton 
and Neutron, by R. Hofstadter, F. Bu- 
miller and M. R. Yearian, Rev. Modern 
Phys., vol. 30, no. 2, Part 1, April 1958, 
pp. 482-489. 

Direct Cycle Nuclear Turbojet Power 
Plants, by Eleanor L. Semple and William 
C. Cooley, Aero/Space Engng., vol. 17, 
Aug. 1958, pp. 30-36. 


Aerophysics, Astrophysics 


Hemispherical Distribution of Cosmic 
Rays at 25° Geomagnetic Latitude, by 
P. S. Gill and A. N. Mitra, Il Nuovo 


Cimento, vol. 9. no. 3, Aug. 1, 1958, pp. 
400-411. 

Cosmic Rays in the Earth’s Magnetic 
Field, by P. Rothwell, Phil. Mag., vol. 3, 
Sept. 1958, pp. 961-970. 

The Reflexion of Radio Waves from a 
Stratified Ionosphere Modified by Weak 
Irregularities, by M. L. V. Pitteway, 
Proc. Roy. Soc., vol. 246A, no. 1247, 
Aug. 26, 1958, pp. 556-562. 

Sputnik as a Tool for Securing Geodetic 
Information, by Louis Gold, J. Franklin 
Inst., vol. 266, no. 2, Aug. 1958, pp. 103- 
107. 

Density and Velocity Distribution of the 
Interstellar Gas, by C. van de Hulst 
Rev. Modern Phys., vol 30., no. 3, July 
1958, pp. 913-930. 

Twenty-one-cm Studies of Some Inter- 
stellar Clouds, by R. D. Davies, Rev 
Modern Phys., vol. 30, no. 3, July 1958, 
pp. 931-939. 

Validity of the Cloud Model of the 
Interstellar Medium, by Bertram Donn, 
Rev. Modern Phys., vol. 30, no. 3, July 
1958, pp. 940-942. 

On the Formation of Interstellar Gas 
Clouds, by S. A. Kaplan, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 943- 
951. 

Microstructure of the Galactic Magnetic 
Field, by K. Serkowski, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 952- 
954. 

Gross Dynamics of the Interstellar 
Medium, by E. N. Parker, Rev. Modern 
Phys., vol 30, no. 3, July 1958, pp. 955- 
965. 

Aerodynamic Dissipation, by H. E. 
Petschek, Rev. Modern Phys., vol 30, no. 
3, July 1958, pp. 966-974. 

Magnetohydrodynamic Dissipation, by 
L. Biermann and A. Schluter, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 975- 
978. 

Thermal Inelastic Collision Processes, 
by M. J. Seaton, Rev. Modern Phys , vol. 
30, no. 3, July 1958, pp. 979-1008. 

On the Cooling of Interstellar Matter, 
by E. Schatzman, Rev. Modern Phys., 
vol. 30, no. 3, July 1958, pp. 1009-1011. 

On the Condensation of Interstellar 
Material, by E. Schatzman, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1012- 
1014. 

On the Gravitational Instability of a 
Medium in Nonuniform Rotation, by N. 
Bel and E. Schatzman, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1015- 
1024. 


Internal Kinematics of the Planetary 
Nebulae, by O. C. Wilson, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1025- 
1032. 

Comparison of Gas and Radiation 
Pressure in the Problems of the Dynamics 
of Planetary Nebulae, by H. Zanstra, Rev. 
Modern Phys., vol. 30, no. 3, July 1958, pp. 
1030-1033. 


Local Density Variations in Planetary 
Nebulae, by M. J. Seaton, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, p. 
1034. 

Internal Motions in the Orion Nebula, 
by G. Munch, Rev. Modern Phys., 
vol. 30, no. 3, July 1958, pp. 1035-1041. 


Kinematics of the Filaments in the 
Crab Nebula, by G. Munch, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1042- 
1046. 

On the Nature of the Emission of the 
Crab Nebula, by I. 8S. Shklovsky, Rev. 
Modern Phys., vol 30, no. 3, July 1958, p. 
1047. 
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Cygnus Loop and Some Related Nebu- 
losities, by R. Minkowski, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1048- 
1052. 

Bright Rims in Diffuse Nebulae, by 
Stuart R. Pottasch, Rev. Modern Phys., 
vol. 30, no. 3, July 1958, pp. 1053-1057. 


On the Stability of Ionization Fronts, by 
}’. D. Kahn, Rev. Modern Phys., vol. 30, 
no. 3, July 1958, pp. 1058-1061. 


On the Propagation and Structure of 
Ionization Fronts, by F. 
Rev. Modern Phys., vol. no. 3, July 
1958, pp. 1062- 1068. 


Collision of Two Highly Ionized Clouds 
cf Gas, by F. D. Kahn, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1069- 
072. 


Gas-dynamical Effects of Star Forma- 
‘ion, by M. P. Savedoff, Rev. Modern 
~*hys., vol. 30, no. 3, July 1958, p. 1073. 

Gas Dynamics of Galaxy Collisions, 
by M. P. Savedoff, Rev. Modern Phys., 
vol. 30, no. 3, July 1958, pp. 1074-1076. 


Examples of Gas Motion and Certain 
Hypotheses on the Mechanism of Stellar 
Outbursts, by L. I. Sedov, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1077- 
1079. 

Some Exact Solutions of the Equations 
of Magnetohydrodynamics when Both 
Self-attraction and Magnetic Fields Are 
Present, by G. C. MeVittie, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1080- 
1083. 


On the Mathematical Structure of 
Turbulence, by J. Bass, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1084— 
1086. 

Distortion of a Toroidal Field by Con- 
vection, by D. W. Allan and E. C. Bullard, 
Rev. Modern. Phys., vol. 30, no. 3, July 
1958, pp. 1087-1088. 

Shock Waves in Magnetogasodynamic 
Turbulence, by 8S. A. Kaplan, Rev. 
Modern Phys., vol. 30, no. 3, July 1958, p. 
1089. 


Possible Hydromagnetic Simulation of 
Cosmical Phenomena in the Laboratory, 
by Winston H. Bostick, Rev. Modern 
Phys., vol. 30, no. 3, July 1958, pp. 1090- 
1094. 


Russian Technical Articles’ 


On the Stability of Permanent Rotations 
of a Body Around a Fixed Point, by V. V. 
Rumyantsev, Prikladnaya Matematika 7 
Mekhanika, vol. 21, no. 3, 1957, pp. 339- 
346. (In Russian.) 

On Bodies of Minimal Resistance in a 
Supersonic Gas Flow, by M. N. Kogan, 
Prikladnaya Matematika i Mekhanika, 
vol. 21, no. 2, 1957, pp. 207-212. (In 
Russian.) 

Interaction of a Shock Wave and the 
Boundary Layer in the Neighborhood of 
the Leading Edge of a Plane Plate for 
High Supersonic Velocities and With the 
Consideration of Radiation, by G. A. 


*The editors of Technical Literature 
Digest are making a systematic search for 
pertinent Russian articles, both in the orig- 
inal Russian and in translated form. For 
the balance of 1959, these will appear in a 
separate section with this heading, mainly 
to draw the readers’ attention to the new 
listing. After that, they will appear item 
by item under the proper subject headings. 


Kulonen, Vestnik Leningradskogo Uni- 
versiteta, Seriya Matematiki, Mekhaniki i 
Astronomii, no. 7 (2), 1958, pp. 172-188. 
(In Russian.) (Abstract in PB 141037T-1, 
Office of Technical Services, U. S. Dept. 
Commerce, $0.50. 

Form of Bodies for Minimum Wave 
Resistance Studied, by G. I. Kostychev, 
Izvestiia Vysshikh Uchebnykh Zabedeniy, 
Aviatsionnaya Tekhnika, no. 2, 1958, pp. 
9-15. (In Russian.) (Abstract in PB 
131891T-9, Office of Technical Services, 
U. 8. Dept. Commerce, p. 97, $2.75 per 
Issue. ) 

On Motions with Homogeneous De- 
formations in Magnetic Hydrodynamics, 
by A. G. Kulikovskiy, Doklady Akademii 
Nauk SSSR, vol. 120, no. 5, June 11, 
1958, pp. 984-986. (In Russian.) 


On Radiation Cooling of Air. I. Gen- 
eral Description of the Phenomenon and 
the Weak Cooling Wave, by Ya. B. 
Zel'dovich, A. S. Kompaneyets and Yu. P. 
Rayzer, Zhurnal Eksperimental’noy 1 
Teoreticheskoy Fiziki, vol. 34, no. 5, 1958, 
pp. 1278-1287. (In Russian.) 

The Kinetic Theory of Magnetohydro- 
dynamic Waves, by K. N. Stepanov, 
Zhurnal Eksper. i Teoret. Fiziki, vol. 34 
no. 5, 1958, pp. 1292-1301. (In Russian.) 


The Effect of Noise Fluctuations on the 
Accuracy of an Automatic Tracking Sys- 
tem With a First Order Astaticism and a 
Passband Controlled by the Input Signal, 
by Yu. M. Kazarinov, Yu. A. Kolomen- 
skiy and R. I. Smirnov, Jzvestiia Vysshikh 
Uchebnykh Zavedeniy-Priborostroyeniye 
no. 2, 1958, pp. 3-12. (In Russian.) 

Use of Doppler Effect to Determine the 


Orbital Parameters of Artificial Earth 
Satellites, by V. A. Kotel’nikov, V. M. 


Dubrovin, V. A. Morozov, O. N. Rzhiga 
and A. M. Shakhovskoy, Radiotekhnikai. 
Electronika, no. 7, July 1958, pp. 872-881 
(In Russian.) 


Silicon Photocells as Solar Energy Con- 
verters, by V. S. Vavilov, Atomnaya 
Energiya, vol. 4, June 1958, pp. 571-575. 
(In Russian.) 


The Effective Boundary Conditions in 
Neutron Diffusion Theory, by G. A. Bat’ 
and D. F. Zaretskiy, Atomnaya Energiya, 
vol. 4, June 1958, pp. 510-519. (In Rus- 
sian.) 

The Energy Spectrum of Neutrons from 
a Pulsed Source in a Heavy Moderator 
with Constant Path Length, by M. V. 
Kawarnovskiy, Atomnaya Energiya, vol. 
4, June 1958, pp. 539-546. (In Russian.) 


Shock Losses in a Divergent Nozzle for 
Supersonic Flows with Pressure Gradi- 
ents, by L. P. Volkova and M. Ia. Iudel- 
ovich, Akad. Nauk, SSSR, Izvestiia, Otdel. 
Tekh. Nauk, 4, April 1958, pp. 67-72. 
(In Russian. ) 

Stability of Shock Waves in Relativistic 
Hydrodynamics, by V. M. Kontorovich, 
Soviet Physics-JETP, vol. 34 (7), July 
1958, pp. 127-132. 

Energy Balance in Transition to Turbu- 
lent Flow, by B. A. Phidman, Akad. Nauk 
SSSR, Izvestiia, Otdelenie Tekh. Nauk, no. 
8, Aug. 1958, pp. 139-142. (In Russian. ) 

Plane Problems in Magnetohydrody- 
namics, by G. S. Solitsyn, Soviet Phys.- 
JETP, vol. 34 (7), no. 3, Sept. 1958, pp. 
473-476. 

On the Theory of Plasma Waves in a 
Degenerate Electron Liquid, by V. P. 
Silin, Soviet Phys-JETP, vol. 34 (7), no. 
3, Sept. 1958, pp. 538-539. 
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